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The  objective  of  this  dissertation  was  to  analyze  the  effect  of 
Inertlal  oscillators  In  energy  systems.  Traditionally  defined 
in  differential  equations  for  hardware,  generalized  systems  diagramming 
is  used  to  extend  the  knowledge  of  inertlal  mechanisms  to  environmental 
systems.  Equations  and  diagrams  are  presented  together,  first  in 
definition  and  review,  then  as  elements  that  are  incorporated  in 
ecosystem  patterns. 

Starting  with  the  simplest  configurations,  Inertlal  elements  were 
combined  Into  oscillators  and  these  into  ecosystem  configurations.  Each 
model  uses  computer  simulations  to  characterize  responce  in  the  time  and 
frequency  domains.  A  unit  oscillator  mechanism  was  included  as  a  clock 
mechanism  in  ecosystem  prototypes.  Potential  for  this  k'nd  of  a 
mechanism  was  considered  in  competition,  energy  selectivity  (the  ability 

xiii 


to  aquire  energy  input  in  a  selective  manner  through  use  of  an  'nertial 
oscillator),  diversity,  and  survival  of  systems.  For  each  model, 
diagrams,  equations,  time  simulations,  and  frequency  responses  are 
given . 

The  results  indicate  that  the  development  of  an  internal 
oscillation,  accomplished  for  this  research  by  an  autocata 1 yt ic  inertial 
oscillator,  increases  the  efficiency  of  energy  aquisition  within  a  range 
of  acceptable  frequencies,  increases  power  transfer  for  the  system  to 
other  systems  within  the  ecosystem,  increases  ecosystem  diversity,  and 
provides  the  system  with  a  greater  competitive  advantage  within  its 
niche.  The  most  competitive  system  is  produced  as  a  result  of  its 
ability  to  generate  multiple  internal  rhythms  that  it  then  employs  to 
better  time  changes  in  its  environment.  Competitive  coexistance  and 
exclusion  are  produced  by  combinations  of  system  configurations  emolying 
internal  inertial  oscillators. 
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INTRODUCTION 

This  research  is  a  direct  outgrowth  from  an  understanding  of  the 
operation  of  some  simple  oscillator  systems.  Knowledge  of  the  energy 
and  power  input  and  transfer  characteristics  of  simple  oscillators  lead 
to  the  idea  that  many  if  not  all  systems  employ  a  combination  of 
inertial  components  to  produce  internal  oscillations.  In  many 
configurations  that  use  inertial  elements,  inertial  oscillators  are 
formed.  Since  backforce  actions  are  known  in  the  behavior  of  biological 
and  social  systems  not  involving  hardware,  inertial  elements  and 
oscillators  may  exist  in  environmental  systems  or  may  be  available  for 
ecological  engineering  designs  of  environmental  systems  of  the  future. 
Therefore,  this  dissertation  is  a  study  of  the  role  of  inertial  elements 
and  inertial  oscillators  in  basic  configurations  of  ecosystems. 

Traditionally  defined  in  differential  equations  for  hardware, 
generalized  systems  diagramming  helps  the  extention  of  existing 
knowledge  on  inertial  mechanisms  to  environmental  systems.  In  the 
following  sections  equations  and  diagrams  are  presented  together,  first 
in  definitions  and  review,  then  as  elements  that  are  incorporated  in 
ecosystem  patterns.  The  simplest  inertial  element  and  its  equations  are 
shown  in  Figure  1.  An  important  characteristic  of  this  simple  inertial 
component  is  the  integral  dependency  of  the  flow  (F)  through  its 
pathway.  The  flow  through  this  inertial  component  is  dependent  upon  the 
value  of  its  inductance  (L)  and  the  integral  of  the  potential  energy  (S) 
forcing  the   flow  through  the  pathway.     The  energy   stored   in  the 
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Figure  1:  Fundamental  inertial  component.  The  figure  shows  the 
simplest  inertial  component.  Inertial  backforce  is 
generated  by  the  component  as  a  result  of  its  induced 
magnetic  field.  The  equations  describe  the  flow  through 
the  inertial  component  and  the  energy  stored  by  the 
component  in  its  magnetic  field. 
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magnetic  field  around  the  element  is  related  to  the  inductance  of  the 
component  (L)  and  the  square  of  the  flow  (F).  It  is  the  energy  stored 
in  the  magnetic  field  of  the  inductor  that  opposes  any  change  in  flow 
through  the  element.  The  opposition  to  changes  in  the  flow  through  the 
element,  by  the  fluctuating  magnetic  field,  produces  the  inertia! 
backforce  generated  by  the  element. 

In  this  study  system  configurations  with  inertia!  mechanisms  and 
oscillators  were  analyzed  and  related  to  the  design  of  environmental 
systems.  The  energy  circuit  model  simulations  and  plots  explain  and/or 
define  the  concept  of  energy-selective  power  utilization  and  transfer. 
First,  two  simple  oscillator  models  are  presented,  one  an  electrical 
res i stance- i nductor-capac i tor  circuit,  and  the  other  a  mechanical  mass- 
spring-damper  model.  Second,  an  autocata I yt i c  inertia!  oscillator  is 
shown.  The  autocata I yt ic  inertia!  oscillator  model  provides  a  means  of 
introducting  an  interna!  control  behavior  into  energy  circuit  systems 
which  in  turn  produce  energy-selective  behavior  within  models  of  energy 
circuit  producers  and  consumers.  Next,  the  simulations  for  models  of 
energy-selective  competition  are  shown.  This  series  of  models  analyzes 
the  advantages  or  disadvantages  of  energy-selective  behavior  for 
individual  producers  in  competitive  situations.  The  model  simulations 
also  aid  in  developing  an  insight  into  the  overall  system  benefits  from 
an  individual  system's  energy-selective  behavior.  The  last  series  of 
model  simulations  examine  the  effects  of  including  inertia!  oscillators 
into  a  salt  marsh  located  in  New  Zealand.  The  model  was  previously 
calibrated  with  biomass  data  and  simulated  seasonal  productivity 
(Campbell,  1984).  The  Kaitun  marsh  simulations  are  the  result  of 
modifications  to  the  original  Kaituna  marsh  model.   Energy   selective 


producers  are  substitutea  'nto  the  marsh  model  in  place  of  the  or-qinal 
non-energy-selective  producers  in  order  to  analyze  the  effects  of  energy 
selective  behavior  witnin  a  known  system.  Following  the  "nsert'on  of 
energy-selective  behavior  into  the  Kaituna  marsh  model,  the  effect  of 
competition  "s  shown  by  the  addition  of  competing  energy-se I ect ' ve 
producers  into  the  modified  marsh  model. 

Previous  Study  of  Electrical  and  Mechanical  Inertial  Mechanisms 

The  ability  to  act  with  a  behavior  that  produces  a  selec"*"'ve 
acquisition  of  input  energy  relies  upon  the  ability  of  the  system  to 
produce  internal  osc ' I  I  at 'ons .  The  simple  electrical  and  mechan'cal 
oscillators  provide  an  example  of  the  concept.  Electrical  circuits  "^hat 
exhibit  their  own  internal  inertial  oscillat'ons  are  called  bandpass 
oscillators  or  filters,  blocking  filters,  tuners,  receivers,  or  simply 
oscillators.  One  exarriple  of  a  mechanical  osc' Mating  system  's  an 
automobile's  suspension  system. 

Figure  2  shows  the  electrical  bandpass  oscillator,  first  w"th 
electrical  diagramming  and  second  with  energy  circuit  diagramming.  '!"he 
differential  equations  for  both  diagrams  are  presented.  The  inertial 
generated  oscillation  is  the  result  of  a  transfer  of  energy  internally 
between  the  two  components  of  the  system,  the  capacitor  (C)  and  the 
inductor  (L).  Interestingly  the  period  of  the  free-running  oscillation 
is  directly  related  to  the  two  components  (C)  and  (l)  and  not  to  the 
rate  of  flow  of  the  electrons. 

The  frequency  setting  coefficient  k3  is  a  constant  and  may  be 
calculated  usiny  two  different  methods.  The  first  method  is  the  most 
direct  and  is  solely  dependent  upon  the  inductance  (L)  and  capacitance 
(C)   of   the  circuit  (see  equation  6).    The  second  method   relies   upon 


Figure  2:  The  electrical  'nertial  oscillator.  The  figure  shows  both 
the  electrical  diagram  and  the  energy  circuit  for  this 
inertial  oscillator.  The  equations  describe  the  rate  of 
change  of  the  charge  (storage  Q)  In  the  oscillator.  The 
frequency  setting  coefficient  k3  Is  a  constant  for  flow 
through  patway  F3  and  Is  calculated  from  the  system 
component  values.  Coefficient  k3  may  also  be  calculated 
from  the  free-running  period  of  the  inertial  oscillator. 
Figure  (a)  is  the  electrical  diagram.  Figure  (b)  is  the 
energy  circuit  diagram. 
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being  able  to  discern  the  free-running  period  of  the  circuit.  Equation 
6  also  provides  the  calculations  for  period  method  to  determine  k3. 
Which  ever  method  is  employed  to  calculate  frequency  setting  the  end 
result  is  that  k3  is  a  constant  and  not  a  van, able  flow  coefficient. 

Figure  3  shows  the  mechanical  mass-spring-damper  inertial 
oscillator.  The  figure  presents  the  mechanical  diagram  and  the  energy 
circuit  diagram  with  their  accompanying  differential  equations.  The 
mechanical  oscillator  is  much  like  the  electrical  oscillator  and  with 
the  exception  of  the  system  components  and  units  has  differential 
equations  that  are  much  alike. 

The  frequency  setting  coefficient  (k3)  for  the  mechanical 
oscillator  may  also  be  calculated  using  two  methods,  and  is  a  constant 
and  not  a  van  I ab I e . 

The  previous  figures  clearly  show  that  the  basic  design  for 
development  of  an  internal  oscillation  in  simple  inertial  oscillating 
systems  is  the  combination  of  two  components  that  act  'n  a  reciprocal 
manner.  These  components,  one  potential  and  the  other  kinetic  combine 
to  transfer  energy  between  themselves  and  as  a  result  provide  their 
system  with  an  internal  inertial  oscillation.  Furthermore,  the  simple 
oscillator  systems,  such  as  the  electrical  and  mechanical  oscillators, 
maximize  their  power  transfer  in  the  most  efficient  manner  when  the 
period  of  their  internal  oscillations  exactly  match  those  of  their 
external  supply.  When  the  period  of  the  system's  internal  oscillation 
is  longer  or  shorter  than  that  of  its  energy  source,  the  power  transfer 
and  power  input  to  the  system  decreases. 


Figure  3:  The  mechanical  inertial  oscillator.  The  figure  shows  the 
mechanical  diagram  and  the  energy  circuit  model  for  the 
mechanical  inertial  oscillator.  The  equations  describe 
the  rate  of  change  of  the  mass  position  (x).  The 
coefficient  k3  is  a  constant  for  setting  the  period  of  the 
system  oscillation  and  is  calculated  from  the  system 
components  or  from  knowing  the  system  free-running  period. 
Figure  (a)  is  the  mechanical  diagram  and  (b)  is  the  energy 
c  ircu  it  d  i  agram  . 
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Previous  Studies  of  Inertial  Mechanisms  !n  Environmental  Systems 

Alexander  1973,  studied  the  cycles  of  order  and  d^sorGer  producing 
a  quantitative  theory  based  upon  energy  analysis  for  global  cycles, 
surges,  and  disasters.  Storns,  earthquakes,  urban  f'res,  flocds, 
volcanic  eruptions,  wars,  and  forest  fires  were  rncdelea.  Order  and 
disorder  were  recognized  as  component  parts  of  a  cycle  that  emoloyed 
disorder  to  recycle  ordered  material  to  stimulate  new  growth  .vith'n  the 
system.  A  global  model  was  employea  by  Alexanaer  (1973)  to  'nd'cafo  the 
effects  of  the  cascading  of  catastophic  events  within  a  web  system. 

One  particular  model  developed  by  Alexander  (1973),  the  Exponenf'al 
Surge  Model,  was  employed  to  simulate  actual  disaster  processes.  ^he 
model,  produced  a  pulsed  self-started  disorder,  thereoy  recycling  orged 
to  disorder  and  producing  a  new  growth  period. 

Campbell  1984,  investigated  the  role  of  ecosystem  structure  'n 
absorbing  or  filtering  external  energy  inputs  to  the  ecosystem.  ^he 
organization  of  system  filter  properties  was  shown  to  prov'.:e  the 
ecosystem  with  the  ability  to  adapt  to  a  variety  of  "nput  frequencies. 
Turnover  time  for  the  structural  storages  w'th'n  a  system  when  matched 
to  the  input  energy  frequencies  usually  determined  the  power  capture  of 
the  system.  Furthermore,  Campbell  discovered  that  the  h?gh  frequency 
cutoff  of  the  individual  components  in  an  ecosystem  and  their  particular 
resonant  peaks  were  determined  by  the  turnover  time  of  the'r  structure. 
However,  these  oscillators  do  not  involve  inertial  elements. 

In  his  book  Systems  Ecology:  An  Introduction  Odum  1983,  sresents  a 
number  of  models  that  utilize  potential  and  kinetic  components.  Odum 
also  developed  the  energy  symbol  for  inertial  response — a  bearded  tank 
(Figure  5)  (Odum,  1967).  The  energy  circuit  model  shown  ^n  Figure  4  was 


Figure  4:  Odum's  diagram  of  the  Van  Der  Pol  oscillator.  The  model 
utilizes  a  kinetic  storage  as  a  backforce  to  produce  the 
system  oscillation.   Odum  1983,  pages  176-77. 
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drawn  by  Odum  as  a  translation  of  the  Van  Der  Pol  oscillator.  An 
oscillation  is  developed  in  the  system  because  of  the  presence  of  an 
opposing  interial  backforce.  The  system  does  not  contain  a  potential 
storage,  only  a  storage  of  velocity  energy,  but  an  oscillation  is 
produced  by  the  interaction  of  an  inertial  backforce,  a  frict"onal 
force,  and  a  driving  force. 

Odum,  1983  further  describes  the  control  effect  of  pulses  on 
timers.  Using  a  model  of  a  mechanical  chain  which  ends  in  a  timing 
pulse,  he  shows  that  timing  devices  require  a  reference  pulse.  The 
pulse  mechanism  has  a  kinetic  property  in  the  oscillating  wheel. 
Feedback  of  this  timing  pulse  provides  control  to  the  slower  cogwheels 
in  the  mechanical  chain.  The  conclusion  states,  "If  most  systems 
develop  chains  and  develop  high  quality  with  pulsing  delivery  of  the 
feedbacks,  a  clock  is  automatically  provided  so  that  the  systems  can 
utilize  timing  functions"  (Odum,  1983,  p.  284). 

Internal  Timing  vs  External  Timing  of  Inertia!  Oscillators 

Whether  the  autocatal yt ic  oscillator  is  environmentally  dependent, 
environmentally  sensitive,  or  completely  autonomous  has  an  effect  on  the 
way  the  rhythm  is  produced.  In  fact,  the  major  point  of  contention  about 
biological  clocks  centers  around  their  utilization  of  or  dependency  upon 
persistent  environmental  signals.  The  two  points  are  best  explained  by 
the  use  of  two  different  types  of  man-made  clocks.  The  internal 
autonomous  biological  clock  is  likened  to  an  alarm  clock.  The  alarm 
clock  contains  all  the  components  necessary  to  produce  internal  timing. 
The  operation  of  the  clock  is  completely  autonomous,  it  keeps  time  by 
its  own  mechanism.  The  external  clock  is  best  described  by  the  sun 
dial.   The  sun   dial  has  an  internal  mechanism  but  the  timing   is   the 
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result  of  an  external  signal.    Therefore  the  sun  dial  acts  only   as   a 
receptor . 

In  this  study  the  autocata I yti c  inertial  oscillator  is  an 
autonomous  internal  oscillator  with  a  free-running  rhythm.  The 
oscillator  produces  rhythmic  activity  that  will  persist  for  some  time 
without  any  environmental  cues.  The  development  of  energy-selective 
behavior  may  require  that  a  system  possess  the  ability  to  sense  changes 
in  its  environment.  With  this  in  mind,  the  autocata I yt ic  inertial 
oscillator  senses  its  environment  and  its  activity  responds  to  changes 
in  that  environment.  The  importance  of  possessing  an  internal 
oscillator  that  has  its  own  free-running  rhythm  and  is  capable  of 
sensing  and  changing  with  its  environment  is  studied  in  this  research. 
The  ability  to  aquire  energy  in  a  selective  manner  may  depend  (amoung 
other  characteristics)  upon  a  system's  ability  to  generate  an  autonomous 
internal  rhythm  which  senses  changes  in  its  environment  and  responding 
to  those  changes. 

Definition  of  Energy  Selectivity 

The  "principle  of  energy  selectivity"  could  be  stated  as  follows. 
When  a  system  is  subjected  to  a  cyclic  energy  input  from  its  environment 
the  system  will  best  insure  its  ability  to  obtain  energy  by  developing 
an  autonomous  Internal  rhythmic  behavioral  guide  that  is  capable  of 
sensing  its  environment.  The  closeness  of  the  system's  internal  rhythm 
to  that  of  its  primary  environmental  cue  determines  the  system's  energy 
acceptance  ability,  and  combined  with  other  environmental,  internal  and 
system  factors  its  ultimate  survivability. 

Energy-selectivity  may  play  an  important  role  in  the  ability  of  a 
system   to  maximize  its  input  energy  and  power.   Since  the  ability  to 
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aqu"re  energy  and  power  has  a  direct  bearing  on  a  system's  ablity  to 
transfer  power  to  the  next  larger  system  energy-selective  behavior  would 
also  play  an  important  role  in  the  power  transfer  amoung  systems.  If  a 
system  can  act  in  an  energy-selective  manner  that  activity  may  help  the 
system  define  its  niche  within  the  next  larger  system.  One  hypothesis 
for  this  research  is  that  energy-selective  behavior  provides  a  system 
with  a  means  of  maximizing  its  input  power,  and  as  a  result  a  means  of 
maximizing  its  power  transfer  to  the  next  component  or  larger  system. 
Within  biological  systems,  the  ability  to  act  in  an  energy  selective 
manner  provides  the  system  with  a  means  for  al locating  input  energy 
among  the  various  subsystems.  The  ability  to  allocate  input  energy  may 
allow  the  formation  of  a  greater  variety  of  niches,  thereby  increases 
system  diversity  and  maximizing  system  power. 

Energy-selectivity  behavior  is  an  ability  possessed  by  a  system  as 
the  result  of  its  generation  of  internal  rhythmic  activity  which  allows 
the  system  to  define  a  range  of  energies  it  will  accept,  all  other 
energies  are  rejected. 


METHODS 

Starting   with  the  simplest  configurations,   inertial  elements  were 

combined    into   oscillators   and   these    into   various   ecosystem 

configurations.    For   each   model   computer  simulations   were   used   to 

characterize  response   in   the  time  and   frequency   domains.   A   unit 

oscillator  mechanism  was   included  as  a  clock  mechan'sm   in   ecosystem 

prototypes.    Potential   for  this  kind  of  a  mechanism  was  considered   in 

competition,   energy   selectivity,   diversity,   and  survival  of  systems. 

For   each  model,   diagrams,   equations,   time  simulations  and   frequency 

responses  are  given. 

The   following   models  were  developed  or   utilized   to   investigate 

energy  selectivity  in  biological  systems: 

1.  a  model  of  an  el ectr ical /mechan ica I  oscillator; 

2.  a  model  for  an  autocata I yt ic  inertial  oscillator; 

3.  a  model  for  an  energy-selective  producer; 

4.  models     for    competition    between    energy-selective 
producers; 

5.  a  model  of  the  Kaituna  Salt  Marsh  - 

a)  the  original  model  was  simulated  to  check  the  normal 
model  operational  characteristics, 

b)  the  model   was   then  modified  to   account   for   the 
diurnal  solar  input, 

c)  the  model   was  modified  to   Incorporate  an   energy 
selective  producer, 

d)  the   model   was   modified   to   Investigate   energy 
selective  competition. 
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Overall  Research  Procedure 

Five  general  sections  were  employed  to  ascertain  the  effects  of 
energy-selectivity  on  the  power  transfer  of  systems  with  cyclical 
external  forcing  functions. 

Phase  I 


The  first  phase  of  this  research  was  devoted  to  understanding  the 
conceptual  similarities  in  simple  oscillator  systems.  Special  emphasis 
was  placed  toward  understanding  what  characteristics  of  these  simple 
oscillators  could  be  used  to  simulate  the  system  rhythmic  activity.  Two 
e 1 ectr i cal /mechan ica 1  models  were  tested.  The  electrical  and 
mechanical  models  were  simulated  repeatedly  (with  any  necessary 
modifications  to  the  models)  until  the  simulation  results  compared  with 
known  system  responses.  The  two  original  choices  were  (a)  a  series 
approach   to  the  production  of  an  electrical   oscillation   and   (b)   a 

parallel  tank  circuit  or  bandpass  approach  to  the  production  of  an 
electrical  oscillation.  Ultimately,  the  tank  circuit  oscillator 
provided  the  the  best  opportun i ty  for  expansion  into  a  model  for  the 
simulation  of  an  autocata 1 yt i c  inertial  rhythm. 

Phase  I  1 


Phase  two  was  used  to  produce  a  phase  sensing,  time  interpreting 
autocatal yt ic  inertial  oscillator  (a  system  rhythm  generator).  The 
development  of  such  a  unit  was  deemed  necessary  because  the  simple 
addition  of  a  linear  oscillator  to  a  producer  or  consumer  unit  within  a 
model  would  provide  only  a  rhythmic  linkage  and  not  phase  adjustment  or 
age-sequence  timimg.    Many  models   for   the  autocatal ytic   inertial 
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oscillator  were  tested.  However,  the  final  oscillator  model  (presented 
in  the  results  section)  is  capable  of  sensing  source  phase  shifts  and 
adjusting  to  those  sh'fts.  The  autocata I yt i c  inertia!  osc'llator  will 
allow  for  the  adjustment  of  its  internal  rhythm  and  oscillator  activity 
range  (the  time  the  oscillator  will  function  'n  constant  conditions); 
and  can  produce  a  variety  of  source  sensing,  internal,  and  age  related 
timimg  functions. 

Phase  I  I  I 


The  third  phase  of  this  analysis  was  used  to  s'mulate  variations 
of  a  model  for  a  simple  energy-selective  producer.  The  model 
incorporated  a  autocata I ytic  inertial  oscillator  into  a  producer  un't 
thus  allowing  the  unit  to  develop  an  energy  acceptance  niche  within 
which  it  could  or  would  accept  energy  from  one  or  more  cyclic  energy 
sources.  The  unit  and  load  power  for  the  producer  was  calculated  during 
the  simulation.    The  model  design  was  based  on   the  assumption  that  the 

producer's   load  (i.e.   its  output  flow)  was   the  unit's  connection  to  a 
more  complicated  system. 

Phase  IV 

Phase  four  was  employed  to  simulate  a  number  of  variations  for  a 
competition  model.  These  models  simulated  the  competition  between  two 
energy-selective  producers  for  the  use  of  an  external  cyclic  energy 
source.  The  source  rhythm  and  the  internal  behavioral  rhythms  of  the 
competitors  were  varied  in  successive  computer  simulations  to  analyze 
the  effect  of  energy-selectivity  on  competition  and  input  energy  niche 
development.  Source  energy  supply  was  also  designed  to  investigate  the 
effects  of  various  combinations  of  source  fundamental  rhythms,   i.e.  the 
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addition   or   multiplication  of  two  source  rhythms  to  form  one   variable 
energy  supply. 

Phase  V 

Phase  five  was  employed  to  simulate  the  existing  Kaituna  salt  marsh 
mini-model  developed  by  Odum,  Knox,  and  Campbell,  1983.  The  model  was 
chosen  because  it  provided  a  proven  mini-model  with  cyclic  external 
forcing  functions.  The  model  had  already  been  evaluated  and  was  of 
reasonable  proportions  to  enable  the  testing  of  energy-selectivity 
without  excessive  complexity.  Furthermore,  the  model  provided  a 
reasonable  framework  for  the  analysis  of  energy-selective  behavior 
within  a  webbed  system. 

Energetics  Diagramming 

The  methodology   employed   to   design   or   construct   the   models 

presented   in   this   dissertation   was  developed  by  Alexander   et  al., 

(1980).   The  procedure  is  as  follows: 

Step  1 .   Construct  a  system  boundary. 

Step  2.    Identify    the   energy   flows   across    the    system 
boundary. 

Step  3.    Organize  the  system's  components, 

Step  4.   Connect   the   system  components   to   the  external 
sources, 

Step  5.   Diagram  the  subsystems. 

Step  6.    Evaluate  of   assign   values   to   the   energy   model 
d  iagram, 

Step  7.    Translate  the  energy  diagram  into  its   corresponding 
differential  equations, 

Step  8.    Simulate  the  model,  and 

Step  9.   Validate  the  model's  response  or  check  response   to 
see  if  it  meets  expectations. 
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Identification  of  the  System  Boundary 

Systems  that  are  open  exist  In  a  hierarchy.  To  differentiate 
between  the  storages,  flows,  and  sources  for  one  system  from  that  of 
another — a  boundary  must  be  defined.  Boundary  definition  for  the  models 
employed  in  this  research  was  based  upon  restricting  complexity  in  model 
design.  The  concept  of  power  utilization  and  transfer  is  complex  enough 
without  complicating  the  simulations  with  large  cumbersome  models.  The 
mini-model  approach  to  analysis  requires  the  use  of  a  simple  model  to 
illustrate  a  particular  concept. 

I  dent! f [cation  of  Energy  Flows  Across  the-Boundary 

Alexander  et  al.,  (1980)  have  developed  a  method  for  positioning 
energy  sources  around  a  system  boundary.  First,  all  pertinent  external 
sources,  with  flows  across  the  boundary,  are  identified.  The  source 
with  the  lowest  energy  quality  is  placed  in  the  lower  left  corner  of  the 
diagram  outside  of  the  boundary.  Next,  the  remaining  sources  are  then 
placed  around  the  boundary  in  a  clockwise  direction  by  increasing  energy 
qua  I ity . 

Organization  of  the  System  Components 

Energy  language,  as  developed  by  H.T.  Odum  (1976),  is  made  up  of  a 
set  of  symbols,  shown  in  Figure  5.  These  symbols  are  employed  to 
construct  energy  models,  and  are  representative  of  component  functions 
which  repeatedly  occur  In  all  systems. 

Figure  5a  is  a  tank  or  storage  and  Is  used  to  identify  system 
state  variables.  All  storages  have  a  particular  turnover  time  and 
capacity  to  store  energy.  The  turnover  time  of  a  tank  varies  depending 
upon  the  particular  pathway  flow  controlling  the  turnover.   Examples  of 
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turnover  time  include  input  turnover  and  output  turnover.  Turnover 
time  can  be  found  by  dividing  the  tank  storage  by  an  input  or  output 
pathway  flow.  A  tank  with  a  storage  of  1000  Cal  and  an  output  pathway 
flow  of  50  Cal/hr  would  have  a  turnover  time  of  20  hours.  The  capacity 
of  a  tank  is  fixed  and  can  best  be  described  by  a  city  water  tank.  The 
capacity  for  any  tank  's  determined  by  its  design,  and  for  discussion 
let  the  capacity  of  a  hypothetical  city  water  tank  be  5000  gallons. 
The  turnover  of  this  tank  if  it  were  half  full  (2500  gallons)  and 
draining  at  250  gal/hr  would  aga'n  be  10  hours,  but  the  capacity  would 
remain  5000  gallons. 

Figure  5b  is  the  symbol  for  an  external  source.  This  type  of  a 
source  provides  a  constant  potential  which  drives  the  flow  along  its 
associated  pathway.  The  greater  the  system  requirement  the  larger  the 
f  low. 

Figure  5c  is  a  energy  sink  and  represents  total  energy  degradation 
of  a  system.  The  energy  entering  the  system  sink  is  no  longer  usable  by 
the  system. 

Figure  5d  is  a  flow-limited  source.  A  flow-limited  source 
provides  flow  energy  only.  One  example  of  a  flow  limited  source  is  the 
solar  energy  provided  to  the  earth.  While  this  energy  may  increase  or 
decrease  at  any  given  moment,  the  total  energy  supplied  can  never  exceed 
the  value  of  J,  and  the  amount  of  energy  available  for  use  at  any  given 
time  is  Jr.  The  only  option  available  to  a  system  dependent  upon  this 
type  of  source  (if  it  wants  to  get  more  energy)  is  to  become  more 
adapted  at  capturing  the  energy  as  it  passes.  This  source  differs  from 
the  previously  described  source  because  it  requires  the  system  to  work 
harder  or   to  become  more  efficient  in  order  to  obtain   energy.   The 
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source  in  Figure  5b  is  capable  of  supplying  as  large  a  flow  as  the 
system  requires,  and  as  a  result  can  lead  to  explosive  growth  if  not 
limited  by  some  other  means. 

Figure  5e  is  a  bearded  tank  and  stores  energy  as  a  result  of  an 
acceleration,  i.e.  kinetic  energy.  The  energy  stored  in  this  type  of  a 
tank  is  due  to  a  change  In  pathway  flow.  The  beard  on  the  tank  was 
derived  by  Odum  as  a  result  of  the  combination  of  two  symbols — one  from 
energetics  and  the  other  from  electronics.  The  two  symbols  were  the 
energy  tank  in  Figure  5a  and  the  electrical  symbol  for  an  inductor  (see 
Figure  la,  in  the  previous  section). 

Figure  5f  is  a  transformation  or  work  gate.  The  interaction  of 
two  or  more  state  variables  produce  a  flow  of  energy  (k*A*B).  The 
interaction  process  is  defined  by  the  mathematical  symbol  or  equation 
placed  inside  the  symbol. 

Figures  5g  is  a  energy  flow  sensor.  The  flow  of  energy  In  pathway 
Fl  Is  sensed  and  utilized  to  produce  or  generate  some  other  funct'on 
within  the  system.  Photoper lod i sm  relies  on  a  sensor  to  determine  the 
amount  of  light  available. 

Figure  5h  Is  a  group  symbol  for  a  production  unit.  The  producer 
mixes  low  quality  energy  with  high  quality  energy  to  produce  an 
Intermediate  quality  energy. 

Figure  51  Is  another  group  symbol,  the  consumer.  The  consumer 
symbol  defines  a  unit  that  Is  self-contained  and  feeds  control  energy 
back  to  other  parts  of  the  system. 

System  components  are  arranged  within  an  energy  circuit  diagram  by 
energy  quality.  Inside  the  system  boundary  the  lowest  quality  component 
Is   located   in  the  upper  left  corner  of  the  diagram.   Components   with 


Figure  5:  Energy  circuit  symbols.  Symbol  (a)  storage  tank  or  state 
variable,  (b)  constant  potential  source,  (c)  energy  sink, 
(d)  flow  limited  source,  (e)  bearded  tank  which  stores 
kinetic  energy  and  is  capable  of  producing  intertial 
backforce,  (f)  work  gate  or  energy  transformation,  (g) 
flow  sensor,  (h)  producer  group  symbol,  (I)  consumer  group 
symbol  (odum,  1976  and  1983). 
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higher   quality  are  then  located  in  a  descending  order  toward  the   lower 
left  corner  of  the  diagram. 

Connection  of  System  Components  With  External  Sources 

System  components  are  connected  to  an  outside  source  when  it  is 
determined  that  there  exists  a  flow  of  energy  from  the  source  to  the 
component.  Solar  energy  to  plants  (producers)  would  constitute  a 
pathway  connection  from  a  source  of  external  solar  energy  to  a  producer 
within  the  system  boundary. 

Diagramming  the  Subsystems 

Each  group  symbol  represents  a  subsystem  within  the  model,  and  "s 
constructed  in  the  same  manner  as  the  overall  model.  Each  subsystem  is 
developed  by  combining  variations  of  the  remaining  symbols  to  produce 
subsystem  functions.  Mathematical  notations  are  placed  inside  of  the 
Interaction  symbols  used  in  the  subsystems  and  in  the  overall  system  to 
Identify  their  specific  operational  characteristics. 

Location  of  energy  symbols  within  the  subsystems  is  determined  by 
their  respective  energy  quality.  The  lowest  quality  storage  is  located 
in  the  upper  left  corner  within  subsystem,  while  the  highest  quality 
storage  is  located  in  the  lower  right  corner  within  the  group  symbol. 
Interaction  or  transformation  symbols  are  generally  located  to  the  right 
of  the  tank  or  tanks  that  It  supplies. 

Evaluation  of  the  Energy  Circuit  Model 

The  first  and  second  laws  of  thermodynamics,  conservation  and 
degradation  of  energy,  are  employed  to  evaluate  energy  circuit  models. 
Since  energy  can  be  neither  created  nor  destroyed,  energy  interactions 
must   balance.    Therefore,   the   energy  flows  entering   an   Interaction 
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symbol  must  be  equal  to  the  energy  flows  leaving  tha+  symbol.  The 
second  law  requires  that  all  storages  and  interaction  have  an  energy 
pathway  to  the  system  sink.  These  flows  represent  storage  depreciation 
and  work  process  losses. 

Translation  of  the  Energy  Circuit  Model  into  its  Corresponding 
Differential  Equations 

Energy   circuit  models  can  be  expressed  in   nonlinear   differential 

equations.    For   example,   the  model  in  Figure  6   contains  a  producer, 

with   an  autocatal ytic  inertial  oscillator  controlling  the  rate  of   the 

producer's  biomass  production.   The  model  may  be  translated  as  follows: 

Jr  -  J/{1+kO*l ) 

dP/dt  -  k5*Jr*/P/-k6»Jr/P/-k4*P-k3*  /P 

dB/dt  -  kl»Jr*l-k2*B 

Variables  (P)  (K)  and  (3)  are  called  "state  variables"  because  they 
determine  the  state  of  the  system.  Constants  kl  ,k2  ,k3 ,  . .  .k7  do  not 
change  during  simulation,  and  are  proportionality  constants  of  flow 
which  determine  the  percent  flow  in  their  respective  pathways.  Flows 
entering  a  state  variable  are  positive  (kl*Jr*l)  while  flows  leaving  a 
state  variable  are  negative  (-k2*B).  The  interaction  in  Figure  6  is 
multiplicative  and  is  indicated  as  such  by  the  "X"  located  in  the  center 
of  the  symbol . 

Balancing  the  flows  around  an  interaction,  as  required  by  the  first 
law  of  thermodynamics,  is  accomplished  by  ensuring  that  each  pathway 
equation  contains  the  same  number  and  type  of  storages  and/or  external 
sources.   For   example,   each   pathway  in  Figure  6   that  touches   the 


Figure  6:  Model  for  Translation  of  diagram  into  equations.  The 
energy  circuit  model  is  translated  into  the  following 
equations: 


Jr  =  J/(1+kO*l ) 

dP/dt  -  kl*Jr*/P/-k6»Jr»/P/-k4*P-k3»/P 
dB/dt  =  kl*Jr*l-k2*B 
NOTE:  The  "/P/"  is  the  absolute  value  of  storage  P 


29 


1    r        2 

dP/dt=kcJ   /P/-k,J    "D'-k   D-k,'? 
5    r  6    r'     '       4        3 


30 

interaction   symbol   has   storages   "Jr"  and   "I"   in   their   respective 
equations,  (see  previous  differential  equations). 

Simulation  of  the  Model 


Energetics  models  may  be  simulated  with  an  analog  or  digital 
computer.  Although  the  analog  computer  is  capable  of  ' nstantaneous 
simulation,  the  necessary  scaling  and  physical  setup  requ'red  for  a 
model  requires  a  lengthy  effort.  Therefore,  simulation  of  the  models 
for  this  research  was  accomplished  with  a  digital  computer  using  the 
Continuous  System  Modelling  Program  (CSMP).  The  models  were  simulated 
on  the  Northeast  Regional  Data  Center's  main  frame  computers  located  'n 
the  Space  Sciences  building  on  the  University  of  Florida  campus, 
Gai  nesv  i I  I e. 

The  CSMP  language  employs  "call  functions"  or  "routines"  for 
integration  of  differential  equations.  The  program  allows  the  selection 
of  the  integration  time  interval  (DT)  and  length  of  tiine  for  the 
simulation.  Example  equations  for  the  rate  of  change  of  the  system 
state  variables  have  already  been  presented.  The  integration  of  that 
rate  of  change  produces  the  levels  for  the  tanks  in  the  model  at  any 
time  during  the  simulation. 

Va! Idatlon  of  the  Model 


Validation  of  complex  systems  models  constructed  from  sets  of 
nonlinear  differential  equations  cannot  be  accomplished  by  statistical 
methods.  Therefore,  other  methods  of  validation  must  be  employed.  Most 
often  models  are  constructed  from  nonlinear  equations  that  attempt  to 
reproduce  the  characteristics  of  the  system.  This  approach  is  not 
without   its   faults.   Namely,   a  variety  of  models  can  be  produced   to 
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reconstruct  a  g'ven  set  of  data  or  conditions.  Construction  of  models 
with  no  underlying  hypothesis  or  theory  often  produces  models  which  only 
attempt  to  replicate  the  data  and  as  a  result  significantly  reduce  the 
credibility  of  the  model. 

To  avoid  that  problem  the  models  constructed  for  th's  analysis  have 
supporting  hypotheses  as  the  basis  for  their  operation  and  design.  The 
model  for  the  autocata I yt i c  inertial  osc'l I ator  is  based  upon  the  theory 
that  system  rhythmic  activity  is  a  function  of  a  nonlinear  oscillators 
constructed  of  potential  and  kinetic  components.  Testing  of  th 
oscillator  model  was  accomplished  by  analysis  of  existing  data  about  th 
operational  phenomena  for  system  rhythmic  activity.  The  model  of  th 
autocata I yt ic  inertial  oscillator  was  coastructed  based  upon  theory  and 
not  to  rep  I  icate  a  set  or  sets  of  known  data. 

Models  for  the  biological  producers  were  constructed  based  upon  the 
established  energy  theories  and  modelling  design  algorithms  developed  by 
H.  T.  Odum,  1976  and  1983.  The  analysis  of  energy-selective  behavior 
was  then  tested  via  the  simulation  of  these  hypothetical  producers. 
Finally  energy-selective  subsystems  were  added  to  an  operational, 
evaluated,  and  previously  simulated  models  of  a  salt  marsh  ecosystem. 
The  ultimate  test  for  these  subsystems  was  a  comparison  of  the  model 
established  operational  characteristics  vs.  its  new  operational 
characteristics  with  the  addition  of  energy-selective  behavior. 
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RESULTS 

Results  and  analyses  are  presented  in  four  sections:  concepts, 
models  for  biological  units,  models  for  competition  between  biological 
units,  and  models  for  a  the  Kaituna  salt  marsh,  located  in  New  Zealand. 
Each  individual  section  presents  the  model (s),  accompanying  simulation 
plots,  and  applicable  discussion  or  explanation  of  the  results. 

Concepts 

The  general  concept  for  energy-selectivity,  as  a  result  of  an 
i ntr i ns ic-autocatal ytic-rhythm ic-f unction,  results  from  the  extention 
of  concepts  from  the  electrical  and  mechanical  inertial  oscillators. 
Once  the  genera!  principles  for  these  simple  systems  had  been 
ascertained  the  development  of  a  model  for  the  autocata I ytic  inert'al 
oscillator  was  undertaken. 

Electrical  and  Mechanical  Oscillator  Models 

System  components  for  the  electrical  tank  circuit  oscillator  are 
shown  in  Figure  7a.  The  inertial  oscillator  contains  a  capacitor, 
inductor,  and  a  source  and  load  resistance,  Rs  and  Rl  respectively.  The 
corresponding  energy  circuit  diagram  is  shown  in  Figure  7b.  While  the 
electrical  diagram  indicates  component  connectivity,  the  system  storage 
information  for  charge  or  potential  energy  in  the  capacitor  (C) , 
magnetic  flux  or  kinetic  energy  around  the  inductor  (L),  and  circuit 
total   energy  are  not  shown.   The  inadequacy  of  this  type  of   diagram 
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Figure  7:  The  electrical  tank  circuit  oscillator.  Diagram  3)  an 
electrical  diagram  of  this  bandpass  type  of  electrical 
oscillator,  b)  an  energy  circuit  diagram  for  tne  tank 
oscillator.  Note  the  indication  of  the  storages  of  energy 
in  the  energy  language  diagram,  an  important  aspect  of 
this  circuit  which  is  not  shown  in  the  electrical  diagram. 
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dQ/dt=l/R^*(S-(Q/C))-(l/R^*C)*Q-l/(L*C).-Q 
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i2=(VL*C).-Q 
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caoacitor  {C) 
C=caDaci tance  (Farads) 
L=induc:ance  (Henries) 
fQ=intrinsic  frequency  of  the 
oscillator 


/R  *(s-(g/c))--, 


1/(R^*C)*Q 


dO/dt=l/R^*(S-(Q/C))-l/(R^*C)*0-l/(L*C).-q 
dO/dt=k^*(S-(Q/C))-k2*Q-k3*.-Q 

Energy  in  Q=Q^/(2*C)   Energy  in  K=(  L/2)*(  1/ (  L*C)*.-Q)^  =  q'^/ (3*L*C) 
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arises  from  the  lack  of  descr i pti veness ,  i.e.  the  storage  information  is 
left  to  the  reader's  imagination.  On  the  other  hand,  the  energy  circuit 
diagram  clearly  indicates  component  connectivity  as  well  as  the  storage 
of  energies  within  the  system.  Charge  potential  is  stored  in  tank  (Q) 
and  flux  kinetic  energy  is  stored  in  tank  (K).  Total  circuit  energy  is 
represented  by  tank  (E).  The  beard  on  tank  (K)  indicates  a  back  pressure 
effect  on  the  flow  of  electrons  within  its  supporting  pathway  r3  and 
results  from  magnetic  field  interaction  with  pathway  electron  flow. 

Table  1  shows  component  values  and  flows,  storage  initial 
conditions,  description,  and  equations  for  the  electrical  tank  circuit 
oscillator.  Table  1  also  presents  the  equations  for  pathway 
coefficients,  and  coefficient  values  for  the  same  oscillator.  Since 
the  model  is  designed  to  calculate  Its  own  pathway  coefficients  the 
equations  shown  in  Table  1  are  the  equations  within  the  model  that 
calculate  the  flow  coefficients.  The  initial  conditions  for  tanks  (Q), 
(K),  and  (E)  are  assumed  to  be  zero.  However,  If  tank  (Q)  has  an 
Initial  charge  then  the  initial  conditions  for  tanks  (K)  and  (E)  are 
calculated  by  the  model.  The  two  equations  necessary  to  accomplish  this 
are  shown  In  the  table. 

One  interesting  characteristic  of  these  models  's  the  integral 
dependency  of  flow  F3  on  tank  (Q).  This  phenomenon  results  from  the 
fact  that  tank  (Q)  provides  the  potential  force  that  drives  current 
through  pathway  F3,  while  tank  (K)  opposes  any  change  in  the  flow  of 
current  through  pathway  F3 .  The  result  produces  an  integral  dependent 
flow  In  pathway  3.  The  build  up  of  magnetic  flux  around  the  inductor 
acts  as  a  storage  of  kinetic  energy,  and  the  transfer  of  energy  between 
tanks  (.Q)    and  (K)  produces  the  system's  Internal  oscillation. 
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Figure  3  's  a  plot  of  the  electrical  inertial  osc'llator's 
storages  vs  time.  The  system  is  driven  by  a  source  potential  which  also 
varies  with  time.  The  important  conceptual  characteristics  shown  by  the 
p I ot  include: 

1)  the  oscillating  nature  of  (Q)  and  (K)  with  time, 

2)  the  90  degree  or  -  /4  radian  lag  of  (K)  to  (0), 

3)  the  effect  of  heavy  damping  on  the  system  storages  as  the 
source  supply  is  interrupted,  and 

4)  the   constant  amplitude  of  (Q)  and  (K)  as  long  as   source 
supp I y  is  ava i  I ab I e. 

Figure   9a    shows   the   mechanical   inertial   oscillator.    Three 

components  make  up  the  system;   they  are  the  mass,   spring,  and  dampers. 

Within   this   system   potential  energy  Is  stored  by   the   spring,   whil 

kinetic  energy  is  stored  by  the  position  of  the  mass.   The  correspondin 

energy  circuit  diagram  is  shown  in  Figure  9b.     As  with  the   electrical 

inertial   oscillator,   the   energy   diagram  more  clearly   indicates   the 

storages   for  potential  energy  (X)  as  shown  by  the  displacement   of   tne 

system's  mass,  kinetic  energy  (K),  and  total  energy  (E).   Table  2   shows 

the  storage  and  flow  initial  cond'tions,  component  values,  a  descriptive 

comment,    and   the   equations   for   the  mass-spr'ng   oscillator.     An 

intriguing   characteristic   of  these  two  inertial  oscillator  systems   is 

seen   by  a  comparison  of  their  pathway  equations.    Mote   the   identical 

equations   for   system   pathway   coefficients,   with   the   exception   of 

coefficient   units   the  system  equations  are  exactly  the   same.    Simply 

stated,   the  system  characteristics  are  identical.    System  behavior   Is 

identical,   and  since  the  same  behavioral  characteristics  are  present  in 

many   simple   Inertial   oscillators  the  behavior  is   repeated   in   other 

systems   and  may  be  conceptually  repeated  in  more  complex  systems.   One 


Figure  8:  Plot  of  the  electrical  oscillator  system  storages  vs. 
time.  Both  the  storage  of  potential  and  kinetic  energies 
oscillate  at  a  frequency  determined  by  the  pathway  flow 
coefficient  k4 .  The  determination  of  k4  is  always 
(2.0»TT»f  )  where: 
I 

Lo=3.1415 

fl=frequency  of  intrinsic  oscillation. 
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Figure  9:  The  mechanical  mass,  spring,  damper  oscillator.  Diagram 
a)  the  mechanical  diagram  indicating  system  components  and 
positional  relationship,  b)  the  energy  circuit  diagram 
indicating  system  components,  energy  storages,  and 
component  positional  relationship. 
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x=change  in  mass  oosition 
m=mass 

f  =intrinsic  frequency  of 
the  osci 1 lator 


k  =spring  constant 


dx/dt=l/D^*(S-X*k^)-(kg/D2)*X-(k^/m 
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important  concept  to  learn  from  these  inertial  oscillators  is  that 
behavior  or  identical  function  can  result  from  the  combination  of  a 
variety  of  different  components. 

Power  transfer  and  frequency  selectivity  are  important  concepts 
fundamental  to  understanding  the  operational  characteristics  of 
inertial  oscillators.  Power  transfer  (the  power  available  to  the 
load  from  the  source)  for  both  the  electrical  and  mechanical  inertial 
oscillators  is  governed  by  a  set  of  identical  behavioral 
characteristics.  Figures  10a  through  lOf  show  the  important  power 
transfer  characteristics  of  the  electrical  inertial  oscillator.  Figures 
lOd  to  lOf  indicate  that  for  the  same  source  potential  the  system  power 
is  distributed  differently  with  decreased  loading.  The  load  power,  the 
power  that  is  transfered  from  the  source  to  the  load,  is  maximized  when 
the  source  resistance  equals  the  load  resistance  (figure  lOe).  When 
balance  occurs,  source  power  (the  power  lost  within  the  source  as  a 
necessity  for  the  transfer  of  energy  from  the  source  to  the  load)  equals 
load  power,  producing  a  50*  transfer  of  power  form  source  to  load. 
Balance  in  this  context  describes  the  ratio  of  the  load  resistance  to 
the  source  resistance,  a  perfectly  balanced  system  occurring  when  the 
source  and  load  resistances  are  equal.  Figures  10a  through  10c  indicate 
that  the  potential  energy  stored  within  the  system  increases  as  the 
loading  is  decreased.  With  less  energy  required  by  the  load  the  system 
is  capable  of  storing  more  energy. 

Figures  11a  through  llf  also  indicate  that  the  power  available  to 
the  load  is  both  frequency  and  balance  dependent.  Figures  lid  through 
llf  show  that  the  power  transfered  from  the  source  to  the  load  is 
maximized   when   the  energy  delivered  is  of  the  same   frequency   as   the 


F  i  gure   1 0  : 


Resistive  power  transfer  characteristics  tor  the  bandpass 
electical  oscillator.  The  greatest  power  transfer  occurs 
when  the  resitance  of  the  source  equals  the  resistance  of 
the  load.  Diagram  (a)  oscillator  potential  and  kinetic 
storages  for  larger  source  resistance,  (b)  oscillator 
potential  and  kinetic  storages  for  equal  source  and  load 
resistance,  (c)  oscillator  potential  and  kinetic  storages 
for  larger  load  resistance,  (d)  source  and  load  power  for 
larger  source  resistance,  (e)  source  and  load  power  for 
equal  source  and  load  resistances,  (f)  source  and  load 
power  for  larger  load  resistance. 
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F I gure   1  1  : 


Frequency  power  characteristics  for  the  bandpass 
oscillator.  The  greatest  power  transfer  occurs  when  the 
frequency  of  the  source  equals  the  frequency  of  the 
intrinsic  freqency  of  the  oscillator.  Diagram  (a) 
oscillator  potential  and  kinetic  storages  for  source 
frequency  greater  than  the  oscillator  frequency,  (b) 
oscillator  storages  for  source  frequency  equal  to  the 
oscillator  frequency,  (c)  oscillator  storages  for  source 
frequency  lower  than  the  oscillator  frequency,  (d)  source 
and  load  power  for  source  frequency  greater  than 
oscillator  frequency,  (e)  source  and  load  power  for  source 
frequency  equal  to  oscillator  frequency,  (f)  source  and 
load  power  for  source  frequency  lower  than  oscillator 
frequency . 
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Inertial  oscillator's  internal  frequency.  When  the  source  frequency  is 
faster  or  slower  than  the  intrinsic  frequency  of  the  inertial  oscillator 
the  oscillator  operates  at  the  source  frequency.  Comparison  of 
Figures  11a  through  lie  indicate  that  as  the  source  frequency  decreases, 
i.e.  increasing  period  the  power  transferred  to  the  load  reaches  a 
maximum  when  the  source  frequency  equals  the  oscillator  intrinsic 
frequency,  f^igures  12a  and  12b  summarize  the  results  of  the  previous 
plots  of  load  power  with  changing  load  resistance  (Ri)  and  with  changing 
source  frequency  (fs),  respectively.  The  greatest  power  transfer  for 
changing  load  resistance  (RI)  occurs  when  the  source  resistance  (Rs)  and 
load  resistance  (RI)  are  equal.  The  greatest  power  transfer  for  changes 
in  source  frequency  (fs)  occurs  when  the  source  frequency  equals  the 
inertial  oscillator's  intrinsic  frequency  (fo).  When  the  two 
frequencies  are  equal  the  system  is  said  to  be  in  resonance,  and  that 
particular  frequency  is  called  the  resonant  frequency  (fr).  Moreover,  at 
resonance  the  circuit  operates  most  efficiently  with  the  least  source 
power  lost  for  the  most  load  power  delivered  (Figure  12b).  Ultimately, 
the  maximum  power  transfer  from  source  to  load  occurs  when  the  system  is 
at  resonance  and  is  in  balance.  The  maximum  storage  of  energy  within 
the  system  occurs  when  the  load  resistance  is  large  (the  least  loading) 
and  when  the  period  of  the  source  energy  is  long. 

The  same  power  transfer  characteristics  apply  to  the  mechanical 
inertial  oscillator,  however,  the  components  are  different.  Within  the 
mechanical  oscillator  maximum  power  transfer  occurs  when  the  damping  of 
the  source  equals  the  damping  of  the  load,  and  when  the  system 
oscillation  equals  the  source  oscillation.  Regardless  of  the  system 
(electrical  or  mechanical)  the  behavior  for  both  systems  is  the  same. 


Figure  12:  Summary  of  power  characteristics  for  the  bandpass 
oscillator.  Plot  (a)  source  and  load  power  with  changing 
load  resistance,  (b)  source  and  load  power  with  changing 
source  frequency.  The  greatest  power  transfer  occurs  when 
source  and  load  resistance  and  frequency  are  equal. 
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The  preceding  paragraphs  and  figures  have  presented  tne  concepts 
for  power  transfer  and  energy  selectivity  with  a  limited  reliance  upon 
mathematics.  However,  the  electrical  inert'al  oscillator  lends  'tself 
readily  to  a  mathematical  presentation  of  these  concepts.  figure  13 
again  shows  the  electrical  and  energy  circui+  diagrams  for  the  tank 
oscillator  (often  called  a  bandpass  oscillator  because  of  its  property 
of  passing  only  a  range  of  frequencies  to  the  load).  Equations  1 
through  11  show  the  mathematical  relationships  for  the  system,  in  both 
languages.  The  range  of  frequencies  that  are  passed  to  the  load  is 
often  limited,  for  practical  use,  to  a  specific  subset  of  frequencies 
called  the  circuit  bandwidth.  The  bandwidth  depends  solely  upon  the 
load  resistance  (Rl),  the  inductor's  (L)  value,  and  a  constant  equal  to 
2.0*  TT  .  Useful  power  transfer  occurs  for  frequencies  that  del  i  ver 
potentials  greater  than  or  equal  to  70.7?  of  the  source  potential 
(Ekeland,  1981).  The  frequencies  that  deliver  these  potentials  are 
within  the  bandwidth  of  the  inertia!  oscillator  circuit.  The  high  and 
low  cutoff  frequencies  for  the  bandwidth  are  calculated  from  tne  values 
of  the  system  components.  Equations  1  and  3  describe  the  calculations 
necessary  to  determine  the  bandwidth  of  the  circuit  shown  in  Figure  18. 
The  accompanying  equations  4  and  11  indicate  the  relationship  of  circuit 
components  necessary  to  calculate  resonant  frequency.  Resonant 
frequency  for  an  inertial  oscillator  defines  the  center  of  its  bandwidth 
because  it  is  the  frequency  that  transfers  the  most  energy  and  power. 

The  concept  of  maximum  power  transfer  as  a  result  of  energy 
selectivity  is  established  for  the  electrical  and  mechanical 
oscillators.  The  previous  paragraphs  and  figures  describe  the  general 
characteristics  of  this  phenomenon.   The  following  sections  describe  the 


F  igure   1 3 : 


The  diagram  shows  the  mathematical  relationships 
for   the  bandwidth  of  the  electrical  oscillator.   Diagram 

a)  the  electrical  circuit    and  corresponding  equations, 

b)  the  energy  circuit  and  corresponding  equations.  The 
equations  describe  the  same  bandwidth  because  the  two 
diagrams  are  mathematically  identical.  The  advantage  of 
the  energy  diagram  is  the  clearity  of  presentation  for 
components  and  storages  not  as  clearly  indicated  by  the 
electrical  diagram. 
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extension  of  these  concepts  into  the  development  of  a  autocata I ytic 
oscillator,  and  through  the  use  of  the  biological  oscillator  present  the 
concept  of  energy  selectivity  for  biological  systems. 

Autocatal ytic  Inertia!  Oscfllator  Model 

Electrical  and  mechanical  inertial  oscillators  exhibit  some  of  the 
important  characteristics  required  for  biological  rhythmic  behavior, 
however  they  do  not  self-adjust  nor  do  they  exhibit  growth  with  time. 
They  do  possess  an  Internal  timing  between  their  potential  (P)  and 
kinetic  (K)  components,  and  exhibit  time  relationships  between  their 
internal  components  and  an  external  source.  While  these  qualities  are 
necessary  for  biological  oscillators  they  fall  far  short  of  the 
complexity  generated  by  the  addition  of  nonlinear  autocata I yt  Ic  behavior 
within  the  Inertial  oscillator  itself. 

The  autocatal ytic  inertial  oscillator  incorporates  a  non-linear 
Internal  feedback  within  the  oscillator,  and  as  a  result,  produces  a 
phase  adjusting  oscillator  capable  of  sensing  environmental  sources  and 
remaining  tuned  to  those  sources.  Timing  functions  within  the 
autocata I ytic  inertial  oscillator  also  grow  with  time  and  might 
therefore  be  utilized  to  Initiate  action  with  Increased  age. 

Figure  14  is  an  energy  circuit  diagram  of  the  autocatal ytic 
inertial  oscillator  developed  for  this  study.  Equations  1  and  2 
present  the  mathematical  relationship  of  the  inertial  oscillator  system 
components.  The  /  /  indicates  the  use  of  the  absolute  value,  and  is 
necessary  because  the  energy  flows  within  the  Inertial  oscillator  are 
capable  of  completely  reversing  their  directional  flow.  Environmental 
sources  normally  have  some  constant  energy  about  which  the  available 
energy  varies.   However,  the  inertial  oscillator  senses  rhythmic  change 


Figure  14:  The  autocata I ytic  clock.  Intrinsic  rhythm  is  made  self- 
adjusting  by  the  addition  of  internal  feedback. 
Furthermore,  the  growth  with  time  of  the  potential  and 
kinetic  energy  storages  might  provide  an  explanation  for 
c I ock- i n i ti ated  age  functions  i.e.  reproductive  behavior. 
Timing  initiators  also  include  nocturnal  activity, 
seasonal  timed  activity,  and  environmentally  sensed 
initiation  of  particular  processes  i.e.  flowering  and 
fruiting. 


56 


dP/dt=k^*J^*/P/-k^*J^*/P/-k2*P-k2*;P 


j^=j/(i+kQVP/: 

S=A*sin(.^t) 


k,=Frequency  setting  coefficient  =  j  ■ 
.  =2*-*f. 

0       1 

f=osci 1 lator' s  intrinsic  frequency 


57 

and   as   a  result  the  /  /   is   required  to  ensure  that   flow   directions 
within  the  inertial  oscillator  remain  in  proper  phase  relationship. 

Autocatal yt ic  inertial  oscillator  timing  functions  include  five 
useful  timing  initiators.  First,  the  summation  of  potential  and  kinetic 
energies  within  the  inertial  oscillator  can  be  utilized  as  a  product  in 
an  energy  transformation  or  work  gate  function.  Figure  15  is  an  energy 
diagram  depicting  an  autocatal yt ic  inertial  oscillator  as  part  of  a 
energy  transformation  or  work  function.  The  system  is  energy-selective 
as  a  result  of  producer's  inertial  oscillator  activity,  and  this 
particular  concept  is  presented  by  simulation  plots  in  the  next  section. 
The  model  is  provided  to  familiarize  the  reader  with  the  location  of 
components  in  an  energy-selective  production  module.  Source  energy  to 
the  system  is  flow  limited,  and  the  system  has  a  recycle  of  (3)  through 
(N).  Flow  F2  indicates  that  the  recycle  of  nutrients  requ'res  some 
energy  loss  from  the  pathway.  Figure  16  shows  the  concept  for  the 
second  timing  initiator.  The  timing  of  the  potential  and  kinetic  tanks 
(P)  and/or  (K)  might  be  utilized  as  time  indicators  or  process  function 
initiators.  For  example,  the  hunting  initiator  for  nocturnal  predators 
or  the  tidal  feeding  initiator  for  clams,  oysters,  and/or  mussels.  For 
the  third  timer,  the  timing  combination  of  tanks  (P)  and  (K)  combine 
with  an  ability  to  sense  a  source  rhythm  to  begin  or  initiate  process 
functions.  Examples  of  this  characteristic  could  be  the  flowering  of 
plants  or  the  thickening  of  animal  furs  during  winter.  Figure  17  is  a 
model  depicting  the  concept  of  timing  with  a  source  and  internal 
inertial  oscillator  combination  as  the  timer.  Fourth,  the  amplitude  of 
tanks  (P)  and/or  (K)  are  employed  to  initiate  some  age  process  or 
function.   One  example,   for  this   type  of  Inertial  oscillator  activity 


F  igure   1 5  : 


Autocata 1 ytic  clock  as  part  of  a  workgate.  The  system 
uses  energy  selectivity  to  acquire  input  energy  and  to 
generate  biomass.  The  diagram  and  fol lowing  plots  are 
presented  only  to  familiarize  the  reader  with  the  use  of 
the  autocatal ytl c  clock. 
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Figure    16:   Timing    initiator    for    potential   component  of    the 

autocatai yt ic  clock.   This  particular   timing  initiator 

might   be  employed   by   a  biological   unit  to  initiate 
nocturnal  activity. 
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Figure  17:  Biological  timimg  using  the  source  rhythm  as  well  as  both 
the  potential  and  kinetic  components  of  the  autocata I yt ic 
c I ock . 
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might  be  the  initiation  of  a  particular  reproductive  behavior  which 
results  from  the  autocata  I  yt  ic  growth  of  tani<  (P)  and/or  (K).  When  one 
or  both  of  the  tanks  reach  a  predetermined  level,  the  autocata I yt ic 
inertial  oscillator  could  initiate  the  production  of  hormones  that  In 
turn  produce  timed  reproductive  behavior.  Figure  18  is  a  model  for 
the  age  process  initiator  of  the  autocata I yt ic  inertial  oscillator. 
The  fifth  timer  results  from  the  combination  of  some  or  all  of  the 
previous  functions  and  could  be  utilized  to  simulate  complex  forms  of 
behavior.  The  model  shown  in  Figure  19  uses  all  the  previous  timing 
functions  to  produce  a  complex  timing  behavior  capable  of  timing  both 
diurnal  as  well  as  annual  rhythms  within  the  environment  and  produces 
autocata I yt Ic  inertial  oscillator  behavior  that  can  take  advantage  of 
that  knowledge. 

Autocata I yt I c  inertial  oscillator  phasing  with  an  external  source 
plays  an  important  role  In  the  behavior  of  an  individual  system  within  a 
given  environment.  Therefore,  the  phase  relationship  of  the  potential 
(P)  and  kinetic  (K)  of  the  Inertial  oscillator  in  relation  to  the  source 
is  Important.  If  an  individual  biological  system  is  to  remain  in  tune 
with  Its  environment  it  must  adapt  to  changes  in  source  phasing. 

Phase-shift  experiments  are  employed  to  test  the  self-phasing  of 
biological  systems.  These  experiments  substitute  laboratory  phase- 
shifted  signals  for  normal  environmental  rhythms  and  examine  the  effects 
of  the  new  false  signal  on  the  system.  As  an  example,  Webb  (1950) 
studied  a  group  of  fiddler  crabs  (Uca)  to  analyze  phase  shifting  In  the 
crab's  color-rhythms.  Webb  substituted  darkness  during  normal  light, 
and  light  during  normal  darkness  producing  a  12-hour  phase  change  In 
the  crab's  color-rhythm.   The  crabs  subjected  to  the  phase-shifted- 


gure  Id:  Biological  timing  as  a  result  ot  the  ^nplitude  of 
autocatal yt ic  clock  rhythms.  This  tim'ng  mignt  be  used 
for  reproductive  behavior.  The  timing  is  a  response  to 
the  age  of  the  biological  unit. 
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Figure  19:  Biological  timing  as  a  comb i tat  ion  of  the  four  pro-  'ous 
timing  initiators.  This  type  of  timing  produces  complex 
behavior  by  timing  the  environment  as  wel I  as  the 
Individual  timing  characteristics  of  the  unit's 
autocatalytc  clock. 
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signal  altered  their  intrinsic  rhythms  to  self-align  with  the  new 
(albeit  false)  environmental  signal.  Therefore,  any  model  of  a 
biological  rhythm  generator  must  be  able  to  self-phase  with  its 
env  i  ronmenta I  . 

Figure  20  is  a  plot  showing  the  autocata I yt ic  'nertial  oscillator's 
ability  to  self-phase  when  presented  with  a  phase-shifted  external 
signal.  The  simulation  begins  with  a  source  supply  that  "s  sine- 
phased.  The  autocatal yt ic  inertial  oscillator  is  in  tune  with  the 
source,  and  storages  (P)  and  CK)  are  increasing  and  In  normal  phase 
relation.  Between  time  2.0  and  4.5  the  source  rhythm  Is  eliminated. 
The  Inertial  oscillator  oscillation  continues  with  the  same  phasing, 
however,  the  levels  in  tanks  (P)  and  (K)  decrease.  This  decrease  is  the 
result  of  depreciation  losses  via  pathway  F2,  and  Is  required  by  the 
second  law  of  thermodynamics.  Beginning  at  time  4.5  the  source  Is 
reinstated  with  an  cosine-phased  oscillation  I.e.  a  90  degree  or  it /4 
radian  shift  from  the  original  source  rhythm.  Oscillator  potential  (P) 
and  kinetic  (K)  components  experience  a  significant  perturbation  wlth'n 
the  next  two  cycles,  after  which  time  the  autocata I yt ic  inertial 
oscillator  self-tunes  to  the  source's  new  phasing.  Within  the  next 
seven  cycles  the  inertial  oscillator  completely  re-establishes  an  in- 
phase  relationship  with  the  source. 

Frequency   setting   and   damping  of  Internal  rhythms  also   play   an 

important   role  In  an  system's  ability  to  maintain  rhythinic  behavior   or 

to   efficiently  utilize  Its  internal  rhythmic  activity  In  a   competitive 

situation.   The  frequency  setting  coefficient  k3  (in  the  model  shown  in 

2 
Figure   14)   Is  always  equal  to  uj  .    In  radian  units    is  2.0*Tr  »fo, 

where   ( fo)   is  the   intrinsic  frequency  of  the   Inertial   oscillator. 


gure  20:  Plot  indicating  the  ability  of  the  autocata I yt ic  clock  to 
self-adjust  to  changes  in  phase  of  an  environmental 
source.  The  damping  during  source  elimination  results 
from  a  second  law  loss  of  energy  via  pathway  F3,  and  is  a 
requirement  for  all  systems — including  internal  clocks. 
Full  recovery  as  a  result  of  the  phase  shift  'ntroduced 
into  the  source  required  seven  oscillations  or  periods. 
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Coefffcient   k3   is   extremel/   important   to   the  operation   of    the 

autocata 1 yt ic  inertial  oscillator  because  it  controls  the  period  of   the 

autocata I yt ic   inertial   oscillator's  internal  oscillation.    The   exact 

method  by  which  biological  oscillators  set  their  rhythms  is   still   not 

clearly   understood,   however   the  ability  to  identify  these  periods   of 

biological  rhythms  does  exist.    The  periods  for  intrinsic  rhythms   are 

known  for  many  plants  and  animals.    Within  humans  internal  free  runn'ng 

biological   rhythms   exhibit  a  period  range  from  19  hours  to  25.8   hours 

and   were  identified  from  cave  and  bunker   studies   (Pa Imer , 1 976) .   The 

following   passage   is   an  example  of  how   biological   system   intrinsic 

rhythms  may  be  identified. 

Over  a  decade  ago,  a  young  speleologist  became 
obsessed  with  the  idea  that  It  was  scientifically 
Important  for  him  to  live  alone,  underground,  sans 
clocks,  for  two  months.  He  chose  an  inhospitable 
cave  In  the  '^rench  Alps  where  he  lived  at  the  375- 
foot  level.  Here  the  temperature  hovered  at  32 
degrees  (F),  the  relative  humidity  remained 
unchanged  at  100!?,  and  the  darkness  was  complete 
save  for  a  small  battery-powered  light  In  his  tent. 
Each  time  he  awoke,  ate,  or  prepared  to  retire  he 
called  over  a  field  telephone  to  a  surface  camp, 
where  his  lonely  words  and  the  times  of  his  calls 
were  recorded.  He  claimed  that  the  inexorable  cold 
and  dampness  reduced  his  body  temperature  to  less 
than   97   degrees   (F),    and   he   was   constantly 

threatened  by  avalanches  and  cave-Ins  still,   he 

held  out  for  the  sake  of  science  and  whiled  away  his 
time  writing  a  best-selling  novel  about  his 
subterranean  adventures,  anxiety,  and  building 
I i  b I  do. 

Throughout  his  underground  stay,  he  tried 
mentally  to  keep  track  of  the  passage  of  time  on  the 
surface.  When  the  men  in  the  surface  camp  informed 
him  on  September  14  that  his  experiment  was  over,  he 
thought  it  was  only  August  20.  His  subjective 
judgement  of  the  passage  of  time  had  been 
exceedingly  sluggish;  mentally,  he  had  lost  25  days! 
A  major  problem  in  his  estimation  of  the  passage  of 
time  was  the  fact  that  he  commonly  Imagined  that  he 
was  taking  a  short  siesta  after  his  midday  meal, 
where  in  reality  he  had  been  awake  for  16+  hours  and 
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his  "nap"  lasted  approximately  8  hours.  However, 
his  living  clock,  as  evaluated  by  the  times  of  his 
retiring  and  awaking  phone  calls  to  the  surface,  had 
ignored  his  mental  confusion  and  guided  his  body 
functions  all  the  while,  measuring  off  periods  of 
activity  and  sleep  that  total  just  longer  than  a 
day:  24  hours  and  30  minutes  on  average  (Palmer, 1975 
p.  125)  . 

The  above  passage  illustrates  the  intrinsic  rhythm  setting  and 
damping  properties  of  internal  biological  rhythm  generators  (living 
clocks  as  described  by  Palmer).  The  frequency  of  biological  oscillators 
is  determined  under  constant  conditions,  and  many  experiments  are 
directed  solely  at  determining  these  rhythms.  Interestingly,  one  way  to 
find  the  intrinsic  rhythm  of  an  electrical  or  mechanical  inertial 
oscillator  is  to  remove  the  oscillator  from  its  energy  supply  and  time 
the  period  of  its  intrinsic  oscillation. 

The  duration  of  oscillator  activity,  the  amount  of  time  an 
autocatal yt ic  inertial  oscillator  rhythm  remains  active  when  isolated 
from  its  energy  cue  by  placing  it  into  constant  conditions  is  also 
important.  The  ability  of  a  biological  system  to  adapt  to  interuptions 
in  energy  supply  is  critical.  With  heavy  damping  an  autocata I ytic 
inertial  oscillator  loses  its  rhythm  (in  constant  conditions)  within  a 
short  period  of  time.  However,  with  light  damping  the  autocatal ytic 
inertial  oscillator  could  maintain  a  rhythmic  pattern  for  long  periods 
of  time — as  in  the  above  passage  for  two  months  or  longer.  Damping 
determines  the  longevity  of  the  autocatal ytic  inertial  oscillator's 
activity  when  environmental  conditions  become  constant.  One  example  of 
an  interruption  in  an  environmental  source  supply  is  a  drought  in  a 
norma  I  I y  wet  area . 

With  these  characteristics  in  mind,  the  simulation  of  biological 
rhythmic  activity  by  an  autocata I yt ic  inertial  oscillator   requires  that 
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the  internal  autocata I yt i c  mertial  oscillator  characteristics  be 
adjustable.  Coefficients  k3  and  k2  satisfy  these  requirements.  Pathway 
coefficient  k3,  sets  the  intrinsic  rhythm,  while  coefficient  k2  adjusts 
the  dampir^.  Damping  of  the  autocata I yt ic  inertial  oscillator's  free- 
running  rhythm  is  the  direct  result  of  the  depreciation  losses  within 
the   inertial  oscillator,   and  in  the  model  is  accounted  for  by   pathway 

F2. 

Biological  Models  with  Autocatal ytic  Inertia!  Oscillators 

The  previous  section  presented  results  for  the  energy  circuit 
autocatalytic  inertial  oscillator.  Development  of  this  oscillator  is 
an  important  initial  step  in  the  analysis  of  energy  selectivity  and  its 
effect  on  the  power  maximization  and  power  transfer  within  biological 
systems.  The  following  model,  plots,  and  figures  illustrate  the  effect 
of  energy  selectivity  on  power  transfer  and  biomass  generation  for  a 
producer  with  an  internal  autocatalytic  inertial  oscillator.  The  model 
Is  Intended  to  Illustrate  the  effects  of  energy  selectivity  on  a 
hypothetical  Diological  system. 

Figure  21  Is  an  energy  diagram  for  an  energy-selective  producer. 
Table  3  presents  the  storage  Initial  conditions,  the  pathway 
coefficients,  a  descriptive  comment,  and  the  equations  for  the  model  in 
Figure  21.  The  production  of  biomass  (pathway  Fl )  Is  regulated  by  the 
activity  of  the  producer's  internal  autocatalytic  inertial  oscillator. 
Since  the  autocatalytic  Inertial  oscillator  is  energy  selective,  the 
work  gate  or  energy  transformation  function  is  also  energy-selective. 
Pathway  F2  Is  the  autocatalytic  Inertial  oscillator's  sensor  and 
provides  information  to  the  inertial  oscillator  about  the  rhythmic 
changes   that  occur   in   Its  environment.   Pathways  F6   and  F7   are 


Figure  21:  Energy  diagram  for  an  energy-selective  producer.  The 
clock  produces  an  internal  rhythm  required  for  the 
production  of  unit  biomass.  Equations  1  through  4 
describe  the  component  relationships  in  mathematical 
I anguage. 
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depreciation  pathways.  A  linear  depreciation  is  provided  for  by  pathway 
F6  and  represents  the  normal  depreciation  losses  within  the  inert'al 
oscillator.  These  losses  include  physical  losses  such  as  cellular 
breakdown.  Pathway  F7  represents  the  deprec'at'on  losses  experienced 
by  the  autocata I yt ic  inertial  oscillator  however  these  loss  are  opposed 
by  a  kinetic  inertial  backpressure.  K'netic  inertial  backpressure  is 
felt  as  a  byproduct  of  the  flow  in  pathway  F7. 

The  total  energy  of  the  autocata I yti c  inertial  oscillator  is  the 
sum  of  the  inertial  oscillator's  potential  and  kinetic  energies,  and  's 
employed  as  the  work  gate  regulator.  When  the  source  rhythm  is 
outside  of  the  producer's  range  of  energy  acceptance  the  producer 
rejects  the  available  energy  and  production  of  biomass  decreases. 

An  energy-selective  producer  increases  its  biomass  or  net 
production  within  a  defined  range  of  energies,  i.e.  a  b'ological  energy 
acceptance  niche.  The  term  bandwidth  is  used  to  describe  the  energy 
acceptance  range  for  the  electrical  and  mechan'cal  inertial  oscillators. 
The  analogy  for  bandwidth  in  ecological  systems  is  MICHE.  Figure  22a 
shows  an  energy  acceptance  niche  for  a  producer  with  a  free-runn"ng 
internal  biological  inertial  rhythm  having  a  24  hour  period  of 
oscillation.  The  width  of  an  energy-acceptance  niche  could  be  employed 
to  determine  whether  the  producer  is  a  specialist  or  simple  general  ist. 
The  models  in  Figure  22b  show  the  relative  frequency  relationships  for 
the  simulations  employed  to  produce  the  energy-acceptance  niche  In  the 
Figure  22a.  The  source  rhythms  are  depicted  in  the  circular  source 
symbols  while  the  producers'  intrinsic  rhythms  are  depicted  in  the 
bullet  symbols  used  for  energy  circuit  producers.  The  figure  clearly 
indicates  that  the  most  biomass  is  generated  when  the  source  and 
producer  rhythms  are  equal. 


Figure   22: 


Plot  indicating  the  niche  width  of  an  energy  selective 
producer  with  an  intrinsic  rhythm  of  one  day.  The  niche 
might  be  called  an  energy  acceptance  niche  since  it 
defines  the  range  of  frequencies  or  rhythms  within  which 
the  producer  accepts  energy  input.  This  figure  was 
produced  by  plotting  the  steady  state  biomass  for  five 
simulations  of  the  model  in  figure  21.  For  each 
simulation  the  internal  rhythm  of  the  producer  was  held 
constant  at  1  period  per  day,  while  the  source  rhythm  was 
varied  during  the  simulations.  The  source  periods  ranged 
from  0.01  to  100  oscillations  per  day.  The  mini-diagrams 
in  (b)  are  only  to  aid  in  distinguishing  between  the 
relative  period  changes.  Since  the  change  in  source 
period  varied  by  order  of  several  magnitudes  plotting 
these  rhythms  for  presentation  was  impractical. 
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Energy  Selective  Blomass  Generation 

This  section  shows  the  results  of  the  analysis  for  the  energy- 
selective  generation  of  biomass.  The  producer  simulations  show  that 
energy  transfer  for  the  generation  of  producer  biomass  occurs  at  the 
frequencies  within  the  energy-acceptance  niche  of  the  producer. 
Frequencies  outside  of  the  producer's  energy-acceptance  n'che  generate 
no  biomass. 

The  plots  "n  Figure  23a  show  the  effect  of  energy-selectivity  on  a 
producer's  biomass  for  source  energies  less  than,  equal  to,  and  greater 
than  the  producer's  intrinsic  rhythm.  The  plots  clearly  show  the 
producer's  capability  for  accepting  or  rejecting  rhythmic  energies. 
Biomass  attains  its  greatest  level  when  the  producer's  free-running 
autocata I yt ic  inertial  oscillation  equals  the  frequency  of  the  available 
energy  (plot  line  a).  When  the  available  source  energy  oscillates  with  a 
rhythm  outside  of  the  producer's  acceptance  range  the  producer  rejects 
the  energy  (plot  lines  b  and  c)  .  Figure  23b  shows  the  producer's 
internal  autocatal yt ic  inertial  oscillator  generates  a  rhythmic  behavior 
within  the  producer  that  is  either  in  or  out  of  step  with  the  source 
rhythm,  thereby  selectively  accepting  or  rejecting  the  opportunity  to 
utilize  available  energy.  The  three  mini-models  show  the  relative 
period  relationships  for  the  source  and  intrinsic  producer  rhythms. 
Whether  the  former  or  the  latter  is  true,  the  producer  by  its 
autocatal ytic  oscillator-regulated  activity  defines  its  acceptance  niche 
within  which  it  accepts  energy.  Interestingly,  biomass  generation  about 
the  producer's  fundamental  rhythm  appears  relatively  symetrical. 


Figure  23:  Simulation  plots  and  relative  frequency  diagram  for  the 
energy-selective  producer  in  Figure  21.  Plot  line  Ca) 
shows  the  growth  in  producer  biomass  for  source  energy 
inside  the  producer's  energy-acceptance  niche.  Plot  lines 
(b)  and  (c)  show  the  decrease  in  biomass  for  source 
energies  at  frequencies  outside  of  the  producer's  energy- 
acceptance  niche.  The  greatest  biomass  generation  occurs 
for  supply  energy  at  the  free-running  rhythm  of  the 
producer's  intrinsic  rhythm. 
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Power  Transfer  for  an  Energy  Selectfve  Producer 

Figures  24a  and  24b  show  the  producer's  unit  power  and  power 
transfer  for  the  three  plots  in  Figure  23.  The  maximum  power 
transferred  by  the  producer  to  the  remainder  of  the  system  (via  F3) 
occurs  when  the  rhythms  of  the  producer  and  the  system  are  matched.  For 
all  source  rhythms  outside  the  producer's  energy-acceptance  niche  power 
and  energy  are  rejected  and  power  transfer  is  minimized.  Producers  that 
develop  acceptance  niches  too  far  removed  from  their  environmental 
signals  do  not  possess  the  ability  to  pass  on  an  appreciable  amount  of 
power  from  an  external  source  to  the  remainder  of  the  system.  A 
producer  capable  of  attracting  more  energy  and  power  and  transferring 
more  power  to  the  next  component  in  a  system  is  also  contributing  more 
toward  the  power  maximization  of  the  system  as  a  whole.  This  concept 
plays  an  important  role  in  the  "Maximization  of  Power"  as  introduced  by 
Lotka,  1922  and  later  refined  by  Odum,  1976. 

Essentially,  the  principle  of  maximum  power,  describes  a  system 
structure  that  gets  the  most  power  for  available  energy  input  while  at 
the  same  time  utilizing  that  input  power  most  efficiently.  Efficient 
power  transfer,  therefore,  becomes  an  important  characteristic  for 
stable  system  behavior.  Were  systems  to  reward  selfish  producers  or 
consumers  the  system  would  be  relegating  itself  to  receiving  less  energy 
from  components  lower  on  the  food  chain  (than  would  be  expected  by  the 
Principle  of  Maximum  Power),  and  as  a  consequence  would  support  a  more 
I  imi ted  food  chai  n  . 

One  way  for  a  system  to  maximizes  its  power  input  while  insuring 
stability  over  a  wide  range  of  inputs  is  to  develop  producers  and 
consumers   that  create  or  occupy  specific  spaces  or  niches.   A  system 


Figure   24: 


Plot  of  power  transfer  for  the  energy-selective 
producer.  Power  transfer  to  the  load  (the  rest  of  the 
system)  via  pathway  F3  reaches  a  maximum  when  the  internal 
rhythm  of  the  producer  equals  the  external  oscillation  of 
the  source. 
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structure  composed  of  units  that  provide  the  most  power  transfer  to  the 
next  unit  in  the  food  chain,  with  the  most  efficient  use  of  power,  not 
only  generates  more  diversity  but  also  provides  for  the  maximum  number 
of  units  with  the  most  efficient  use  of  input  energy.  Simply  stated, 
the  best  way  to  maximize  system  power  is  to  provide  a  means  of  selecting 
the  greatest  number  of  most  efficient  biological  units.  Diversity, 
created  by  components  that  are  energy-selective,  produces  both 
stability  and  maximum  power  transfer.  A  second  inspection  of  the  plots 
in  Figures  24a  and  24b  indicates  that  energy  input  and  power  transfer 
to  the  next  component  within  the  system  (via  pathway  F3)  are  most 
efficiently  accomplished  when  the  rhythm  of  the  source  energy  is  the 
same  as  the  producer's  fundamental  or  free-running  rhythm. 

Figure  25a  through  25c  shows  the  biomass  and  power  levels  for  an 
energy-selective  producer  with  an  internal  period  much  larger  than  the 
rhythmic  period  of  the  available  source  energy  (Figure  25d).  While  the 
outcome  of  this  plot  indicates  death  for  the  producer  the  initial  peak 
in  biomass  level  might  indicate  an  early  success'onal  producer.  The 
unit  exhibits  an  ability  to  attract  energy  rapidly  but  for  only  a 
relatively  short  period  of  time.  The  result  of  such  behavior  favors  the 
early  development  of  the  unit  but  provides  the  unit  with  a  short 
existence,  therefore  a  large  amount  of  biomass  is  generated  early  in 
system's  development.  A  behavioral  pattern  that  allows  units  to 
generate  large  quantities  of  biomass  quickly  but  provides  only  a  short 
existence  could  provide  the  system  with  a  means  of  readying  an 
inadequate  environment  for  future  more  complex  biological  activity. 


Figure  25:  Plot  of  biomass,  unit  power,  and  power  transfer  vs.  time 
for  an  energy-selective  producer  with  an  intrinsic  rhythm 
much  larger  than  the  available  energy  oscillation. 
Producing  units  which  exhibts  this  type  of  behavior  (the 
rapid  generation  of  biomass  with  no  longevity)  are  often 
called  early  successionai  forms. 
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Competf tlon  Between  Energy  Selective  Units 

Within  this  section  are  the  figures,  plots,  and  tables  indicating 
the  effects  of  energy-selective  competition  on  a  variety  of  hypothetical 
biological  units.  The  intent  is  to  present  a  theory  for  competition 
based  upon  the  energy-selective  behavior  exhibited  by  one  or  all  of  the 
competing  units.  The  following  nine  simulations  present  various 
combinations  of  energy-selective  producers  in  competition  with  each 
other  for  limited  resource. 

Figure  26  is  an  energy  diagram  depicting  competition  between  two 
energy-selective  producers.  Table  4  presents  the  storage  initial 
conditions,  a  descriptive  comment,  pathway  coefficients,  and  equations 
for  the  model.  This  diagram  is  used  for  the  majority  of  the  nine 
simulations,  and  is  only  repeated  within  the  simulation  sections  were 
the  model  is  modified  for  that  particular  simulation.  Each  simulation 
represents  a  change  in  either  the  energy  selectivity  of  one  or  both  of 
the  competing  units  or  a  change  in  the  source  supply  rhythmic  pattern. 
Each  unit  employs  an  autocata I yt i c  internal  oscillator  that  governs  the 
unit's  behavior.  The  inertia!  oscillators  sense  the  environment  and 
produce  an  internal  activity  which  defines  the  units'  energy  acceptance 
n  iche. 

The  following  description  of  the  model  is  general  and  each 
simulation  incorporates  some  new  combination  for  competition  between  the 
energy-selective  units  within  the  model.  However,  the  description 
provides  a  general  undertsand i ng  of  the  model's  design. 

Producer  #1  (the  producer  with  the  shortest  free-running  internal 
period)  a  specialist  in  comparison  to  unit  #2.  However,  the  two 
producers   are  identical  in  all  respects  except  for  coefficients  k7   and 


Figure  26:  Energy-selective  competition.  The  model  has  two  energy- 
selective  units  in  competition  with  one  another.  The 
models  is  employed  to  simulate  a  series  of  competitive 
situations.  In  all  situations  the  model  shows  the 
producer  with  the  shortest  period  for  its  internal  rhythm 
in  the  upper  left  corner.  Thus  the  specialists  is  in  the 
upper  left  corner  and  the  general ist  is  'n  the  lower  right 
corner . 
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k16  which  set  the  free-running  internal  rhythms  for  their  respective 
autocatal ytic  inertial  oscillators.  The  source  rhythmic  pattern  is  a 
composite  cycle  having  a  diurnal  and  annual  component  to  its  available 
energy.  Each  producer,  as  a  result  of  its  internal  behavior,  senses  the 
source's  rhythmic  pattern  and  selectively  regulates  the  production  of 
unit  blomass  based  upon  its  acceptance  or  rejection  of  available  energy. 
The  two  producers  compete  for  available  energy,  however  their  internal 
behavioral  characteristics  determine  which  producer  is  more  capable  of 
utilizing  part  or  all  of  the  energy. 

Power  transfer  for  both  units  is  determined  by  the  amount  of  energy 
flow  per  unit  time  along  four  pathways  FO,  F8,  F3,  and  FIO.  The  energy 
flow  per  unit  time  (power)  via  pathway  FO  is  the  power  required  by  unit 
#1 .   Pathway  F8   provides   the  power  demanded  by   producer  42.  The 

summation  of  energy  within  pathways  F3  and  FIO  determines  the  amount  of 
power  transfers  from  each  unit  to  its  respective  load,  i.e.  the 
remainder  of  the  system.  The  diagram  shows  pathways  F3  and  FIO  as  open 
ended  and  indicates  that  these  pathways  transfer  energy  and  power  to 
other  components  or  to  a  more  complex  system.  The  producer  that 
acquires  the  most  input  power,  utilizes  that  power  most  efficiently, 
and  transfers  the  greatest  amount  of  power  to  its  load  (the  rest  of 
the  system)  survives  in  competition,  and  successfully  fills  a  niche 
within  the  system. 

Simulation  1;   Competitive  Exclusion  of  a  Simple  Generallst  by  a 
Spec lal Ists 

This   simulation  shows  the  effects  of  energy-selective  competition 

for   input   rhythmic   energy  between  a  specialist   (the   short   rhythmed 

producer)  and  a  generallst  (the  long  rhythmed  producer).   The  model   in 
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Figure   27a   shows  the  relative  period  relationships  between  the   source 

and  the  competing  producers.     The  niche  plots  in  Figure  27b  indicates 

the   relative   (for   visual   aid   only)   or   graphical   range  of    the 

competitors'  energy-acceptance  niches.    The  source  energy  is   available 

only   at   the  intrinsic  rhythm  of  the  specialist  producer   and   provides 

that  unit  with  a  competitive  advantage. 

Figure   28   shows  the  simulation  plots  for  competition  between   two 

energy-selective  producers.    Producer  ^1,   the  specialist,   employs   an 

intrinsic  free-running  rhythm  period  of  1  day.   The  general ist,  unit  ^2, 

utilizes  an   intrinsic  rhythm  with  a  period  of  50  days.    The  specialist 

clearly   outcompetes  the  general  ist  reaching  a  maximum  biomass  level   of 

6 
8.0  X  10   Cal  . 

The  competition  between  these  two  producers  results  in  the 
competitive  exclusion  of  the  general  Ist  by  the  specialist.  The  initial 
success  of  the  general  ist  is  attributable  to  the  longer  period  of  its 
internal  free-running  behavior.  Units  with  longer  internal  free- 
running  rhythms  initial ly  begin  to  store  a  greater  amount  of  energy, 
thereby,  producing  an  initial  rapid  increase  in  their  biomass. 
However,  the  consequence  of  this  behavior  over  a  long  period  of  time 
results  in  the  exclusion  of  unit  with  the  longer  internal  period  because 
the  production  of  biomass  is  not  supported  by  the  interaction  of  source 
energy  with  internal  behavior.  The  two  rhythmic  energies  which  combine 
to  generate  biomass  for  the  long  period  producer  (its  autocata I ytic 
inertial  oscillator  and  the  source)  counteract  each  other  more  often 
than  they  complement  one  another,  leading  to  the  exclusion  of  this 
producer . 


Figure  27:  Simulation  1  system  relative  period  and  energy-acceptance 
niche  diagrams  (a)  system  relative  periods  (b)  system 
relative  energy-acceptance  niches.  The  diagram  shows 
only  the  relative  system  periods  and  niches  not  the  actual 
periods  and  energy-acceptance  niche  curves.  The  d'agram 
is  presented  to  aid  in  understanding  the  intrinsic  rhythms 
of  the  competing  producers.  The  short  period  producer  has 
a  rhythm  of  1  day  while  the  long  period  producer  has  a 
rhythm  of  50  days.  The  source  energy  is  cyclic  at  a 
period  equal  to  that  of  the  short  period  producer  and  as  a 
result  the  short  period  producer  links  well  with  the 
source. 
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Figure  28:  Plots  for  competition  between  a  specialist  and  a  simple 
general ist.  Plot  (a)  the  amount  of  biomass  is  greater  for 
unit  #1  because  it  specializes  in  the  frequency  of  the 
energy  source,  (b)  the  power  of  unit  ,^1  is  greater  than 
that  of  unit  til  because  unit  #2's  rhythmic  behavior 
produces  activity  too  far  removed  from  that  of  the  source, 
(c)  power  transfer  for  unit  'ri\  is  clearly  greater  than 
that  of  unit  ti2  and  therefore  unit  ti\  is  of  more 
importance  to  the  next  larger  system. 
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Figure  28b  and  28c  indicate  the  effects  of  ccmpetitron  on  the  DOwer 

demand  and   transfer  for  both  units.   The  specialist  receives  more  power 

from   the  source  and  as  a  result  transfers  more  power  to  the  system   via 

pathway  F3 .    Figure  23b  indicates  that  power  demand  for  specialist  un^t 

5 
it\       reaches   a  maximum  at  9.8  x  10   Cal/year,   while  the  level   for   the 

5 
general  ist   decreases  to  1 .0  x  10   Cal/year.    The  power  demand   by   the 

specialist   levels   and   thereafter  remains  relatively  constant   at   its 

maximum  rate  of  consumption.    Figure  28c  shows  that  the  power   transf-'T 

by   unit  #1  far  exceeds  that  transfered  by  unit  H2    and  reaches  a  maximum 

4 
of   7.0  X  10   Cal/year,   or  approximately  1.0%   of  the  power  utilized   by 

the   unit.    The  power  transfer  for  the  aeneralist  reaches  a  maximum  of 

4 
2.5   X   10    Cal/year   within  the  first  months   of   the   s'mula+ion   and 

4 
thereafter  decreases  to  a  min'mum  level  of  1 .0  x  10   Cal/year. 

If   the  specialist  producer  were  a  member  of  a  more  complex   system 

its  greater  power  transfer  would  make  it  more  useful  to  the  next   level 

within  the  system,   however  for  this  particular  mini-model  the  remainder 

of  the  system  is  not  modelled.   The  effects  of  power  transfer  within  web 

systems  are  presented  later,   for  the  moment,   it  is  (nore  beneficial  for 

the   system   as  a  whole  if  units  lower  on  the  food  chain   transfer   more 

energy/power  to  the  next  level  in  the  food  chain. 

Simulation  2:  Competition  Between  Identical  Producers 

This  simulation  was  generated  to  analyze  tne  effects  of  energy- 
selective  competition  between  two  identical  energy-selective  producers. 
The  two  competitors  compete  for  a  source  of  energy  that  has  a  rhythmic 
pattern  that  is  identical  to  the  internal  rhythms  of  both  the  producers. 

Figure's  29  and  31  show  the  model  and  simulation  plots  for  energy 
selective  competition  between  producers  which  are  alike.    The  producers 


Figure  29:  Energy  diagram  of  competition  between  identical  producers. 
The  diagram  shows  the  model  for  competition  between  like 
producers.  The  two  producers  are  shown  side-by-side 
because  they  are  the  same. 
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are   arranged   directly  beside  one  another  to  indicate  their   identical 

behavior.   Both   producers  employ  the  same  free-running  "nternal   period 

for    their   autocata I yt ic   inertial   oscillation,    and    the    pathway 

coefficients  for  all  pathways  and  inertial  oscillators  are  identical  for 

both   producers.   The  relative  period  relationships  ("igure  30a)   are 

shown   for   the   source  and  both  competitors  and  clearly   indicates   the 

equality   between  the  competing  units.    Figure  30b  indicates   that   the 

producers'  energy-acceptance  niches  are  identical. 

Clearly,   the  competition  is  harmful  to  both  units,  and  the  biomass 

of  both  units  is  greatly  reduced  by  the  competition.   Biomass  levels  for 

5 
both   units  reach  a  maximum  of  4.0  x  10   Ca I .    These  levels  represent  a 

decrease  of  approximately  50'^  from  levels  obtainable  by  each   individual 

producer   without   the  presence  of  competition.    However,   both   units 

survive  and  coexist  within  an  Identical  niche.    It  should  be  noted  that 

the   units   are   exactly   identical  in   all   respects.    Were   the   flow 

coefficients   for  the  various  pathways  within  the  producers   to   differ, 

while   their  autocata I yt I c  inertial  oscillators  remained  identical,   one 

or  the  other  would  benefit  more  from  the  competition.   The   ultimate 

effect  of  the  competition  would  then  depend  on  the  competitive  advantage 

provide  by  other  aspects  of  the  unit's  Internal  composition. 

Figures  31b  and  31c  show  the  power  required  by  the  individual  units 

and  the  power  transferred  through  these  units  to  the   system.   Power 

5 
demand  by  both  units  reaches  a  maximum  of  5.0  x   10   Cal/year.   Both 

4 
units  transfer  4.1   x  10  Ca 1  /year  via  pathways  F3   and   FIO.   Power 

transfer   to  the  system  remains  high  or  unaffected  because  both   units 

transfer  an  identical  amount  of  power  to  the  remainder  of  the  system. 


Figure  30:  Simulation  2  system  relative  period  and  energy-acceptance 
niche  diagrams  (a)  system  relative  periods  (b)  system 
relative  energy-acceptance  n'ches.  The  producers  have  a 
rhythm  of  1  day  as  does  the  source  energy.  3oth  producers 
I  ink  we  I  I  with  the  source  and  compete  equal  ly  wel  I  tor  the 
ava  i I ab le  energy. 
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Figure  31:  Plots  for  competition  between  identical  energy-selective 
producers.  Plot  (a)  biomass.  Plot  (b)  unit  power.  Plot 
(c)  power  transfer.  The  producers  act  identically  and 
therefore  one  can  not  outcompete  the  other.  The  overal I 
power  transfer  to  the  system  remains  high  because  each 
unit  transfers  1/2  the  power  transfer  of  an  individual 
unit. 
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The  producers  exhibit  identical  biomass  production  and  power 
transfer  because  their  intrinsic  behavior  is  identical.  Although  no  two 
producers  are  identical  in  every  respect  the  model  and  simulation  serve 
to  show  the  effects  of  intraspeci f ic  competition  in  its  purest  form. 

Simulation  5;  Competitive  Exclusion  of  a  Specialist  by  a  Simple 
General  1 st 

This  simulation  was  generated  to  analyze  the  effects  of  energy- 
selective  competition  between  two  energy-selective  producers.  The  model 
was  employed  to  investigate  the  competitiveness  of  a  general  ist  that  is 
in  tune  with  its  environment.  Unlike  the  previous  simulation  the 
specialist  does  not  link  well  with  the  environment  but  the  general  ist 
does  link  effectively  to  the  environment.  The  producers  have  energy- 
acceptance  niches  generated  by  their  internal  rhythmic  behavior.  The 
source  rhythm  oscil I ates  at  the  same  frequency  as  the  long  rhythmed 
producer  or  general ist  and  the  relative  period  relationships  between  the 
source  and  the  two  producers  are  shown  in  Figure  32a.  Relative  period 
diagrams  are  used  here  to  indicate  that  there  is  a  marked  difference 
between  the  periods  of  the  system  components — the  diagrams  do  not  shown 
the  actual  relationship  between  the  rhythm  periods  because  the 
difference  is  often  in  orders  of  magnitude  an  could  not  be  plotted  in  a 
resonable  manner  for  presentation.  A  graphical  representation  of  the 
relative  energy-acceptance  niches  (also  distorted  and  intended  for 
visual  aid  only)  for  the  two  competing  producers  is  shown  in  Figure  32b. 
The  diagram  indicates  that  the  source  input  energy  is  available  at  the 
same  frequency  as  the  internal  rhythm  of  the  general ist  producer  unit. 

Figures  33a  through  33c  show  the  effects  of  energy-selective 
behavior   for  two  producers.    The  general ist  employs   the   longest 


Figure  32:  Simulation  3  system  relative  period  and  energy-acceptance 
niche  diagrams  (a)  system  relative  periods  (b)  system 
relative  energy-acceptance  niches.  The  diagram  shows 
only  the  relative  system  periods  and  niches  not  the  actual 
periods  and  energy-acceptance  niche  curves.  The  short 
period  producer  has  a  rhythm  of  1  day  while  the  long 
period  producer  has  a  rhythm  of  20  days.  The  source 
energy  is  cycl  ic  at  a  period  equal  to  that  of  the  long 
period  producer  and  as  a  result  the  long  period  producer 
links  well  with  the  source. 
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intrinsic  rhythm  of  the  two  competing  producers,   and  easily  outcompetes 

the  specialist  because  the  free-running  rhythm  of  the  general  ist  matches 

the  period  of  the  source.   The  specialist  unit  develops  such   a  rapid 

internal  behavior  it  produces  an  extremely  narrow  niche  outside  of  which 

it  rejects  available  energy.   On  the  other  hand,  the  genera! ist  develops 

an   energy   acceptance  niche  close  to  that  of  the  available  energy   and 

therefore  accepts  energy  from  the  source.   The   specialist's   behavior 

results  in  its  exclusion  by  the  general  ist.   The  biomass  plots  in  Figure 

33a   indicate  that  the  general ist  easily  outcompetes   the   specialist, 

6 
producing  a  maximum  biomass  level  of  8.0  x  10   Cal/year. 

Rejection  of  energy  by  the  specialist  does  not  imply  that  the  unit 
willingly  accepts  its  own  demise.  A  unit  that  does  not  accept  an 
available  energy  does  so  because  of  physical  and/or  psychological 
reasons.  Two  examples  of  physical  rejection  include:  seed  rejection  by 
birds  with  inadequate  bill  size  and  the  lack  of  photosynthet ic 
production  for  solar  radiation  with  wavelengths  outside  of  the  vis'ble 
light  range.  An  example  of  psychological  rejection  is  the  abstinence 
from  the  ingestion  of  food  at  the  noon  meal  when  dieting. 

Power   requirements   and   transfer  characteristics   for   the  two 

producers  are  indicated  in  Figures  33b  and  33c.   The  general ist  producer 

5 
acquires   the   most   power   from  the  source  (9.8   x   10    Cal/year)   and 

4 
transfers  more  power   (7.8  x  10   Cal/year)  to   the  remainder  of   the 

system.   The  plot  of  unit  power  (Figure  30b)  also  indicates  that  the 

competition  was  relatively  easy  for  the  general ist.   Within  a  few  months 

it  captures  a  majority  of  the  available  energy  while  at  the   same  time 

competitively  excluding  the  specialist  producer. 


Figure  33:  Competitive  exclusion  of  a  specialist  by  a  simple 
general ist.  Plot  (a)  biomass  levels  for  the  competing 
units  indicates  that  unit  #1  easily  outcompetes  unit  #2, 
(b)  unit  power  for  unit  /^1  reaches  its  maximum  rapidly 
because  it  has  a  greater  competitive  edge,  (c)  power 
transfer  for  unit  #1  is  significantly  greater  than  unit  #2 
and  as  a  result  unit  //I  is  of  greater  importance  to  the 
next  larger  system. 
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Simulation  4;  Competf tf on  of  two  Producers  for  a  Limfted  Energy  Within 
Their  Niche-Acceptance  Ranges  but  not  at  Either  Producers'  Free-Running 
Rhythm 

This  section  describes  the  results  for  the  computer  simulgtion  of 
energy-selective  competition  between  two  competing  producers  for  input 
energy  at  a  frequency  which  does  not  oscillate  af  either  producers 
intrinsic  rhythm.  The  diagram  in  Figure  34a  shows  the  relative  period 
relationship  for  the  source  and  competitors'  rhythmic  oscillations.  The 
energy-acceptance  niches  for  the  two  competing  producers  are  shown  in 
Figure  34b,  and  are  graphically  represented  to  depict  the  special  ist- 
generalist  concept  for  short  and  long  period  energy-selective  units. 

Figures  35a  through  35c  show  the  results  from  the  fourth  simulation 
of  the  competition  model  shown  in  Figure  26.  The  simulation  indicates 
the  effects  of  competition  at  the  extremities  of  two  producers'  niches. 
Producer  #1  (the  specialist)  employs  an  intrinsic  free-running  rhythm 
with  a  period  of  1  day.  Producer  #2  (the  general ist)  uses  a  intrinsic 
rhythm  with  a  free-running  period  equal  to  20  days,  and  the  source 
provides  a  yearly  rhythm  composed  of  periods  equal  to  2.5  days  or  52 
hours.  To  simulate  the  effects  of  a  extended  day-night  cycle  each  2.5 
day  period  generates  no  source  energy  supply  for  1/2  of  the  period  (1 .25 
days).  The  source  energy  rhythmic  period  (2.5  day  cycle)  is  therefore 
closer  to  the  internal  rhythm  behavior  of  the  specialist. 

Figure  35  indicates  that  the  specialist  producer  outcompetes  the 
general Ist,  although  the  competition  's  intense.  The  specialist  is 
competing  for  energy  at  the  edge  of  Its  niche  width  and  can  not  rapidly 
suppress  the  general ist  unit.  Although  the  general ist  uses  a  longer 
internal  period  for  behavior  and  competes  well  during  the  early  stages 
of  competition   it  also  competes  at  the  limits  of   its   niche.   The 


F  i  gure   34 : 


Simulation  4  system  relative  period  and  energy-acceptance 
niche  diagrams  (a)  system  relative  periods  (b)  system 
relative  energy-acceptance  niches.  The  short  period 
producer  has  a  rhythm  of  1  day  while  the  long  period 
producer  has  a  rhythm  of  20  days.  The  source  energy  is 
cyclic  at  a  period  of  2.5  days  and  is  not  cyclic  at  the 
intrinsic  period  of  either  producer.  The  source  period  is 
closer  to  the  rhythm  of  the  specialist  producer  and  as  a 
result  the  specialist  is  able  to  acquire  more  of  the 
ava  i I ab le  energy . 
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competitive  edge  exhibited  by  the  specialist  results  because  the   source 

provides   energy   within  its  narrow  niche  width.   The  general  ist  on   the 

other  hand,    exhibits  the  same  type  of  early  behavior  exhibited  by   the 

individual   producer   in  Figure  25,  i.e.  rapid  biomass  generation   as   a 

result   of  early  competitiveness.    The  behavior  of  the  general  ist   unit 

makes   it   compete   well  for  short  periods  of  time  but  not  in   the   long 

term.    This  type  of  behavior,   rapid  early  biomass  generation  and  short 

term  competitiveness  (perhaps  because  of  energy  selectivity),    Is  often 

described  as  an  early  successional  stage. 

Figure   35b   shows   the  unit  power  demand  for   the   two   producers. 

During   the  two  year   period  of   simulation   the  general  ist   producer 

exhibits  a  higher  unit  power  requirement  reaching  a  maximum  level  at  the 

5 
beginning  of  the  simulation  at  9.8  x  10   Cal/year  but  unit  power  rapidly 

reduces  as  the  competition  progresses.   Early  in  the  simulation  the  unit 

power  for  the  general  ist  exhibits  a  noise  or  rapid  fluctuation  of   power 

Input.   This   fluctuation  results  from  the  acceptance  of  smaller  period 

energies   generated  by  the  competitive  Interaction  of  the  two  units   and 

serves   as   an  indication  of  the  ability  of  energy   selective   units   to 

utilize   energies   other  than  their  primary  energies  within  their   niche 

range.   The  unit  power  requirements  of  specialist  Initially   decreases 

because  of  the  rapid  albeit  short  lived  growth  In  the  biomass   for   the 

general Ist.   The  behavior  of  specialist  provides  the  unit  with  a   long 

term  strategy   for  competitive  success,   at  least  for  this  competitive 

scenario.   The  general ist  transfers  more  power  to  the  system  during  the 

4 
short   term   (reaching  a  maximum  of  5.5  x   10   Cal/year,   however,   the 

competitive  advantage  of  specialist  begins  to  provide  the  system  with   a 

more  stable  power  transfer  for  the  long  term. 


Figure  35:  Competition  at  energy  niche  extremities.  Plot  (3) 
biomass  levels  for  the  two  competitors,  the  specialist 
clearly  outcompetes  the  general  ist  'n  the  long  term.  Plot 
(b)  unit  power  the  specialist  eventually  attracts  more 
power  from  the  source.  Plot  (c)  power  transfer  'n  wh'ch 
the  general  ist  transfers  more  power  initially,  however  the 
specialist  transfers  more  power  in  the  longterm.  These 
plots  reflect  the  fact  that  the  source  energy  is  at  a 
rhythm  more  closely  coupled  to  the  specialist,  but  not 
coupled  at  the  specialist's  intrinsic  rhythm. 
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Sfmulatlon  5;  Competition  of  Sfmilar  Producers  With  Overlapping  Niches 

In  this  section  the  effects  of  energy-selective  competition  between 

two   producers  with  close  intrinsic  rhythmic  patterns  is  analyzed.    The 

intent   of  the  simulations  was  to  determine  if   competitive   coexistance 

might   occur  as  a  result  of  overlapping  energy-acceptance   niches.    The 

diagram   in   Figure  36a  shows  the  relative  period  relationships  for   the 

oscillations  of  the  source  and  the  two  competi tiors .    Figure  36b   shows 

the  relative  energy-acceptance  niches  for  the  two  competi tiors ,  and  also 

indicates   that  the  source  input  energy  is  at  the  same  frequency  as  that 

of   the  specialist  or  short  rhythmed  producer.    The   diagram   indicates 

that  the  energy-acceptance  niches  for  the  two  competitors  overlap   one 

another   and  as  a  result  the  competition  for  input  energy  is  between  two 

producers  that  have  a  differential  ability  to  acquire  the  same  energy. 

Figures  37a  through  37c  present  the  effects  of  competition  between 

two   energy-selective   producers  with  intrinsic  behavior   patterns   that 

have   periods  that  are  only  slightly  different  from   one   another.    The 

producers   could  be  either  specialists  or    general ists  since  the  concept 

in  this  case  Is  a  relative  term.   The  free-running  period  of  producer  ^1 

is   1   day  and  the  Internal  period  for  producer  #2  is  1.05   days.   Both 

units  survive  the  competition  producing  a  competitive  coexistence  where 

both   units  experience  decreased  biomass  production  as  a  result  of   the 

6 
competition.   Biomass  levels  reach  6.3  x  10  Ca  I   for  unit  ,^1   and  1.9  x 

6 
10   Ca  I  for  unit  ^'2.    Producer  ff]    acquires  the  most  unit   power  (7.9  x 

5  4 

10   Cal/year)   and   transfers  more  power  (7.0  x   10   Cal/year),   via 

pathway  F3.   Producer  H2   obtains  or  requires  less  unit  power  from  the 

5 
source   (2.5  x  10  Cal/year)  than   does  unit  #1  and  also  transfers   less 

4 
power   (2.1  x  10  Cal/year),  via   pathway  FIO. 


Figure  36:  Simulation  5  system  relative  period  and  energy-acceptance 
niche  diagrams  (a)  system  relative  periods  (b)  system 
relative  energy-acceptance  niches.  The  short  period 
producer  has  a  rhythm  of  1  day  while  the  long  period 
producer  has  a  rhythm  of  1.05  days.  The  source  energy  's 
cyclic  at  a  period  of  1  day  and  is  cyclic  at  the  intrinsic 
period  of  short  producer.  The  source  period  is  also  close 
to  the  rhythm  of  the  long  period  producer  and  as  a  result 
the  long  period  producer  is  able  to  compete  better  for  a 
portion  of  the  available  energy. 
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Figure  37:  Competition  of  similar  producers  with  overlapping  energy 
niches.  Plot  (a)  biomass.  Plot  (b)  unit  power.  Plot  (c) 
power  transfer.  The  plots  indicate  that  competitive 
coexistance  occurs  when  competition  is  within  both 
competitors'  energy  acceptance  ranges. 
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Some  interesting  aspects  of  this  particular  simulation  include:  the 

quick   allocation   of   input  power  and  the  increased   diversity   of   the 

system.   Figure  37  indicates  that  this  type  of  competition  provides  both 

units  with  a  stable  relationship,   and  does  so  rapidly.    Power  transfer 

for  both  units  stabilizes  as  their  biomass  levels  stabilize.   The  system 

receives   the   same   amount  of   power   transfer   as   it   did   from   the 

simulations   of  this  model  that  produced  a  competitive  exclusion  of   one 

unit;   but   the   system   experiences  a  higher   diversity.    Producer   #2 

experiences   a  small   growth   and   therefore   is   considered   to  be   a 

subordinate  species.   However,  if  the  intrinsic  rhythm  of  producer  H2    is 

further   removed   from  that  of  producer  #1  and   the   source,   but   still 

within  the  competitive  range  of  unit  #2's  niche  biomass  production  would 

decrease.   The  resulting  competition   then  places  producer  It!      to   the 

classification  of  a  rare  species. 

Simulation  6:  Competition  for  Flow  Limited  Energy  Compossed  by  the 
Addition  of  Two  Fundamental  Cycles  About  a  Supporting  Mean  Energy  Level 

The  purpose  of  this  simulation  was  to  analyze  the  effects  of 
competition  for  energy  that  is  available  from  the  same  source  but  that 
has  frequencies  that  may  be  employed  by  more  than  one  competing  unit. 
The  intent  of  the  simulation  was  to  discover  if  the  addition  of  a  second 
component  rhythm  to  the  source  might  allow  the  source  to  support  a 
second  energy-selective  producer  In  a  competitive  environment.  Also 
could  the  addition  of  a  second  rhythm  to  the  source  al low  for 
competitive  coexistance  of  producers  shown  in  earlier  simulations  to 
competitively  exclude  when  only  one  source  rhythm  was  available? 

Figure  38a  shows  the  relative  period  relationships  for  the  source 
and  the  two  competing  producers.   The  diagram  does  not  show  the   actual 
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period  relationships  and  is  intended  only  as  an  aid  to  discern  the 
relative  differences  in  frequency  for  the  system  oscillations.  The 
energy-acceptance  niches  for  the  two  competitions  are  shown  In  figure 
38b.  The  diagram  also  indicates  that  the  source  rhythm  is  a  compination 
of  the  producers'  internal  rhythms  and  as  a  result  the  competitors  each 
attempt  to  acquire  energy  from  the  same  source  only  at  different  energy 
f requenc  ies . 

Figures  39a  through  39c  indicate  the  effects  on  biomass,  unit 
power,  and  power  transfer  for  two  competing  producers  that  are  in 
competition  for  source  energy  that  is  flow  limited  and  composed  by  the 
addition  of  two  energy  cycles.  The  source  energy  is  compossed  of  two 
rhythms.  The  first  rhythmic  period  is  equal  to  1  day,  while  the  second 
source  period  is  equal  to  20  days.  The  two  competing  producers  each 
have  only  one  internal  free-running  rhythmic  period.  Producer  //I,  the 
specialist,  uses  an  intrinsic  rhythmic  pattern  with  a  period  equal  to  1 
day.  Producer  H2,  the  general ist,  uses  an  intrinsic  period  equal  to  20 
days.  Consequently,  both  producers  have  the  ability  to  link  or  interact 
with  one  of  the  source's  fundamental  energy  components. 

Figure  39  clearly  indicates  that  the  specialist  successfully 
outcompetes  the  general ist — even  though  the  general ist  possess  the 
longer  internal  rhythm.  The  figure  also  indicates  a  fluctuation  in  the 
biomass  of  the  specialist  at  a  frequency  of  20  days.  The  fluctuation 
indicates  that  the  specialist  feels  the  longer  period  energy,  although 
it  may  act  only  to  produce  the  fluctuation.  Previous  simulation  plots 
indicate  that  the  energy-selective  specialist  could  not  accept  longer 
period  energy.  However,  the  shorter  period  intrinsic  behavior  of  the 
specialist  producer   provides  that  unit  with  a  competitive  advantage. 


Figure  38:  Simulation  6  system  relative  period  and  energy-acceptance 
niche  diagrams  (a)  system  relative  periods  (b)  system 
relative  energy-acceptance  niches.  The  short  period 
producer  has  a  rhythm  of  1  day  while  the  long  period 
producer  has  a  rhythm  of  20  days.  The  source  energy  is 
cyclic  at  a  periods  of  both  1  and  20  days  and  is  therefore 
cyclic  at  the  intrinsic  periods  of  both  producers.  The 
short  period  rhythm  supplies  more  energy  to  the  short 
period  producer  than  does  the  long  period  rhythm  for  the 
long  period  producer.  The  greater  energy  input  for  the 
short  period  producer  makes  it  a  better  competitor. 
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Figure  39:  Competition  for  additive  energy  frequencies  by 
competitors  that  each  match  to  one  of  the  souce 
frequencies.  Plot  (a)  biomass.  Plot  (b)  unit  power.  Plot 
(c)  power  transfer.  The  plots  indicate  that  the 
competitor  with  the  shorter  period  'nternal  rhythm 
outcompetes  the  unit  with  the  longer  period  internal 
rhythm.  Therefore,  when  competition  occurs  for  two 
rhythms  from  the  same  energy  source  (with  the  same  energy 
quality)  the  competitor  that  links  to  the  shorter  period 
energy  receives  and  transfers  the  most  power  and  survives 
the  competition. 
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That   is   not  to  suggest  that  a  short  period   source   is   advantageous, 

however   for   the   same  energy  quality,  source's   with   shorter   periods 

provide  more   power.    The   short   period   intrinsic   behavior   of   the 

specialist   unit   allows  it  to  draw  more  energy  and  power  (of   the   same 

quality)   from   the   source  than  does  the   intrinsic   behavior  of   the 

general  ist.   The  competition  results  in  the  competitive  exclusion  of  the 

genera  list. 

Figures   39b  and  39c  show  the  unit  power  requirements  and   transfer 

capabilities  for  the  two  producers.   The  specialist  clearly  obtains  more 

5 
power   from  the  source  3.0  x  10   Cal/year  and  as  a  result  transfers  more 

4 
power   3.0   x   10   Cal/year  to  the   system.    The   general ist   producer 

5 
rejects  input  power  in  the  long  term  obtaining  1.0  x  10   Cal/year  and  as 

4 
a  result  transfers   less   power  1  .0  x   10   Cal/year   to  the   system. 

Although   the   general ist   exhibits   the   early   competitiveness   shown 

previously  for  long  period  competitors,  and  suffers  the  same  fate  as  the 

Its   previous  counterparts — exclusion  by  a  shorter  period,   more  adapted 

compet I  tor  . 

Simulation  7;  Competition  for  a  Flow  Limited  Energy  Compossed  of 
two  Fundamental  Cycles  by  Multiplication  About  a  Mean  Energy  Level 

This   section   describes  the  results  for   competition   between   two 

energy-selective  producers  competing  for  a  dual  rhythmed  input  energy 

from   the   same   source.   The   simulation   is   much   like   the   previous 

simulation   but   the   source   rhythm   construction   Is   different.    The 

difference   in   the  two  source  rhythms  Is  the  reason  for  this   analysis. 

The  intent  of  this  simulation  was  to  discover  if  the  method  used  by  the 

source  to  combine  the  two  or  more  of  its  rhythms  had  any  effect  on   the 

competitive  stance  of  the  competing  producers. 
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The  diagram  in  Figure  40a  shows  the  relative  period  relationships 
for  the  system,  while  the  niche  diadrams  in  Figure  40b  describe  the 
energy-acceotance  niches  for  the  competitors.  Figures  41a  through  41c 
show  the  biomass,  unit  power,  and  power  transfer  ca  ability  of  two 
competing  producers  with  different  fundamental  intrinsic  rhythmic 
behaviors.  The  source  energy  supply  is  a  combination  of  two  rhythms  by 
multiplication  and  are  supported  by  or  around  a  constant  mean  energy 
level.  This  simulation  is  identical  to  simulation  fib  with  the  exception 
of  the  way  the  two  fundamental  source  rhythms  are  combined.  The  first 
source  energy  rhythmic  period  is  equal  to  1  day  with  the  second  rhythm 
having  a  period  equal  to  20  days.  The  two  competing  producers  are 
identical  to  the  producers  in  the  previous  simulation. 

Figure  41  indicates  that  generation  of  a  source  energy  supply  by 
multiplication  of  Its  two  fundamental  frequencies  increases  the 
competitiveness  of  the  simple  general ist.  However,  the  result  is  the 
same  and  the  generallst  is  eventually  excluded.  The  combination  of  two 
energy  rhythms  by  multiplication  produces  more  interaction  rhythms  fhan 
does  combination  by  addition  and  as  a  result  provides  the  general ist 
(which  uses  the  longer  intrinsic  free-running  behavior  pattern)  with  a 
better  competitive  stance.  The  more  by-product  rhythms  generated  the 
more  possible  energy  there  is  available  for  producers  with  wider  niches. 
Since  the  general ist  employs  the  wider  of  the  two  niches  it  obtains 
more  energy  than  in  the  previous  simulation.  The  reason  the  general ist 
can  not  outcompete  the  specialist  results  from  the  fact  that  the 
specialist  obtains  more  energy  from  the  source's  daily  rhythm  than  does 
the  general ist  from  the  combination  of  the  interaction  rhythms  plus  the 
20  day  fundamental  rhythm.   The  competitiveness  of  the  simple  general ist 


Figure  40:  Simulation  7  system  relative  period  and  energy-acceptance 
niche  diagrams  (a)  system  relative  periods  (b)  system 
relative  energy-acceptance  niches.  The  short  period 
producer  has  a  rhythm  of  1  day  while  the  long  period 
producer  has  a  rhythm  of  20  days.  The  source  energy  is 
cyclic  at  a  periods  of  both  1  times  20  days  and  is 
therefore  cyclic  at  the  intrinsic  periods  of  both 
producers.  The  short  period  rhythm  supplies  more  energy 
to  the  short  period  producer  than  does  the  long  period 
rhythm  for  the  long  period  producer.  The  greater  energy 
Input  for  the  short  period  producer  makes  it  a  better 
competitor.  However,  the  side  frequencies  generated 
because  of  the  nonlinear  multiplication  of  source 
fundamental  frequencies  allows  the  long  period  producer  to 
acquire  more  energy  and  compete  more  effectively  for 
energy  than  in  the  previous  simulation 
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Figure  41:  Competition  for  multiplicative  energy  frequencies  by 
competitors  that  each  match  to  one  of  the  frequencies  of 
the  source.  Plot  (a)  biomass.  Plot  (b)  unit  power.  Plot 
(c)  power  transfer.  As  with  figure  34  the  competitor 
that  matches  with  the  shorter  period  energy  from  the  same 
quality  source  survives  the  competition.  However,  less 
power  and  energy  are  accepted  and  transferred  when  source 
energy  results  from  the  product  of  frequencies. 
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does  however  decrease  the  steady  state  biomass  and  power  transfer  levels 
obtained  by  the  specialist  unit. 

Figures  41b  and  41c  show  the  unit  and  power  transfer 
characteristics  for  the  two  competing  producers.  Initially  the 
general ist  accepts  the  most  unit  power,  however,  as  the  competition 
progresses  the  amount  of  unit  power  it  accepts  decreases.  The 
specialist  increases  its  demand  for  power  during  the  simulation  because 
of  its  intrinsic  behavior.  Since  the  specialist  free-running  rhythmic 
pattern  matches  the  daily  rhythm  of  the  source's  shortest  fundamental 
cycle  it  eventually  obtains  more  power  from  the  source  than  does  the 
general  ist.  Figure  41c  indicates  that  the  specialist  producer  transfers 
more  power  than  does  the  simple  general  ist.  The  intensity  of  the 
competition  and  the  demonstrated  ability  of  the  general ist  for  early 
competitive  success  does  provide  for  a  considerable  power  transfer  by  it 
during  the  simulation.  Ultimately,  the  specialist  transfers  more  power 
to  the  system  and  outcompetes  the  general ist  because  it  possess  the 
ability  to  accept  more  energy  from  the  source. 

Simulation  8:  Competition  Between  a  Two  Rhythm  Complex  Generallst  and  a 
Single  Rhythm  Specialist  for  a  Flow-Limited  Source  Composed  of  Two 
Fundamental  Cycles  by  Addition 

This  simulation  was  completed  to  analyze  the  effects  on  energy- 
selective  competition  for  producers  which  develop  two  or  more  intrinsic 
rhythms.  The  Intent  of  this  simulation  was  to  discover  if  the  addition 
of  a  second  Inertial  oscillator  linked  to  a  second  environmental  cue 
provided  the  producer  with  a  greater  competitive  advantage. 

The  diagram  in  Figure  42a  shows  the  relative  period  relationships 
for  the  system  components.  The  energy-acceptance  niche  diagram  (Figure 
42b)   shows  the  relative  niches  developed  as  a  result  the  energy- 


Figure  42:  Simulation  8  system  relative  period  and  energy-acceptance 
niche  diagrams  (a)  system  relative  periods  (b)  system 
relative  energy-acceptance  niches.  The  short  period 
producer  has  a  rhythm  of  1  day  while  the  long  period 
producer  has  a  rhythm  of  land  365  days.  The  source  energy 
is  cyclic  at  a  periods  of  both  1  and  365  days  and  is 
therefore  cyclic  at  the  intrinsic  periods  of  both 
producers.  The  short  period  rhythm  supplies  energy  of  one 
rhythm  to  the  short  period  producer.  The  short  and  long 
period  rhythms  supply  energy  to  the  long  period  producer. 
The  greater  energy  input  for  the  long  period  producer 
makes  it  a  better  competitor. 
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selective  activity  of  the  competitors.    The  dashed  line  represents   the 
relative  niche  for  the  dual  rhythmed  producer,   and  is  generated  by  the 
combination  of  its  two  internal  rhythms. 

Figure  43a  shows  the  biomass  levels  for  two  competing  producers. 
Producer  H\  (the  specialist)  utilizes  only  one  internal  inertia! 
oscillator  and  therefore  only  one  intrinsic  behavior  pattern.  Producer 
#2  (the  complex  general ist)  uses  two  internal  inertial  oscillators  and 
as  a  result  two  internal  rhythmic  behavior  patterns.  The  complex 
general  ist  may  possess  the  ability  to  attract  or  to  link  to  both 
fundamental  energies  of  the  source. 

The  source  is  composed  of  two  fundamental  rhythms.  The  first 
having  a  period  of  1  day  ,  and  the  second  having  a  period  of  365  days. 
The  figure  clearly  demonstrates  the  general  ist's  competitive  advantage, 
with  both  fundamental  source  rhythms  appearing  in  the  biomass  level  of 
the  producer.  The  specialist,  on  the  other  hand,  utilizes  the  shorter 
period  source  energy  but  is  competitively  exclude  by  the  complex 
general ist  because  Its  single  internal  rhythmic  pattern  does  not  provide 
an  adequate  behavior  to  obtain  as  much  of  the  available  power  as  does 
the  general ist's  dual  rhythm  internal  behavior. 

The  figure  also  indicates  that  early  in  the  competition  the 
specialist  increases  its  biomass  slightly.  This  result,  unlike  the 
early  success  of  the  one-rhythm  general ists  In  previous  simulation 
sections,  may  be  attributed  to  the  fact  that  annual  intrinsic  behavior 
(from  the  second  inertial  oscillator  of  the  general  ist)  does  not  produce 
a  significant  competitive  advantage  in  early  stages  of  the  competition. 
The  control  action  from  the  two  inertial  oscillators  of  the  general  Ist 
initially   resembles   the  controlling  behavior  of   the   short   period 


Figure  43:  Competition  between  a  complex  general  ist  and  a 
specialist.  Plot  (a)  biomass.  Plot  (b)  unit  power.  Plot 
(c)  power  transfer.  The  plots  'ndicate  the  competitive 
advantage  exhibited  by  a  complex  general ist.  The  complex 
general ist  has  two  internal  rhythms,  and  as  a  consequence 
obtains  more  power,  builds  greater  biomass,  and  transfers 
more  power  than  does  the  specialist.  The  complex 
general ist  times  its  environment  better  than  the  single 
rhythm  spec  i  a  I ist. 
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Inertia!  oscillator.  Consequently,  the  two  units  are  competitive 
equals  early  in  the  competition  after  which  time  the  general  ist's  annual 
autocata I ytic  inertial  oscillator  (longer  rhythmic  behavior)  provides  a 
competitive  advantage  that  eventually  excludes  the  specialist.  Figures 
43b  and  43c  show  the  unit  and  load  power  for  the  two  producers.  The 
unit  power  obtained  from  the  source  by  the  general ist  is  the  highest 
initial  power  of  any  of  the  simulations  using  a  dual  source  rhythmic 
pattern,  but  appears  to  reach  a  steady  state  about  a  mean  biomass  level 
approximately  equal  to  that  of  the  specialist  in  simulation  6.  Thus, 
the  combination  of  short  and  long  term  internal  behavior  provides  a 
competitor  with  the  best  competitive  advantage  insuring  survival.  Figure 
43c  indicates  that  power  transfer  from  the  general ist  to  the  system  far 
exceeds  that  of  the  specialist. 

The  advantage  offered  the  system  by  the  dual  inertial  oscillator 
producer  is  i ts  ab i 1 i ty  to  1  i nk  to  both  source  rhythms  or  to  either 
fundamental  rhythm  if  one  or  the  other  is  absent.  Thus,  the  complex 
general  ist  is  better  equipped  to  deal  with  source  energy  that  provides 
energy  or  energies  at  rhythms  with  periods  between  the  two  intrinsic 
patterns  of  the  producer.  The  timing  of  two  or  more  rhythms  al lows  the 
complex  general ist  to  accept  energies  from  a  wide  range  of  input.  When 
those  inputs  are  the  product  of  a  number  of  rhythms  the  complex 
general Ist  also  possesses  the  ability  to  accept  energies  formed  as 
sideband  energies  from  the  Interaction.  The  final  analysis  indicates 
that  the  complex  general Ist  Is  a  formidable  competitor. 

Models  for  a  Selected  Spartlna  Marsh 

The  following  section  contains  new  simulation  plots  for  a  Spartlna 
marsh  previously  simulated  and  presented  by  Odum,   Knox,   and  Campbell, 
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1983.    Additional   simulations  of  the  model  are  presented  by   Campbell, 
1984. 

The  original  mini-model  of  the  Kaituna  salt  marsh  (located  in  New 
Zealand)  consists  of  a  Spartlna  townsendll  storage,  a  storage  of  dead 
matter  as  sediment  carbon,  and  a  nutrient  storage  of  changeable  or 
labile  nitrogen.  The  two  primary  energy  sources  that  to  determine  the 
growth  of  Spart i  na  are  solar  energy  (external  forcing  function)  and 
available  nitrogen  (internal  storage).  Two  other  external  sources,  tide 
and  river  flow,  are  also  part  of  the  model. 

Figure  44  shows  the  mini-model  as  originally  evaluated  (Odum,  et 
al  .,  1983).  Equations  1  through  4  describe  the  mathematical 
relationships  between  the  various  system  components.  Table  5  provides 
the  initial  flow  and  storage  values  from  the  original  model  evaluation. 
The  Appendix  contains  the  CSI^P  computer  programs  used  to  s'mulate  the 
original  and  five  supplemental  Kaituna  mini-models. 

The  following  eight  simulations  of  the  Kaituna  mini-model  show  the 
effects  of  the  addition  of  a  autocata I yt ic  inertial  oscillator  to  the 
Spartlna  townsendii  production  unit.  Each  simulation  with  the  exception 
of  the  original  model  use  an  autocata I yti c  inertial  oscillator  as  a 
behavioral  control  within  the  Spartlna.  Simulations  5  and  6  incorporate 
a  second  producer  unit  that  is  in  competition  with  the 
Spartlna  townsendl I  .  Simulations  5  and  b  are  therefore  theoretical  and 
indicate  the  effects  of  energy  selective  competition  within  a  web 
system . 

Original  Kaituna  Salt  Marsh  Mini-model 

The  intent  of  this  simulation  was  to  duplicate  the  original 
simulations    for    the  Kaituna  marsh   model.    The   system   period 


Figure  44:  The  original  Kaituna  salt  marsh  mini-model.  The  SpartI na 
producer  receives  solar  energy  and  generates  biomass.  The 
biomass  is  recycled  to  the  nutrient  pool  which  in  turn  is 
used  to  produce  more  Spartina  biomass.  A  river  adds 
nutrients  to  the  system,  while  the  tidal  energy  removes 
nutrients  and  dead  organics  from  the  system.  Dead  organic 
material  is  recycled  by  the  action  of  consumers  within  the 
marsh . 
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Figure  45:  The  system  period  relationships  for  the  original  Kaituna 
marsh  model.  The  systam  has  a  annual  source  rhythm  and 
the  Spart i  na  producer  has  no  internal  rhythmic  pattern. 
The  system  is  not  energy-selective  and  the  producer  does 
not  sense  or  time  its  environmental  energy  source. 
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Figure  46:  Simulation  plots  tor  the  original  Kaituna  salt  marsh 
energy  model.  Plot  (a)  dead  organic  matter.  Plot  (b) 
Spartina  biomass  and  system  nutrients  Plot  (c)  Spartina 
unit  power.  Plot  (d)  Spartina  power  transfer.  The  system 
produces  increased  Spartina  biomass  that  exhibits  an 
annual  cycle  as  a  result  of  the  solar  driving  function. 
Dead  matter  increases  as  does  the  system  nutrient  pool. 
The  unit  input  power  for  the  Spartina  producer  increases 
rapidly  after  which  time  it  levels  and  exhibits  an  annual 
pattern  of  input  power.  The  power  transfer  for  the 
SpartI na  module  increases  more  gradually  than  the  input 
power,  however  it  to  levels  and  exhibits  a  slight  annual 
pattern . 
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relationships   are  shown  in  Figure  45.   The  source  has  an  annual  rhythm 

and   the  Spart tna  producer  has  no  intrinsic  rhythmic   pattern.   Figures 

46a   and  46b  show  the  levels  for  the  three  storages  of  the  Kaituna   salt 

marsh  system.   With  a  solar  input  exhibiting  of  a  annual  rhythm,  a  tidal 

input   with  approximately  a  fortnight  rhythm,   and  a  river  input  with   a 

annual   rhythm   the   system  experiences  an  increase   in  the   level   of 

standing  Spartlna,   dead  matter,   and  labile  nitrogen.    Although  labile 

nitrogen   seems   to  be  the  least  affected   storage — reaching   a  maximum 

-2 
level   of  1000  g  N  m   during  the  10  year  time  frame  of  the   simulation. 

The   figures  also  indicates  that  Spart I na  standing  stock  increases   from 

-2  -2 

523  g  C  m   to  1100  g  C  m   ,   while  dead  matter  increases  from  5226  g  C 

-2  -2 

m   to  9200  g  C  m   . 

Figures  46c   and  46d  show  the  unit  power   requirement   and   power 

transfer   capability  of  the  Spart i  na  producer.   The  unit  power  of   the 

Sparttna,   in  solar  Calories,  achelves  an  initial   maximum  power  of  3900 

-2     -2  -2     -2 

Cal   m    day   thereafter  leveling  at  3000  Cal  m    day   .   The  power 

transfer  of   the  Spart I na  producer  increases  gradually  to  a  maximum 

-2    -2 
sustained  level  of  17  Cal  m   day 

Kaituna  Salt  Marsh  Mini-model  with  the  Annual  Input  Solar  Energy  Rhythm 
Constructed  of  a  Series  of  Diurnal  Rhythms 

This  simulation  was  employed  to  introduce  a  diurnal  rhythm  into  the 

solar   source  of  the  Kaituna  model.   The  original   Kaituna  mini-model 

utilizes   a  solar  Input  energy  with  an  annual  rhythmic  pattern  for   its 

primary   external    forcing   function.    The   following   simulation 

Incorporates   into  the  original  solar  input  a  diurnal   rhythm,   thereby, 

producing  a   solar   input  with  a  day-night  cycle  combined  within   the 

annual  pattern.   The  new  rhythm  Is  necessary  because  plants  and  animals 
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exhibit  an  internal  free-running  rhythm  that  is  circadian,  approximately 
24  hours.  Since  the  power  transfer  capabilities  of  energy-selective 
behavior  require  an  internal  free-running  oscillation  and  a  comparable 
external  oscillation  the  Kaituna  model  was  modified  to  accomodate  that 
requirement.  The  annual  rhythm  is  supported  by  an  underlying  or 
fundamental  diurnal  rhythm.  Including  a  diurnal  rhythm  into  the 
original  solar  source  required  that  the  amplitude  and  mean  level  of  the 
original  rhythm  be  divided  by  0.3535  (or  increased  2.857  times)  in  order 
to  provide  the  same  amount  of  energy  as  the  original  solar  input  energy. 

Figure  47  presents  a  series  of  diagrams  which  explain  the  neccesity 
for  an  increase  in  the  amplitude  of  the  new  more  complex  rhythm  in  order 
to  maintain  the  same  energy  input.  Essentially  the  addition  of  a  period 
of  darkness  into  the  new  rhythm  requires  that  the  source  amplitude  be 
increased  to  account  for  the  loss  in  energy  during  darkness.  The  energy 
of  the  new  more  complex  source  rhythm,  at  the  same  amp  I  itude,  is  only 
35.35?  of  the  energy  of  the  original  solar  source — without  the  addition 
of  the  diurnal  rhythm. 

The   system  period  relationships  are  presented  'n   Figure  48   and 

indicate  that  the  source  contains  or  supplies  energy  to  the  system  at 

both  a  diurnal  and  annual  rhythm.   The  Spart?  na  producer  however,  sti I  I 

contains  or  produces  no  intrinsic  rhythmic  pattern.   Figure  49  shows  the 

simulation   plots    for   the  three   storages   within   the   model. 

Interestingly,   the   standing  Spart Ina  biomass   and   the   dead  matter 

decreased   as  a  result  of  the  addition  of   a   diurnal   rhythm.   Their 

-2  -2 

storages  reach  a  maximum  of  980  g  C  m   and  8250  g  C  m   ,   respectively. 

Nutrient  storage  of  labile  nitrogen  remained  relatively  constant  during 

the  simulation  as  in  the  previous  simulation,   reaching  a  maximum  of  975 


Figure  47:    Conversion   of   average  yearly  energy   input   for   diurnal 

energy  input:  (a)  mean  energy  level,  (b)  70.7  perc  ^'^  O' 
the  mean  energy  level  In  (a),  (c)  35.35  percent  of  the 
mean  energy  level  in  part  (a).  To  Insure  that  the  same 
energy  Is  input  by  a  diurnal  rhythm  the  mean  energy  level 
must  be  increased  by  2.8289  times. 
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F  i  gure   48 : 


The  system  period  relationships  for  the  original  Ka'tuna 
marsh  model  modified  to  include  a  diurnal  rhythm  for  the 
available  source  energy.  The  systam  has  a  annual  3nd 
diurnal  source  rhythm  and  the  Spart f  na  producer  has  no 
internal  rhythmic  pattern.  The  system  is  not  energy- 
selective  and  the  producer  does  not  sense  or  time  'ts 
environmental  energy  source. 
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Figure  49: 


Simulation  plots  for  the  Kaituna  salt  marsh  model  with 
the  addition  of  a  diurnal  solar  input.  Plot  (a)  Dead 
organic  matter.  Plot  (b)  SpartI na  blomass  and  system 
nutrient  storage.  Plot  (c)  Spartlna  unit  power.  Plot  Cd) 
Power  transfer  for  the  Spart ina  producer.  Dead  matter 
increases  but  not  as  rapidly  as  in  the  original  model. 
SpartI na  biomass  increases  slight  but  again  not  as  much  as 
in  the  original  model.  Nutrient  recycle  is  slower  than  in 
the  original  model,  and  unit  power  and  power  transfer  are 
a  I  so  decreased . 
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g  N  m   .   The   storage  of  dead  matter   shows  the  greatest  change 

decreasing  ]0%  from  its  previous  level.  One  explanation  for  the  model's 
reaction  to  the  addition  of  the  new  more  rapidly  changing  rhythm  in 
solar  input  could  be  that  the  Spartlna  (as  modelled  and  without  the 
guidance  of  some  internal  oehavioral  characteristic)  acts  as  a  filter  to 
rapidly  cycling  energy.  Therefore,  the  solar  input  energy  does  not 
transfer  as  rapidly  to  the  storage  of  sediment  carbon  as  in  the  previous 
simu I ation . 

All   of   the  following  models  of  the  Kaituna  marsh   'ncorporate   an 
autocata I yt ic  inertial  oscillator  into  the  Spartina  producer. 

Kaftuna  Salt  Marsh  with  the  Addition  of  an  Autocata I yt ic  Oscillator  as 
an  Internal  Behavioral  Guide  for  the  Spartina  Producer 

The   intent  of  this  simulation  was  to  introduce  an   autocatal yt ic 

i  nerti  a  I   osc  i I  I ator  i  nto  the  Sparti  na  producer  within  the  Kaituna  marsh 

model.   Figure  50a  indicates  the  relative  period  relationships  for   the 

modified  Kaituna  model.   The  Spartina  producer  has  an  intrinsic  diurnal 

rhythmic   pattern   and   as  a  result  links  to  that  rhythm  of   the  source 

energy.   The  relative  energy-acceptance  niche  for  the  Sparti  na  producer 

is   shown   in   Figure  50b,   and  is  only  intended  to  provide  a  visual 

representation  of  the  energy-acceptance  n'che.   Figure  51  is  an   energy 

diagram  of  the  Kaituna  salt  marsh  mini-model  including  the  addition  of 

an   autocata I ytic   inertial   oscillator  within  the  Sparti  na  production 

module.    Table  6   indicates  the  component   initial   conditions, 

description,   and  the  equations  for  the  flows  and  storages  in  Figure  40. 

The  biological  oscillator  produces  an  Internal   free-running  circadian 

rhythmic   behavioral   guide   for  the  Sparti na .   The  Intrinsic  rhythmic 

pattern  of  the  oscillator  guides  the  behavior  of  the  Spartina  during  the 


a 


Figure  50:  The  system  relative  niche  and  period  relationships  for  the 
modified  Kaituna  marsh  model  including  one  energy- 
selective  Sparti  na  producer.  The  systam  has  a  annual  and 
diurnal  source  rhythm  and  the  Sparti  na  producer  has  an 
diurnal  internal  rhythmic  pattern.  The  system  is  energy- 
selective  and  the  producer  does  sense  or  time  its 
environmental  energy  source.  Diagram  (a)  shows  the  system 
period  relationships  and  (b)  shows  the  relative  energy- 
acceptance  niche  for  the  producer. 
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Figure  51:  Kaituna  marsh  model  modified  by  the  addition  of  an  energy- 
selective  Spart ina  producer.  The  Spart I na  producer 
senses  the  solar  input  and  interacts  with  that  energy 
source  based  upon  its  internal  timing  ability  or  behavior 
to  acquire  energy.  The  following  plots  are  a  series  of 
simulations  using  this  particular  modification  of  the 
original  Kaituna  marsh  model.  The  internal  oscillator 
represents  the  action  of  a  biological  clock. 
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photosynthetic  production  of  biomass  using  labile  nitrogen  as  a  nutrient 

pool.   No  other   biological   oscillators  are  included   tor  the  model 

simulation,   and  all  forcing  function  rhythms  are  the  same  ones  used  for 

the  previous  simulations  of  the  marsh. 

Figures  52a   and   52b   present  the   levels   for   the  three  marsh 

storages.   Spartina  standing  biomass  increases  significantly  1150 
-2  -2 

g  C  m   when  compared  to  the  previous  simulation  980  g  C  m   ,   and 

compares  more   favorable  to  the  level  of   biomass   for   the  original 

simulation   of   the  mini-model.    Sediment  carbon   also    increases 

-2 
significantly  9200  g  C  m   when  compared  to  the  previous  simulation  8250 

-2 
g  C  m   ,   and  also  compares  well  with  the  original  level  for  the  storage 

dead  matter.   Labile  nitrogen  remains  relatively  constant  decreasing  to 

-2  -2 

a  level  of  900  g  N  m   as  compared  to  1000  g  N  m   .   Interestingly,   the 

addition  of  the  autocata I ytic  inertial  oscillator,  (or  the  addition  of 
an  intrinsic  rhythmic  behavior)  to  the  new  mini-model  of  the  salt  marsh 
compensates  for  the  decreases  in  biomass  and  dead  matter  resulting  from 
the  addition  of  a  diurnal  solar  rhythm  to  the  original  model.  The 
result  can  only  be  explained  as  a  better  coupling  between  the  system  and 
the  solar  input.  Since  the  model  with  the  new  solar  Input  (the  day- 
night  pattern  supporting  an  annual  rhythm)  and  no  autocatal yt ic  inertial 
oscillator  did  not  obtain  as  much  solar  energy  as  did  the  model  with  the 
Inertial  oscillator;  the  addition  of  the  inertial  oscillator  must 
produce  a  better  couple  between  source  and  system.  This  result  Is  a 
significant  indication  of  the  ability  of  energy-selective  units  to 
better  couple  with  sources  that  supply  energies  which  are  cyclic  or  are 
cyclic  about  some  mean  energy  level. 


Figure  52:  Simulation  plots  for  the  modified  Kaituna  marsh  model. 
Plot  (a)  Dead  organic  matter.  Plot  (b)  Sparti  na  biomass 
and  system  nutrients.  Plot  (c)  Sparti  na  unit  power.  Plot 
(d)  Sparti na  power  transfer.  Dead  matter,  and  biomass 
increased  to  approximately  the  same  level  In  the  original 
model  simulation.  The  storage  of  system  nutrients 
increased,  however  not  to  the  level  predicted  by  the 
original  model.  Therefore  the  recycle  of  nutrient  Is 
slower  than  for  the  original  simulation.  The  Sparti  na 
unit  power  and  power  transfer  both  exceeded  the  levels  for 
the  original  simulation  Indicating  the  advantage  offered 
because  of  energy  selective  behavior. 
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Figures  52c   and  52d  show  the  unit  power   requirement  and   power 

transfer   capability  of  the  Spartfna  producer.   The  unit  power  of   the 

Spartina,   in  solar  Calories,  acheives  an  initial   maximum  power  of  4500 

-2     -2  -2     -2 

Cal   m    day   thereafter  leveling  at  3400  Cal  m    day   .   The  power 

transfer  of   the  Sparti  na  producer  increases  gradually  to  a  maximum 

-2    -2 
sustained  level  of  18  Cal  m   day   .    The  increase  in  unit  power  by  the 

addition  of   the  autocata I yt ic  inertial  oscillator  is  ]5%      above  the 

original   models  prediction,   and  the  increase  in  power  transfer  by  the 

Spartina   is  6%      above  the  original   model.   The   simulation   clearly 

indicates   that  energy  selectivity  produces  a  better  couple  between   the 

producer  and  the  source,   and  as  a  result  passes  on  a  greater  amount  of 

power  to  the  system. 

Kaituna  Salt  Marsh-Min i-model  with  Spartina  Production  Module  that  Uses 
an  Autocata I ytic  Oscillator  Producing  a  Niche  Incapable  of  Utilizing 
the  Solar  Input  Energy 

This  simulation  shows  the  enregy-sel ect I ve  property  of  the  Sparti na 

producer   as  a  result  of  the  addition  of   the  autocata; ytic   inertial 

oscillator.   The   system   relative  period  relationships  are  shown   in 

Figure  53a.   The  energy-acceptance  niche  for  the  producer  is  shown   in 

Figure  53b.   The   figure  also  indicates  that  the  two  rhythms  of   the 

source  are  not  at  the  same  period  as  is  the  Intrlnsci  oscillation  of  the 

producer.   Figure  54  shows  the  effect  on  the  salt  marsh's  storages  as  a 

result  of   changing   the  Spartlna's   internal   free-running  behavioral 

characteristics.   The  model   is  Identical  to  the  previous   simulation 

model   with   the  exception  of  the  period  for  the  autocata I yt Ic   inertial 

oscillator's   Internal   rhythm.   The  oscillator's   rhythm   during   the 

previous  simulation  linked  exactly  with  the  diurnal  rhythm  of  the  solar 


Figure  53:  The  system  relative  niche  and  period  relationships  for  the 
modified  Kaituna  marsh  model  including  one  energy- 
selective  Spart [ na  producer.  The  systam  has  a  annual  and 
diurnal  source  rhythm  and  the  Sparti  na  producer  has  an  20 
day  internal  rhythmic  pattern.  The  system  is  energy- 
selective  and  the  producer  does  sense  or  time  its 
environmental  energy  source.  The  producer  rhythm  employed 
for  this  simulation  does  not  allow  the  unit  to  acquire 
available  source  energy  and  therefore  the  unit  does  not 
survive.  Diagram  (a)  shows  the  system  period 
relationships  and  (b)  shows  the  relative  energy-acceptance 
niche  for  the  producer. 
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Figure  54:  Simulation  plots  for  SpartI na  producer  with  an  energy 
acceptance  niche  Incapable  of  utilizing  solar  energy.  Plot 
(a)  Dead  matter.  Plot  (b)  Spartlna  blomass.  Plot  (c) 
Unit  power  for  the  SpartI na  producer.  Plot  (d)  SpartI na 
power  transfer.  Dead  matter,  blomass,  unit  power,  and 
power  transfer  all  decrease  when  compared  with  the 
original  simulation.  Nutrient  storage  within  the  system 
Increases  because  the  consumers  recycle  dead  matter  Into 
nutrients  faster  than  dead  matter  Is  replaced.  The  system 
and  the  Individual  do  not  survive.  The  rejection  of  input 
energy  by  Spartlna  reduces  input  power  and  power  transfer 
thereby  eliminating  the  system. 
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energy  input,   however  the  new  inertial  oscillator  rhythm  is  20  times  as 

long  as  the  previous  rhythm. 

The  plots  In  Figure  54  clearly  indicate  that  if  Spartlna ' s  internal 

rhythm  were  elongated  to  20  days,   instead  of  1  day,   the  producer  could 

not  obtain  adequate  solar  energy  to  support  itself  or  the  remainder  of 

the   system.   Spart i  na   standing   biomass   rapidly  diminishes   and   is 

converted  to  sediment  carbon/dead  matter.   The  storage  of  dead  matter  in 

-2 
turn   rapidly   decreases   to  a  level  of  1000  g  C  m    and   is   in   turn 

converted  by  consumer  action  to  system  nutrients.   The  nutrient  nitrogen 

storage   increase  because  consumers  are  adding  to  its  storage  while  the 

Inadequate  behavioral  activity  of  the  producer  decreases  Its  demand   for 

the   uptake  of  nutrients.   As  a  consequency  the  system  stores   nitrogen 

-2 
and  the  storage  reaches  a  maximum  level  of  7000  g  N  m 

Figures   54c  and  54d  indicate  that  the  power  demand  by  the  Sparti  na 

-2 
decreases   as  a  result  of  an  inadequate  niche  to  a  level  of  700  Ca I   m 

-2 
day   .   Power   transfer  from  the  Spartlna  producer  to  the  remainder  of 

-2    -2 
the  system  decreases  to  a  level  of  2.2  Cal  m   day 

Within  the  model  the  system  relies  upon  the  Spartlna  producer  for 
Its  sustenance  and  therefore  the  system  vanishes.  In  reality  one  would 
expect  some  other  producer  to  fill  the  niche  vacated  by  the  Sparti na, 
but  the  model  does  not  incorporate  that  concept  into  the  simulation. 

The  following  models  address  the  replacement  concept  via 
consideration  of  energy-selective  competition  within  the  Kaltuna  web. 

Kaltuna  Salt  l^arsh  Mini-model  with  Competitive  Exclusion  of  a  Producer 
with  an  Inadequate  Niche 

This   simulation  introduces  energy-selective  competition  into  the 

Kaltuna  marsh  web  with  the  intent  of  analyzing  this  concept  in  a  web 

system.   The   system  period  relationships  are  presented  in  Figure  55a. 


Figure  55:  The  system  relative  niche  and  period  relationships  for  the 
modified  Kaituna  marsh  model  including  two  energy- 
selective  producers  in  a  competition  for  available  solar 
energy.  The  systam  has  a  annual  and  diurnal  source  rhythm 
and  the  producers  have  a  diurnal  and  20  day  internal 
rhythmic  pattern,  respectively.  The  system  is  energy- 
selective  and  the  producers  do  sense  or  time  the 
environmental  energy  source.  The  short  period  producer's 
rhythm  links  with  the  avaiable  energy  rhythm  and  as  a 
result  the  short  period  producer  is  more  competitive  than 
the  long  period  producer.  The  competition  proceeds  to 
exclusion  of  the  long  period  producer.  Diagram  (a)  shows 
the  relative  system  periods  (b)  shows  the  relative  energy- 
acceptance  niches. 
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The  long  rhythmed  producer  clearly  generates  a  rhythm  that  is  not  linked 

to  either  source  rhythm,  while  the  short  rhythmed  producer  links  to  the 
diurnal  rhythm  of  the  source.  Figure  55b  shows  the  energy-acceptance 
niches   for   both  producers. 

Figure  56  is  an  energy  diagram  of  the  competition  modified  mini- 
model  for  the  Kaituna  salt  marsh.  The  model  consists  of  two  competing 
producers.  The  only  difference  between  the  two  competitors  is  their 
internal  rhythmic  behavior.  Spartlna  townsendii,  unit  #1,  has  an 
internal  circadian  period  for  a  behavioral  rhythm.  Producer  1^2,  a 
substitute  unit  for  investigative  purposes,  has  an  internal  behavioral 
rhythm  with  a  period  of  20  days.  No  Input  energy  with  a  period  of  20 
days  is  available  to  the  system.  The  solar  input  supply  energy  has  a 
diurnal  support  rhythm  which  is  supplied  to  both  producers.  The 
internal  rhythm  of  Spartlna  provides  it  with  a  better  linkage  to  the 
solar  source  energy  and  is  therefore  used  by  the  Sparti  na  to  exclude  its 
compet i  t ion . 

Figures   57a  through  57d  clearly  show  the  competitive  advantage  of 

the  Spartlna' s  behavior.   While  Its  competitor  is  only  obtaining  solar 

energy   for  small  portions  of  the  day  Spartlna' s  behavior  allows   it  to 

capture  solar   energy  during  the  entire  daylight  hours.   As  a  result, 

-2 
Spartlna' s   standing  blomass  level  Increases  to  a  mean  of  1000  g  C  m 

while   its  competitor's  blomass  level  diminishes.   The  maximum  blomass 

-2 
level  obtained  by  the  Spartlna  producer  is  1150  g  C  m    and  this  occurs 

during  the   fourth   year  of  the   simulation.   After   the   fourth   year 

Sparti na' s  blomass  follows  an  annual  pattern  about  Its  mean  level  until 

the  end  of  the  10  year  simulation.   The  sediment  carbon/dead  matter   in 

-2 
the  system  increases  reaching  a  maximum  level  of  9500  g  C  m 


Figure  56:  Kaituna  marsh  model  modified  to  include  competitive 
producer.  The  system  blomass  is  evenly  split  between  the 
competitors . 
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Figure  57:  Plots  for  modified  Kaituna  marsh  model  with  competitive 
exclusion  of  a  producer  with  an  inadequate  energy 
acceptance  niche.  Plot  (a)  sediment  carbon.  Plot  (b) 
biomass  and  nutrients.  Plot  (c)  unit  power  for 
competitors.  Plot  (d)  power  transfer  for  the  competitors. 
The  plots  clearly  indicate  the  competitive  advantage 
afforded  the  producer  with  the  acceptance  niche  period 
closest  to  the  source  energy's  period. 


:a) 


ic,:co 


182 


Time    (Years) 


I'j 


2,000 


:b) 


ime  '  "ear'; 


;.oco 


10 


Time   (rears) 


183 


■a 

c 

o 


01 

s- 
<u 

c 


03 


4-> 

q; 

■a 

s_ 

0 

0 

5^ 

"4- 

c 

on 

0 

c 

0 

-4-> 

•  r- 

•  P— 

-M 

4-> 

ro 

o; 

3 

Q- 

cr 

t; 

UJ 

0 

0 

TD 

c 

-C 

(O 

1/1 

i. 

•» 

fD 

c 

Z 

0 

•^ 

rtj 

■(-> 

c 

a. 

3 

•fm- 

■t-J 

i~ 

•  (— 

0 

ra 

I/) 

^ 

a; 

0 

a; 

-C 

#1 

-t-j 

a; 

(C     ■!- 


to 

* 

ai 

c 

cr. 

0 

03 

i. 

■!-> 

0 

re 

-M 

3 

OO 

to 

> 

O) 

re 


^i/ics  —  «  ••■        3        - 

■«  ^     ■•  ^         «-i    ^     ^  •— '    •'^ 

«  o       1  2     ^      .^£     ■      u;     uo 

'— •      ~«  —  -*      _it—  ■*■      * 

C                     :%i  -M.         >  *    ^     ui^      u^ 

^  w        _ 

*         3        u 

c  -^    -: 


— -      —      rs, 


•3 
O 


O       E 


—        u        O 


c       o       ^      -^ 
i^     1=     —     3 


■^ 

•^ 

c 

-D 

v> 

c. 

- 

^ 

>. 

0 

^ 

X 

3 

^ 

3 

u 

0 

c 

Oi 

CJ 

C 

— 

t 

K 

^ 

C 

u 

"■. 

f. 

0 

0 

C 

0 

c 

•-I 

^ 

■— 

^ 

.J 

^ 

^ 

c 

■•— 

gn 

t/* 

1/1 

/B 

'fl 

t/* 

Ik 

z 

— 

^ 

CJ 

Lrt 

>/i 

u^ 

3 

E 

^ 

^ 

c< 

U 

ro 

m 

-^ 

0 

OJ 

Of 

fO 

E 

^ 

c= 

^ 

■c 

> 

c 

L. 

u 

■^ 

> 

-a 

^— 

C 

0 

3 

fc 

(^ 

^ 

kJ 

^ 

c 

0 

a; 

0 

0 

-1 

a: 

►— 

t/1 

CC 

c: 

c: 

LJ 

Z 

z 

z 

VO 

>c 

n^ 

^ 

in 

0 

C^ 

CE 

■^ 

^ 

CN. 

f-^ 

(\, 

w 

tr 

c 

0 

■n 

0 

\o 

CN- 

tsj 

m 

CN( 

(NJ 

tO 

un 

^ 

^ 

^ 

rsi 

•— 

0 

^ 

w 

1/1 


LT) 

O) 

cr. 


O 


(U 


J2 


o 


184 

Comparison  of  the  simulation  plots  indicates  that  the  competition 
between  the  Spart ! na  and  its  longer  rhythmed  competitor  is  relatively 
none  existent.  In  tact  the  system  as  a  whole  reached  approximately  the 
same  levels  for  standing  biomass,  dead  matter,  and  nutrients  as  is 
acheived  without  competition.  This  result  may  be  explained  by  the 
reluctance  of  the  long  rhythmed  producer's  to  provide  any  real 
competition  with  the  Spart I na.  Since  the  longer  rhythmed  producer  could 
or  would  not  capture  solar  input  with  the  same  efficiency  that  the 
Spart ina  captures  the  energy  the  competition  ends  quickly  and  produces 
little  or  no  stress  on  the  Spartina. 

Table  7  shows  the  distribution  of  biomass  and  flows  for  the 
addition  of  the  competing  producer.  The  system  biomass  and  solar  input 
energy  is  shared  between  the  producers  to  start  the  simulation.  Pathway 
flows  to  the  storages  of  sediment  carbon  and  labile  nitrogen  are  halved 
to  maintain  the  same  system  input  to  the  storages.  Therefore,  the  model 
incorporates  a  new  competitive  unit  without  disturbing  the  initial  sense 
of  the  system. 

Kaltuna  Salt  Marsh  i^odifled  Mini-model  to  Include  Competition  Between 
Producers  with  the  Ability  to  Obtain  Solar  Energy  at  Different  Rhythms 

The  competition  model  used  for  this  simulation  is  the  same  model 

used   for   the  previous  simulation,   however  the   period  of   the   long 

rhythmed   competitor   (producer  #2)  is  increased   to  365   days.   This 

simulation   was  employed   to  investigates  the  effects  of   competition 

between  two  competitors  that  utilize  different  fundamental  rhythms   from 

the  same  source.    In  the  competition  section  that  employed  hypothetical 

producers  this  concept  was  analyzed  an  indicated  that  the  short   rhythm 

competitor   should  outcompete  the   long   rhythmed   competitor.    This 


135 
simulation  expands  that  concept  into  web  systems  in  an  attempt  to 
discover  if  web  feedback  has  an  effect  on  the  compet'ton.  Figure  58a 
shows  the  system  period  relationships  for  the  model  simulated.  Source 
energy  is  composed  of  two  rhythms  one  a  circad'an  rhythm  and  the  other  a 
circannual  rhythm.  The  acceptance  niches  for  the  two  competitors  are 
shown  in  Figure  58b. 

Figures  59a  and  59b  show  the  levels  for  the  four  storages  of  the 
salt  marsh  (Spartlna  biomass,  labile  nitrogen,  sediment  carbon,  and  the 
storage  of  biomass  for  the  longer  rhythmed  producer).  Initially  the 
Spart I na   with   the  circadian  internal  rhythm  outcompetes  the  unit   with 

the   internal  annual  behavior  pattern,   and  reaches  a  maximum   standing 

-2 
biomass   of  1050  g  C  m   at  time  equal  3  years  in  the   simulation.    The 

long   rhythmed   (365   day)   competitor  outcompetes   the   short   rhythmed 

producer   during  the  middle  years  of  the  simulation  reach'ng   a  max"mum 

-2 
biomass   level   of  300  g  C  m    at  time  equal  6  years  'n  the   simulation. 

However,  the  Spart  i  na  again  regains  its  dominant  posture  at  time  equal  8 

years.     The   simulation   indicates   that   the   two   producers   coexist 

thereafter  but  with  slightly  less  total  system  b'omass.    The   summation 

of   the  two  producers'  biomass  indicates  that  system  biomass  level  at   a 

-2 
mean   of   1000   g  C  m   .   Dead  matter  within  the   system   increases   and 

-2 
reaches   a  maximum   level  of  9800  g  C  m  — approximately  equal   to   the 

level  obtained  during  the  previous  simulation.    Nutrient  storage  within 

-2 
the   system   remained   relatively   constant  at  950  g  N  m    as   in   the 

previous  five  simulations.    The  competition  produces  a  system  in   which 

the  short  rhythmed  Spartina  townsendff  dominates  for  the  majority  of  the 

time.   The  system  has  more  diversity  while  maintaining  approximately  the 

same  level  of  dead  matter  and  nutrient  storage. 


Figure  58:  The  system  relative  niche  and  period  relationships  for  the 
modified  Kaituna  marsh  model  including  two  energy- 
selective  producers  in  a  competition  for  available  solar 
energy.  The  systam  has  a  annual  and  diurnal  source  rhythm 
and  the  producers  have  a  diurnal  and  annual  internal 
rhythmic  pattern,  respectively.  The  system  is  energy- 
selective  and  the  producers  do  sense  or  time  the 
environmental  energy  source.  The  short  period  producer's 
rhythm  links  with  the  avaiable  energy  rhythm  as  does  the 
long  period  producer.  The  web  feedback  effect  provides 
long  period  feedback  from  the  system  and  as  a  result  the 
long  period  producer  is  more  competitive  and  coexistance 
is  established  between  the  two  producers.  Diagram  (a) 
shows  the  relative  system  periods  (b)  shows  the  relative 
energy-acceptance  niches. 
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Figures  59c   and  59d   show  the  power   requirements   for   the   two 

competitors,   as  well  as,   their  power  transfer  to  the  remainder  of  the 

system.   The   unit  power  for  the  two  competitors   indicates   that  the 

Spart i na  producers  acquires  the  greatest  amount  of  Input  power  reaching 

-2    -2 
a   steady  state  level  of  2500  Cal  m   day   .   The  plot   also   Indicates 

that  the  total   system   input  power  from  the   solar   source   remains 

-2    -2 
approximately  the  same  at  3700  Cal  m   day 

Summation  of  the  power  transfer  for  both  producers  indicates   that 

-2     -2 
the  total  power  transfer  to  the  rest  of  the  system  Is  19  Cal  m    day   , 

and  that  the  Spart 1 na  producer  transfers  the  largest  amount  of  the  two 

un  i  ts . 

Remembering  that  the  competition  simulation  plots  for  energies  from 
the  same  quality  source  in  the  previous  section  utilized  no  feedback 
from  other  components  in  the  system  provides  Insight  Into  the 
competitive  coexistance  that  Is  shown  In  the  plots  In  Figure  45.  Since 
the  longer  period  annual  rhythm  of  the  input  solar  energy  Is  from  the 
same  energy  quality  source  as  is  the  diurnal  rhythmed  energy,  and  since 
the  longer  period  producer  coexists  with  the  Spart I na  producer;  It 
appears  that  long  period  competitive  success  at  the  bottom  of  the  food 
chain  requires  the  feedback  of  a  longer  period  rhythm  from  a  component 
within  the  system — or  for  that  matter  from  a  component  within  the  next 
higher  system.  The  annual  rhythm  of  the  nutrient  storage  provides  that 
longer  rhythm  feedback  for  this  simulation.  The  nutrient  storage 
acquires  an  annual  rhythm  via  the  nutrient  input  from  the  river  source 
and  from  the  Input  from  the  Spart I na  producer. 

Odum  1983,  indicates  that  there  is  an  advantage  to  control 
provided  by  long-period  feedback  from  components  In  or  out  of  the  system 


Figure  59:  Plots  of  competition  between  producers  in  the  Kaituna 
marsh  that  obtain  energies  at  different  frequencies  from 
the  same  energy  quality  source.  Plot  (a)  dead  matter. 
Plot  (b)  biomass  and  system  nutrients.  Plot  (c)  unit 
power  for  the  two  competitors.  Plot  (d)  power  transfer 
for  the  competitors.  The  plot  Indicates  a  coexistance 
between  competitors. 
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in  question.  The  simulation  plots  shown  in  Figure  45  clearly  indicate 
the  advantage  of  longer  period  feedback  in  this  simulation  to  be 
increased  system  diversity. 

Kaituna  Salt  Marsh  Mini-model  with  Competition  Between  an  Energy- 
Selective  Producer  and  a  Non-energy  Selective  Producer 

The  intent  of  this  simulation  was  to  analyze  the  competitve 
strength  of  an  energy-selective  producer  when  competing  with  a  non- 
energy-selective  producer.  The  system  period  relationships  are  shown  in 
Figure  60a.  The  energy-acceptance  niche  of  the  energy-selective 
producer  is  shown  in  Figure  60b.  The  niche  for  the  non-energy-selective 
producer  is  not  shown  because  it  accepts  energy  at  all  frequencies  or 
per  iods . 

Figure  61  is  an  energy  circuit  diagram  of  the  Kaituna  marsh  model 
with  the  addition  of  a  non-energy-selective  competitor.  Figures  62a 
through  62d  present  the  simulation  results.  The  biomass  plots  for  the 
two  producers  clearly  indicate  the  competitive  advantage  offered  by 
energy-selective  behavior.  The  energy-selective  SpartI na  unit 
competitively  excludes  the  non-energy-selective  unit  rapidly   and  with 

little  or   no  difficulty.    Biomass  for  the  energy-selective  producer 

-2 
reaches   a  maximum  level  of  1100  g  C  m   while  the  biomass  of   the  non- 
energy-selective   producer   decreases   to  a  level   of   exclusion.   Dead 

-2 
matter  within  the  system  reaches  a  level  of  9500  g  C  m   as  i n  many  of 

the  previous  simulations,   and  the  storage  of  nutrients  reaches  a   level 

-2 
of   1000  g  n  m   .    The  producer  unit  and  transfer  power  plots  clearly 

indicate  the  dominance  of  the  energy-selective  SpartI na  producer.  The 
plot  of  unit  power  for  the  two  competitors  indicates  that  the  energy- 
selective  producer  quickly  captures  the  available  energy  of  the  source, 


Figure  60:  The  system  relative  niche  and  period  relationships  for  the 
modified  Kaltuna  marsh  model  including  one  energy- 
selective  and  one  non-energy-selective  producer.  The 
system  has  a  annual  and  diurnal  source  rhythm  and  one 
producer  has  a  diurnal  internal  rhythmic  pattern.  The 
system  is  energy-selective  and  one  producer  senses  or 
times  the  environmental  energy  source.  The  energy- 
selective  producer  is  the  better  competitor.  Diagram  (a) 
shows  the  relative  system  periods  (b)  shows  the  relative 
energy-acceptance  niches.  The  non-energy-selective 
producer  acquires  energy  at  all  periods  but  not  as 
effectively  as  the  energy-selective  producer  in  the  range 
or  niche  of  the  energy-selective  unit. 
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Figure  61:  Energy  circuit  model  for  competition  within  the  Kaituna 
marsh  between  an  energy  selective  producer  and  a  non- 
energy  selective  producer. 
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Figure  62:  Plots  for  competition  between  an  energy-selective  producer 
and  a  non-energy-sel ecti ve  producer.  Plot  (a)  sediment 
carbon.  Plot  (b)  biomass  and  system  nutrients.  Plot  (c) 
unit  power.  Plot  (d)  power  transfer  for  the  competitors. 
Clearly  the  energy  selective  producer  prevails. 
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and   as   a  result  obtains  the  most  source  power — capturing  a  maximum  of 

-2     -2 
4500  Ca I   m    day   .   The  producer   incapable  of   rhythmic   behavior 

captures  virtually  no  input  power  and  is  excluded  within  the  first  month 

of   the  competition.   The   plot  of  power  transfer   for   the  producers 

indicates  that  the  energy  selective  producer  transfers  the  most  power  to 

-2 
the  remainder  of  the  marsh  food  chain,   reaching  a  level  of  18  Cal   m 

-2 
day 

Table  8   summarizes   the   system   initial    conditions,    pathway 

coefficients,  and  pathway  coefficient  equations  for  the  simulation.   The 

table  provides  a  description  of  the  flow  or  storage,   and  the  equations 

for  the  model . 

Kaituna  Salt  Marsh  Mini-model  with  a  Two-OscM  I ator  Producer  Competing 
for  a  Limited  Energy  Flow  Composed  of  the  Addition  of  Two  Rhythms 

The  intent  of  this  section  is  to  show  the  effects  on  the  Kaituna 
web  from  the  addition  of  an  Spart I na  producer  with  two  intertial 
oscillators.  Figure  63a  and  63b  show  the  system  period  relationships 
and  energy-acceptance  niche  for  the  producer. 

Figures  65a  through  65d  show  the  results  of  the  Kaituna  mini-model 
simulation  which  contains  the  complex  energy-selective  producer  in 
Figure  64.  The  producer  uses  an  Internal  behavior  composed  of  and 
guided  by  the  Interaction  of  two  biological  Inertlal  oscillators.  The 
first  autocata I yt ic  inertlal  oscillator  produces  a  circadian  rhythmic 
behavioral  pattern,  while  the  second  autocata I ytl c  inertlal  oscillator 
produces  a  circannual  rhythmic  behavior.  The  combination  of  the  two 
inertlal  oscillators  provides  the  producer  with  a  behavioral  pattern 
sensitive  to  both  of  the  rhythms  of  the  source.  The  result,  as 
indicated  from  the  plots  In  the  figures,  provide  some  Interesting  facts. 


Figure  63:  The  system  relative  niche  and  period  relationships  for  the 
modified  Kaituna  marsh  model  including  one  energy- 
selective  producer.  The  systam  has  a  annual  and  diurnal 
source  rhythm  and  one  producer  with  a  diurnal  and  annual 
internal  rhythmic  pattern.  The  system  is  energy-selective 
and  the  producer  senses  or  times  the  environmental  energy 
source.  The  energy-selective  producer  generates  a  niche 
with  two  areas  of  energy-acceptance.  Diagram  (a)  shows  the 
relative  system  periods  (b)  shows  the  relative  energy- 
acceptance  niches.  The  energy-selective  producer  acquires 
energy  at  all  periods. 
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The   addition  of  the  second  oscillator  does  not  significantly  effect  the 

generation  of  producer  biomass,   only  slightly  increasing  it  to  1200  g  C 

-2 
m   .  The  system  receives   no  more  input  solar   power   from   the  two 

oscillator   producer   than   it  did  from  the  same  producer   guided   by  a 

single  internal  biological  oscillator. 

Spart Ina  unit  power,  Figure  65c,  is  the  same  for  the  complex 
producer  as  it  was  for  the  single  inertial  oscillator  producer;  and  the 
power  transfer  by  the  two  inertial  oscillator  producer  to  the  remainder 
of  the  system.  Figure  65d,  is  the  same  as  for  the  single  inertial 
oscillator  or  simple  Spartfna  producer.  That  does  not  however  indicate 
that  the  producer  is  wasting  time  and  presumably  energy  in  developing 
itself  as  a  more  complex  unit.  The  fact  is  the  unit  is  a  general ist. 
With  the  development  or  addition  of  the  new  inertial  oscillator  the  unit 
adds  a  new  dimension  to  its  energy  selective  behavior,  and  benefits  from 
the  new  behavioral  additions  when  placed  into  a  competitive  environment. 
The  results  of  the  simulation  indicate  only  that  the  addition  of  more 
complex  behavior  does  not  hinder  the  producer's  normal  performance. 
However,  the  two  complex  Spart I na  producer  does  reach  steady  state 
conditions  more  rapidly,  and  this  aspect  should  aid  in  competition  with 
other  producers. 

Tables  9  presents  the  initial  condition  for  the  model  storages, 
flow  coefficients,  and  pathway  equations.  No  pathway  flows  were  changed 
during  the  simulation  of  this  model.  The  only  addition  to  the  model  Is  a 
second  autocata I yt Ic  Inertial  oscillator  within  the  Spartlna  production 
modu le. 


Figure   64:   Energy  diagram  of  a  two-clock  producer  positioned  'n 
Kaituna  marsh  web. 
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Figure  65:  Plots  for  the  two-clock  SpartI na  producer  located  in  the 
Kaituna  marsh.  Plot  (a)  sediment  carbon.  Plot  (b) 
biomass  and  system  nutrients.  Plot  (c)  SpartI na  unit 
power.  Plot  (d)  SpartI na  power  transfer.  The  simulation 
plots  indicate  that  the  addition  of  a  second  clock 
produces  no  change  in  system  state  variables.  This 
finding  coincides  with  the  earlier  competition  models  in 
the  previous  section.  However,  earlier  models  indicate 
that  the  producer  is  a  more  competitive  unit  than   is   a 
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DISCUSSION 

The  following  discussion  is  presented  in  three  sections.  First, 
the  discussion  concentrates  on  the  autocata I yt ic  inertial  oscillator 
developed  during  this  research.  Mext,  the  discussion  continues  with 
energy-selective  behavior  and  its  effects  on  niche,  diversity, 
competition,  and  maximum  power.  The  previous  effects  are  presented  with 
the  discussion  concerning  maximum  power  interwoven  throughout  the 
sections  devoted  to  niche,  diversity,  and  competition.  Third,  future 
research  topics  are  described,  including  the  possibility  that  energy 
selectivity  is  a  general  systems  principle,  and  what  future  research 
might  be  undertaken  in  order  to  prove  or  disprove  that  possibility.  The 
possibility  of  the  relevance  of  the  autocatal ytic  inertial  oscillator  as 
a  physiological  clock  is  also  discussed. 

The  Autocatal ytic  Inertial  Oscillator 

Autocatal ytic  inertial  oscillators  provide  systems  with  a  means  of 
selecting  energy  and  power  inputs.  The  use  of  inertial  oscillators  by 
systems  clearly  produces  an  energy-selective  activity  for  that  system, 
and  as  a  result,  produces  power  transfer  for  a  selected  range  of  input 
energ  ies . 

In  the  previous  section  autocatal ytic  inertial  oscillators  were 
shown  to  aid  in  the  development  of  a  system  energy-acceptance  nicne. 
The  system's  energy-acceptance  niche  being  described  or  defined  by  the 
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range  of  Input  energies  a  system  utilizes  for  the  particular  internal 
rhythmic  pattern  it  possess  or  generates.  This  research  has  shown  that 
an  energy-acceptance  niche  for  a  system  may  be  defined  by  the  generation 
of  an  internal  autocata I yt ic  inertial  oscillation.  When  the  period  or 
frequency  of  the  available  energy  is  to  far  removed  from  that  of  the 
system's  internal  autocatal yt ic  inertial  oscillation  then  energy  is 
rejected  and  the  system  must  utiize  other  sources  if  it  is  to  survive. 

Simulations  in  the  previous  section  also  indicate  that  the 
autocata I ytic  inertial  oscillator  plays  a  role  in  defining  the 
competitive  advantage  one  system  has  over  another  system.  The  system 
that  most  closely  duplicates  the  period  of  the  available  source  energy 
survives  in  competition  with  other  single  rhythmed  systems. 
Furthermore,  energy-selective  behavior  is  advantageous  even  when 
competing  with  a  non-energy-selective  system.  If  the  available  energy 
supply  is  rhythmic  and  within  the  acceptance  range  of  an  energy- 
selective  system  that  system  effectively  outcompetes  a  non-energy- 
selective  system.  Competition  between  a  single  rhythmed  system  and  a 
more  complex  general  1st  system  leads  to  the  competitive  exclusion  of  the 
single  rhythmed  system. 

Simulations  have  also  indicated  that  competition  for  available 
energy  can  produce  coexistance  of  competiting  energy-selective  systems. 
The  overlapping  of  energy-acceptance  niches  generated  by  internal 
autocata I yt ic  inertial  oscillators  can  produce  both  competitive 
exclusion  and  coexistance.  When  the  rhythyms  generated  by  the  competing 
systems'  autocatal yt ic  inertial  oscillators  match  closely  by  period  and 
are  closely  matched  with  their  common  energy  source  the  result  Is 
competitive  coexistance.   When   the   Internal   autocatal yt Ic   Inertial 
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oscillations  of   the  competitors  differ  greatly  then   the  result  of 
competition  is  competitive  exclusion. 

External  and  Internal  Timing  of  System  Oscillations 

The  problem  of  external  or  internal  timing  is  important  to  the 
concept  of  energy-selectivity.  If  timing  is  generated  externally  then 
the  generation  of  energy-selective  behavior  is  not  likely  to  occur! 
Without  an  autonomous  Internal  rhythm  to  guide  the  rhythmic  activity  of 
a  system  the  system  merely  reflects  the  rhythm  of  its  environment.  The 
ability  to  be  energy-selective,  that  Is  to  reject  energy  because  it  is 
not  of  the  same  or  closely  related  period,  depends  upon  the  ability  to 
discern  a  different  period  of  fluctuation.  Without  the  ability  to 
generate  an  internal  free-running  rhythmic  activity  a  system  can  only 
reflect  the  rhythm  of  its  environment  and  therefore  cannot  discern  any 
difference  between  Its  rhythmic  activity  and  that  of  the  environment. 

Energy  Selectivity  and  Systems 

Simulation  plots  In  the  results  section  of  this  analysis  Indicate 
that  systems  which  have  the  ability  to  generate  internal  rhythmic 
patterns  compete  for  energy  within  a  given  range  (energy-acceptance 
niche)  of  energy  alternatives,  and  do  so  more  effectively  than  units 
that  cannot  act  to  selectively  choose  their  energy  Input.  Furthermore, 
the  range  of  the  system's  internal  rhythmic  behavior  cannot  vary  to  far 
from  the  frequency  of  Its  primary  energy  source.  If  a  system  generates 
a  rhythmic  Internal  guide  that  is  to  far  removed  from  that  of  its  energy 
source  the  system  will  not  survive.  Within  a  system's  energy-acceptance 
niche  survival  is  determined  by  the  strength  of  Its  competitiveness. 
Factors  effecting  the  competitive  strength  of  a   systems   include  the 
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strength  and  quality  of  the  feedback  available  to  the  system  from  other 
systems  and  components,  the  closeness  of  the  system's  internal  rhythm 
to  that  of  the  rhythm  of  its  primary  energy  source,  the  system's  power 
transfer  capabilities,  i.e.  its  participation  in  the  larger  system's 
maximization  of  power,  and  the  competitive  strength  of  the  systems  with 
which  it  must  compete. 

The  simulations  presented  in  this  study  indicate  that  energy  input 
is  most  efficiently  accomplished  when  the  system's  internal 
autocatal ytic  inertial  rhythmic  activity  is  at  the  same  period  as  its 
primary  energy  source.  Therefore,  the  greatest  power  Input  and  power 
transfer  is  accomplished  at  free-running  internal  rhythm  of  the  system 
when  that  rhythm  is  of  the  same  period  as  that  of  the  system's  primary 
source  energy.  All  source  energy  with  periods  outside  of  the  system's 
energy-acceptance  niche  are  rejected  and  no  power  is  transfered  by  the 
system.  For  source  energy  within  the  system's  energy  acceptance  niche 
but  not  at  the  same  frequency  or  period  as  the  system's  free-running 
oscillation  energy  and  power  are  accepted  but  not  at  the  most  efficient 
rate.  Thus  power  transfer  may  occur  in  energy-selective  systems  at 
other  than  the  maximum  power  transfer  capabilities  of  the  system. 

Transfer  of  power  at  other  than  the  most  efficient  rates  may 
reflect  the  necessity  for  increased  diversity.  A  larger  number  of 
components  transfering  power  at  lesser  rates  may  still  be  able  to 
transfer  the  same  amount  of  power  to  the  next  larger  system  as  that  of  a 
single  highly  efficient  specialist. 

Energy  Selectivity  and  the  Autonomous  Inertial  Oscillations 

The  advantages  for  an  autonomous  inertial  oscillator  can  be  argued 
from  the  point  that  it  may  provide  a  more  diverse  and   health   system. 
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The  ability  of  an  individual  system  to  develop  its  own  autocata I ytic 
Internal  rhythmic  guide  as  a  result  of  its  inertial  components  has  been 
shown  to  be  important  to  the  system  as  well  as  to  the  individual.  If 
all  systems  simply  internalize  external  environmental  signals  then  their 
would  be  less  selectivity  and  less  system  diversity.  The  ability  to 
generate  an  energy-selective  activity  aids  in  the  creation  of  more 
diverse  ecosystems  because  energy  rejected  by  on  component  of  the 
ecosystem  is  energy  available  for  use  by  another  component  of  an 
ecosystem.  While  diversity  is  generated  without  energy-selective 
behavior  the  addition  of  energy-selective  behavior  may  aid  in  the 
production  of  greater  diversity. 

When  individual  systems  produce  different,  yet  close,  internal 
autonomous  rhythms  the  interaction  of  their  internal  rhythms  with  that 
of  the  environment's  rhythms  produces  a  greater  variety  of  systems.  As 
has  been  stated,  the  simple  Internalization  of  external  variables  cannot 
produce  energy-selective  behavior.  Only  the  Interaction  of  an  internal 
autonomous  oscillation  operating  within  an  environment  subject  to  cyclic 
change  can  produce  an  energy-selective  activity  within  Its  system.  That 
Is  not  to  say  that  a  system  that  Internalizes  an  external  signal  could 
not  time  the  environment  but  It  would  not  have  the  ability  to  be  energy- 
selective.  Energy  selective  behavior  requires  the  system  to  produce  Its 
own  autonomous  oscillation!  The  oscillations  studied  for  this  research 
were  Inertial ly  generated  but  the  pr Inc  ip I e  may  be  expanded  to  include 
other  Internally  generated  osclllatons. 
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Energy  Selectivity  and  the  System 

What  does  the  system  obtain  by  selecting  for  system  components 
which  exhibit  energy-selective  behavior?  This  question  is  best  answered 
by  the  "principle  of  maximum  power". 

The  principle  of  maximum  power  is  stated  as  follows.  Systems  which 
obtain  the  most  available  energy,  use  that  energy  most  efficiently,  and 
provide  the  next  higher  system  with  a  means  of  maximizing  its  power 
ultimately  survive.  Energy  selective  behavior  by  individual  subsystems 
result  in  the  maximization  of  system  power  because  energy  selectivity 
produces  greater  diversity.  When  individuals  have  the  capability  of 
selecting  their  energy  Input  the  reciprocal  also  holds,  they  must  as  a 
consequence  of  their  selection  of  some  energies  reject  other  energies. 
The  rejection  of  energy  by  one  system  component  results  in  the 
availablity  of  that  energy  for  use  by  another  system  component. 
Therefore,  energy  rejection  by  one  system  component  releases  energy  for 
use  by  another  component  and  produces  greater  system  diversity. 
Increased  diversity  leads  to  greater  specialization  of  system 
components,  and  energy  specialization  through  the  use  of  energy 
selectivity  results  In  a  more  efficient  utilization  of  the  available 
energy  input  producing  greater  power  transfer  to  system  components 
higher  in  the  system  web  or  chain. 

The  simulation  plots  in  the  results  section  Indicate  that  energy- 
selective  specialization  requires  that  a  system's  energy-acceptance 
niche  be  closer  to  its  primary  energy  frequency  than  if  the  component  is 
a  general  1st.  Power  transfer  Is  for  the  most  part  more  efficiently 
accomplished  by  a  specialist  than  it  is  for  a  general ist.  The  exception 
occurs   when   a  general ist's  internal  rhythm  is  closely  matched  to   its 
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primary  energy  frequency  and  its  selectivity  range  includes  one  or  more 
secondary  energy  frequencies.  Therefore,  energy-selective  activity 
within  an  ecosystem  results  in  more  diversity,  more  efficient  use  of 
available  energy,  and  more  energy  and  power  transmitted  to  components 
higher  in  the  system  web  or  chain. 

Energy  Selectivfty  and  the  Energy  Acceptance  Niche 

The  development  of  an  energy-selective  behavior  contributes 
directly  to  the  establishment  of  a  system  niche.  When  the  free-running 
period  of  the  system's  Internal  rhythm  is  equal  to  and  phased  with  an 
environmental  cycle  the  system  obtains  the  greatest  amount  of  energy 
from  the  environment  and  as  a  consquence  generates  the  greatest  amount 
of  blomass.  If  the  free-running  rhythm  of  a  system  is  not  equal  to  and 
phased  with  Its  primary  environmental  cycle  the  energy  it  accepts  is 
decreased  or  entirely  rejected  and  its  blomass  production  decrease  or 
halts.  The  interaction  of  an  environmental  rhythm  with  that  of  the 
system's  internal  autonomous  rhythm  produces  selective  energy  aquisition 
within  the  system's  energy-acceptance  niche. 

When  a  system  component  fills  a  niche  it  occupies  a  space  in  the 
web  of  the  environment.  The  filling  of  a  energy-acceptance  niche 
signifies  that  the  component  occupying  that  niche  for  both  internal  and 
external  reasons  passes  the  system's  requirement  for  that  space  and 
function . 

When  a  system  possesses  the  ability  to  selectively  accept  or  reject 
energy  it  defines  a  range  of  inputs  It  will  compete  to  obtain.  The 
simple  general ist  will  try  to  obtain  energy  from  a  wider  range  of  inputs 
than  will  a  specialist,  and  when  no  other  system  or  component  offers 
competition   a  simple  general Ist  will  accept  energies  from  a  fairly  wide 
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range  of  alternatives.   However,   a  specialist   is  more  successful   in 

competition   within  its  defined  range  than  is  a  simple  general ist. 

The  complex  general  ist,   on   the  other   hand,   is   an   extremely 

competitive  system  component  and  can  compete  effectively  because  it  has 

developed  the  ability  to  time  its  environment  and  accept  energies  from  a 

variety  of  external  cycles.   The  ability  to  time  a  multitude  of  external 

rhythms   provides   the  complex  general  ist  with  a  competitive  advantage 

when  faced  with  another  competing  system  or  component.   Simulations  have 

shown   that  even  systems  with  the  ability  to  narrow  their   acceptance 

range  to  the  extreme  fail  in  their  attempt  to  outcompete  a  complex 

general ist   if  that  general ist  possesse  the  ability  to  time  the   primary 

rhythm  of  the  specialist.   The  complex  general ist  understands  more  about 

its  environment  than  other   systems,   however   it  can  outcompete   a 

specialist   only   if   it  has  the  ability  to  time  the   specialist's 

fundamental  rhythm.   A  specialist  and  a  complex  general ists  coexist  only 

if  the  complex  general Ist  does  not  possess  the  ability  to  time  the 

specialist's  fundamental  rhythm.   The  coexistence  of  systems,   at  least 

from  an   energy-selective  perspective,   occurs   from   an   ability   to 

selectively  obtain  different  environmental  energies  or  the  same  energy 

at  periods  close  to  both  competitors  fundamental  rhythms. 

Energy  Selectivfty  and  Diversity 

When  energy  input  is  available  the  ecosystem  finds  a  way  to  utilize 
that  energy!  This  is  a  partial  restatement  of  the  principle  of  maximum 
power.  The  utilization  of  available  energies  would  occur  even  if 
energy-selective  behavior  did  not  exist.  However,  for  the  same  input 
energy  available,  energy  selective  behavior  results  In  a  more  efficient 
and   specific   use  of  energy  by  the  components  of   a   system.   Energy 
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selective  system  components   acquire  input  energy  within   a  range  of 

alternatives.  Therefore  selective  behavior  allocates  energy  amoung 
system  components  and  amoung  systems,  and  as  a  result  of  their 
selection  of  energy  about  specific  primary  frequencies. 

When  available  input  energy  is  present  at  a  number  of  frequencies 
energy-selective  system  components  acquire  only  the  energies  available 
within  their  energy-acceptance  niches,  the  remainder  of  frequencies  are 
rejected.  The  rejected  energies  are  then  available  for  use  by  other 
energy-selective  system  components.  This  type  of  energy  allocation 
(between  a  variety  of  energy-selective  components)  increases  the  overall 
diversity  of  the  system. 

The  complication  of  energy  distribution  because  of  selective 
behavior  Is  indicated  by  examining  a  theoretical  system  with  only  one 
energy  input  consisting  of  two  primary  frequencies.  When  specialists 
are  available  to  accept  energy  at  both  frequencies  then  the  input  energy 
will  be  divided  by  the  two  specialists.  When  the  internal 
characteristics  of  the  two  specialists  are  the  same  and  only  their 
internal  rhythms  differ  then  more  energy  is  acquired  by  the  specialist 
with  the  nearest  rhythmic  match  to  that  of  its  specific  source  energy 
rhythm.  When  both  units  match  exactly  to  their  primary  energy  rhythm 
then  the  system  component  with  the  faster  Internal  rhythm  acquires  more 
energy,  but  both  units  survive. 

When  two  specialist  with  different  internal  characteristics  are 
available  to  accept  energy  from  the  source — competitive  exclusion  may 
result,  If  and  only  if,  both  primary  rhythms  of  the  source  are  also 
within  one  specialist's  energy  acceptance  niche.  When  energies  from  the 
two  primary  rhythms  of  the  source  are  allocated  between  the  specialists, 
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one  to  each  specialist  and  with  no  range  overlap  for  the  units,   then  no 

competition  occurs  and  the  two  specialists  coexist  on  the  same  energy 
source.  Moreover,  each  unit  receives  energy  from  one  of  the  two 
separate  rhythmic  components  of  the  same  source  energy.  Energy- 
selective  behavior  allows  system  components  to  exist  on  different  energy 
frequencies  form  the  same  source.  The  concept  might  be  thought  of  as  a 
means  of  allocating  energy  from  the  same  quality  source  amoung  groups 
of  different  systems.  In  other  words,  competition  need  not  occur  to 
allocate  different  periodic  energies  form  the  same  energy  quality 
source. 

Future  Investigation 

This  study  has  addressed  the  concept  of  energy  selectivity  from  an 
internal  autocata I yt ic  inertial  oscillator's  generation  of  intrinsic 
rhythmic  behavior  within  systems.  The  ability  of  systems  to  act  with  an 
inertial  backpressure  from  an  internal  oscillation  may  be  one  method  by 
which  energy-selective  designs  are  accomplished  for  general  ecosystem 
use.  However,  some  important  questions  still  remain  for  future 
investigation.  Is  it  possible  that  the  biological  clock  is  composed  of  a 
potential  and  kinetic  components,  and  if  so  do  these  components  vary 
widely  in  use  and  design?  Furthermore,  do  systems  have  other  ways  to 
utilizing  potential  and  kinetic  components  as  time  keeping  mechanisms? 

The  Autocatalytic  Inertial  Oscillator  and  Design  Characteristics 

Simple  electrical  and  mechanical  oscillators  are  constructed  of 
potential  and  kinetic  components  that  have  identifiable  properties  which 
allow  an  engineer  to  choose  which  component  best  fits  a  particular 
situation.   Ecosystem  oscillators  are  more  complex  and  the  eco legist  or 
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biologist  working  with  the  concept  of  energy-selective  ecosystem 
components  is  left  in  the  dark  with  no  idea  of  what  the  ecosystem's 
potential  or  kinetic  component  design  values  are.  While  it  appears  to  be 
relatively  easy  to  determine  the  free-running  period  for  a  particular 
system,  given  that  the  system  rhythm  generator  does  not  have  a  high 
depreciation,  no  data  are  available  to  determine  what  value  should  be 
given  to  the  potential  or  kinetic  ecosystem  components  which  may  be 
responsible  for  the  rhythm.  Therefore,  research  needs  to  be  conducted 
that  would  allow  the  determination  of  specific  values  for  ecosystem 
Inertial  or  rhythm  generating  components!  If  an  analyst  knew  that  a 
particular  component  acted  as  a  potential  generator  for  a  specific 
system's  rhythmic  behavior,  and  also  observed  the  free-running  rhythm  of 
the  system  it  might  provide  some  clues  as  to  the  nature  of  an  unknown 
kinetic  component.  New  research  should  be  directed  toward  identifying  a 
method  that  could  determine  the  potential  and  kinetic  properties  of 
components  for  individual  system  or  ecosystem  inertial  rhythms.  Once 
standard  values  for  Inertial  or  rhythmic  components  are  identified 
analysis  of  rhythmic  activity  might  be  more  easily  accomplished. 

Energy  Selectivity 

Energy-selective  behavior  in  general  ecosystem  designs  have  been 
clearly  indicated  by  the  plots  and  figures  presented  for  this  study. 
Much  work  has  already  been  completed  but  more  work  remains.  The  ability 
of  energy-selective  system  to  define  an  energy-acceptance  niche  has  been 
shown  quite  convincingly,  however  more  work  must  be  completed  to 
identify  how  much  of  an  effect  energy  selectivity  has  on  the  other 
aspects  of  niche  development. 
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The  Effects  of  Energy  Selectivity  Within  Autocatal ytic  Systems 

The  addition  of  autocata I ytic  behavior  within  energy  circuit  models 
clearly  indicates  the  advantage  of  internal  feedback  In  ecosystem 
design  (Odum  1976  and  1983).  If  that  system  also  utilized  autocata I ytic 
feedback  as  a  pumping  action  within  its  defined  energy-acceptance  range 
it  should  be  an  even  more  competitive  system  component.  Therefore, 
research  should  be  conducted  Into  the  effects  of  energy-selective 
activity  for  autocatal ytic  ecosystem  components. 

Energy  Selectivity  and  Competitive  Strength  In  Constant  Conditions 

Experiments  to  determine  the  competitiveness  of  systems  in  constant 
conditions  might  provide  Insight  Into  the  competitive  strength  of  the 
system,  and  could  be  used  to  establish  ratios  for  energy-selective 
behavior  by  type  of  system.  One  might  expect  that  systems  would  compete 
equally  as  well  In  constant  conditions  but  that  might  not  be  the  case. 
If  a  system  relied  heavily  on  its  rhythmic  Internal  guide  and  evolved  a 
highly  efficient  energy-selective  behavior  it  might  not  fair  well  In  a 
competitive  situation  with  another  energy-selective  or  non-energy- 
selective  system  when  both  were  placed  Into  constant  conditions.  The 
likelyhood  exists  that  energy-selective  behavior  is  stronger  than  non- 
energy-selective  behavior  only  when  the  system  Is  dependent  upon  a 
source  of  cyclic  energy.  The  unit  with  a  strong  energy-selective 
behavior  might  be  expected  to  compete  poorly  In  constant  conditions. 
Therefore,  an  energy  selectivity  index  based  on  competitive  strenght  in 
constant  conditions  might  be  employed  to  define  a  selectivity  Index,  or 
to  determine  if  the  system  possess  any  energy-selective  behavior. 
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Energy  Selectivity  and  the  Effects  of  High  Quality  Pulsing  Feedback 

The  effects  of  high  quality  pulsed  feedback  for  non-selective 
systems  is  well  defined.  Alexander  indicated  that  pulsed  feedback 
increases  the  power  of  a  system  because  the  recyc I e-to-rebu i I d  process 
uses  more  energy  and  power.  Therefore,  pulsed  feedback  in  energy- 
selective  systems  would  be  expected  to  increase  the  power  of  the  system, 
however  this  remains  for  future  research. 

Possible  Relevance  of  Autocatal ytic  Oscfllators  to  Physiological 
Clocks  within  Organisms 

Physiological  clocks  have  been  studied  for  approximately  2000  years 

with   the  earliest  Identification  of  a  biological  rhythm  occuring  2400 

years  ago.   Before  the  birth  of  Aristotle,   an  amateur  naturalist 

observed   that   legumes  had  a  rhythmic  pattern  to  their   leaf  movements 

(Brown,   Hastings,  and  Palmer,  1970).   By  the  beginning  of  the  twentieth 

century,   studies  of  rhythms  in  plants  and  animals  had  brought  about  a 

realization   that  these   internal  rhythms  were  the   result  of  complex 

processes  within  the  biological  unit.  Today  the  operation  of  biological 

clocks   is  still  a  mystery  and  internal  rhythms  have  been  Identified   in 

plants,   terrestrial  and  marine  organisms,   amphibians,  insects,  and  man 

(Brown,  1960;  Brown,  Bennett,  and  Webb,  1954;  Bennett  and  Stanely,  1960; 

Alder,  1970;  Goss,  Dinsmore,  Grimes,  and  Rosen,  1974).   The  major  points 

of  controversy  about  the   internal   operation  of   biological   clocks 

Include:    a)  whether  the  clock  Is  Internally  or  externally  timed  and  b) 

whether  the  clock  is  temperature  compensated. 

Internal  vs  external  timing  relates  to  the  ability  of  the  clock   to 

self-generate  a   free-running   biological  rhythm  or    whether   the  clock 

merely  reflects  the  signals  of  its  environment.   The  simulation  models 


Figure  66:  Energy  diagram  of  a  temperature  compensated  autocatal yti c 
Inertlal  oscillator.  The  clock  appears  to  be  capable  of 
explaining  the  observed  temperature  compensation  by 
c locks. 
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completed   for   this  study  indicate  that  an  internal  oscillation  with   a 

free-runnung  rhythm  is  necessary  for  the  generation  of   energy-selective 

system   behavior.    If   biological   clocks  were   autocatal yt ic   inertial 

oscillators  internal  timing  would  be  advantagous  tor  both  the  biological 

unit   possessing   the  oscillator  and  for  the  system   within   which   that 

component  functions. 

Temperature  compensation  of  biological  clocks  has  been  propossed  to 

explain   the  observed  period  stability  of  biological  clocks  for   varying 

temperatures.   The  effects  of  temperature  on  the  period  of  a  biological 

clock   is  quite  interesting.    Research  indicates  that  biological  clocks 

are  relatively  insensitive  to  temperature  changes  in  their   environment, 

producing   relatively  stable  periods  for  temperatures  within  a  range  of 

o       o 
5  C  to  30  C  (Bennett,  1974)  . 

Three  different  hypotheses  have  been  proposed  to  explain  the 
accuracy  of  the  biological  clock  with  changing  temperature  (Palmer, 
1976).  First,  the  biological  clock  is  constructed  around  a  diffusion 
process  and  diffusion  rates  are  less  sensitive  to  temperature  changes. 
Second,  the  chemical  reactions  of  the  clock  are  temperature  compensated. 
Third,  the  biological  clock  is  externally  timed  and  therefore  is  only 
processing  a  rhythm  and  does  not  need  to  temperature  compensate. 

If  biological  rhythms  were  generated  by  an  autocata I ytic  inertial 
oscillator  temperature  compensation  could  be  accomplished  by  the 
interaction  of  the  flows  within  the  oscillator.  The  model  in  Figure  66 
(although  not  simulated)  would  account  for  the  effects  of  stable  period 
development  with  changing  temperature.  The  flows  all  change  the  same 
amount  with  increased  or  decreased  temperature  and  as  a  result  the 
oscillator  period  remains  stable. 
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The  advantages  for  energy-selective  activity  or  behavior  have 
clearly  been  demonstrated  by  the  simulations  of  this  research.  The 
ability  to  act  in  an  energy-selective  manner  has  been  shown  to  occur 
when  a  system.  In  a  cyclic  environment,  deveiopes  an  internal  free- 
running  rhythm.  The  energy-selective  system  selectively  aquires  input 
energy  more  efficiently  and  transfers  more  power  with  greater  efficiency 
to  the  remainder  of  the  system.  The  ability  to  aquire  input  energy  more 
efficiently  and  to  use  that  energy  more  efficiently,  as  stated  by  the 
principle  of  maximum  power,  provides  the  system  with  a  greater 
competitive  edge.  Therefore,  the  development  of  biological  oscillators 
or  rhythm  generators  should  be  advantageous  to  the  individual  biological 
unit  as  well  as  the  system  it  functions  within. 
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XS=R5IN/( JRC« 10 ) 

X6=R6IN/( JRC* 10) 

XTsRTIN/C  JRC«I  01 

xa=RaiN/io 

X?=GMEaA**2 

S1C=R9IN/X23 

TY=36S 

T  C=  1  .  0 

CMEGAY=2.0*3 . 1*1592  7/TY 

CMEGAO  =  2  .0*3  .  L -1.1  i  ^J  7 /TO 


Ovr.AM  [; 
NC3CWT 


JGfkT 


3i_LAR=  J  :  ,N  ( i.'iE  J-?*  T  :  •*=  ) 

IF     (S.JLA^     .LT.      3.0)      3CLAP=0.0 

jf^  =  sj.N/  (  1 .0  +  x  :»NUT^cLacK) 

Su;.i=  <  ■-»759-K0365*S  I  N  (  CMECAY  «T  tME)  )-«SCLAR 


JRl=SUNy{ 1.0+X5*ABS(Q) ) 
PU=INTGRL(  IPU, FO  ) 
FO=XO«CLCCK*NUT* JR 
PL=INTGRL(  IPL.LCAO) 

LCA0=x2*eic-t-x4»e  lo 

XPRG=  JR  1  «Aas  <  C  ) 

QoaT=x7*xppo-x6*  xPRa-xg«s-x6«a 
a=iNTGriL(  :q.  gcct  j 

S=INTGRL(S1C,G) 
F<3=xg«S 

xpaT=(a*»2)/ (2*c) 

XKlN=(L/2  )  *F9-*»2 

CLCCK=XPaT«-XK  IN 

aDaT  =  X  I  *JR*NUT»CLCC.<-X2«-3  10 

3rG=INTGRL(BICC,aOCT) 

ND0T  =  X2»Q  ra-Xj«JP«NUT«CLCCX-X4a,NUT 

^■-T=  I.iTGflLtMUTC.NOOT) 

term: n-l 

;]^^5-^FINTr.1=3e5  0.0.0ELMIN=2.0E-4  5,0£LT=0.0  1,POC£L  =  10 

^^\^\^    i :  :  .NUT  ,pu  ,PL 

5r:p 

CNJJC3 

/« 
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•/MARSH    JOS     (2002.2002.90.2.0) , PAUL. Z^ICX .CLASS^l 

//STEPl        EX£C       CSMP. 

//  CCMP=, 

//   .  l.K£D= 

INITIAL 

I NCQN    BiaC=523.0 

INCQN    0£AOC=5226.0 

INCON    NUTC=522.0 

INCON    PUI=0.0 

INCON    PH=0.0 

PARAMETER    JC=3450.0 

PARAMETER    fl0IN=2T60.0 

PARAMETER  JRC=6gO«0 

PARAMETER  R1IN=4.0 

PARAMETER  R21N=2.l 

PARAMETER  R3IN=0.2l 

PARAMETER  R4.IN=1,59 

PARAMETER  R5IN=0.4a 
•  PARAMETER  R6rN=0.29 

PARAMETER  R7IN=0.09 

PARAMETER  RaiN=0.39 

PARAMETER  R9IN=0.73 

PARAMETER    RIV=i5.3gE5 

PARAMETER    TI0=1.13Ee 

X0=R0IN/  (JRC*NIJTC) 

X1=R1IN/C JRC*NUTC) 

X2=fl2IN/BiaC 

X3=fl3INy'2I0C 

X^sRo.lN/'OEADC 

X5=flSlN/(DEAOC«TID) 

X6=fl6lN/R4IN 

X7=RT I N/ ( NUT  C«  T  r  0 ) 

xa=RaiN/Riv 

X9=Hg  rN/( JRC'NUTC ) 
aMeGAY=2.0»J  .14.15/365 
0MEGAT=2.0*3.  14.15/14.77 
DYNAMIC 

RIVER=4.26E5>2.12ES»CGS (OMEGA Y« T I  me ) 
TIDE=1.13E6+5.65E5*SIN( CMEGAT^T IME) 
JR=SUN/{  1.0>X0*NUT) 
SUN=3450>22S  0*SIN(OMEGAY-«TIME) 
PU=INTGRL(PUI .PLCXS) 
PU=  I NTGRl.(P1_X  ,FLCV«l.J 
FLCWS=XO*JR«NUT 

Fuaiii.=(  (X2*aici  ♦  (X3»aia)  ) 
aooTsxi  «jR«NuT-x2*3rc-x3«a  ic 
aic=iNTGRL(a loc.aooT) 
ooaT=x2»a ic-x4-»oEA0-xs-«cEAC*r  ide 

OEAO=INTGRU£  DEAOC.OCOT) 

NOaT=xa«Rl VER>X3*SIC*X6*X4*0EA0-X7«NUT*TIDE-Xg*NUT« JR 

NUT=INTGRLCNUTC .NOCT) 
TERMINAL 

TIMER    FINTlM=3e50  .  0  ,0£l_-'4 1  N=2  .  CE-4  5  .  OEl.T=0  .  l.?ROEL=l  0.  0 

METHCO     SIMP 

PRINT  3  10  .C£AO  .NUT.PU.PL 
END 
STOP 

ENOjoa 
/* 
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//XARSH2  JOS  C2002,2002.90,2,0)  .PAUL.  2*ICKtCl.ASS-l 

//STEP  I  £;<EC   CSMP, 
//  CCMP=, 

//  LKHD= 


INXTIAl. 


DYNAMIC 
NCSQRT 

SORT 


INC  ON    BI0C=S23.0 
INCON    OeAOC=S226.0 
INCON    NUTC=523.0 
INCON    PUI=0.0 
INCON    P1_I=0«0 
PAfiAMETEH     JC=3450.0 
PAPAMETER     JRC=690.0 
PARAMETER    RO  IN=2760«0 
PARAMETER     H1IN=A.0 
PARAMETER    R2IN=2.1 
PARAMETER     R3IN=0.21 
PARAMETER    RalN=l«59 
PARAMETER     R5IN=0.4.a 
PARAMETER    R6IN=0.29 
PARAMETER    R7IN=0.09 
PARAMETER     RarN=0.39 
PARAMETER    R9iN=0.73 
PARAMETER    RIV=0«39E5 
PARAMETER    TT0=1.13E6 
XO=flOIN/(  JRC«NUTC) 
Xl=RXIN/( JRC«NUTC) 
X2=R2IN/8iaC 
X3=R3IN/S  IOC 
X<V=flAIN/OEAOC 
X5=HSIN/(DEA0C*TI0) 
X6=R6tN/R4.IN 
Xr=R7IN/ {NUTC«TIO) 

xa=RaiN/Riv 

X9=R9lN/(  JRC*NtJTC) 
aMEGAY=2.0«3 .1415/3  65 
QMEGAT  =  2.0*3  .14.i5/;>*«77 
aMEGAO=2.0*3.  14.1  5 

RIVER=A..26E5-t-2.X3E5«C0S(aM£GAY«TIME) 
TIOE=l.  13E6>5.65E5*SI.N(aMEGAT«T  IME) 

JR=SUN/ (  1  .0>XO»NUT) 

SCLAR=SIN( CM£GAO*T I  ME ) 

IF     (  SOLAR     .i_T.     0.0)     SCI-AR  =  0.0 

SUN=( 9  759+63  6  5*5 IN(CMEGAY»T IME) )* SOLAR 


PU= INTGRL<PUI .FLOWS) 
PL=  INTGRL(PLI  .i=LawL  ) 

FLawS^XO«JR*NUT 

FLC»L=( (X2*aiC)+(X3«3iaj ) 

800T=X1  •JR»NUT-X2*9 IQ-X3«aiC 

aX0=INTGRL(3iaC.aD0T ) 

0  0aT=X2«aia-X4*OEAO-XS*DEAC*T lOE 

DEAO= INTGRLI OEAOC.OOOT) 

NOQT=X3*RIVER+X3*Bia*X6*X4*0EAO-X7*NUT«TT0E-X9*NUT*JP 

NUT=INTGRL(NUTC.NCCT) 

TIMER  FINTIM=365  0.0.0ELMIN=2.  OE-4  5  .  DEL  T=0  .  0  I  .P«^CEL=l  0  .  0 

METHOD  SIMP 

PRINT  aiQ.OEAO.NUT .PU .PL 
END 
STOP 
ENOJOB 
/* 
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/•/MARSH5     JCS     (2002.2002  .  150  ,2iO  ),  PALI.  .ZmCK  fCLASS=l 

//STEPl      excC      cs.-^p. 

//  CC."^P=. 

//  LKED= 

INITIAL 

INCCN     IQl=2.£150 

INCCN     IQ2=2.6150 

INCQ.N    3  ICCl=2cl  .50 

INCCN    5ICC2=2ei.50 

INCCN    JcACC=522g.O 

INCCN    riUTC=5  23.0 

INCCN     IPLII=0.0 
INCCN    IPU2=0.0 

INCCN    IPL1=0.0 

INCCN     IPU2=0.0 

PAPAMtTES     JC=3450.0 

PAPAMETER     JRC=690..J 

PARAMETEP    P0IN=l3a0.0 

PARAMETER     R1IN=2.0 

PARAMETER     R2IN=1.0S 

PARA-»ETER     R3IN=0.105 

PARAMETER    R4.INil.59 

PARAMETER     R3IN=0.4.3 

PARAMETER     R6IN=0,29 

PAPAMETER     R7IN=0.09 

PARAMETER     R8IN  =  0.3'5 

PARAMETER     R9lN=0,3d5 

PARAMETER     R10IN=1330.3 

PARAMETER     .^lirN  =  2.0 

PARAMETER    R12IN=1.0S 

PARA.HETES     =13IN=0.105 

PARA.METER     R19IN-0.3c5 

PARAMETER     R20  IN-=133G.O 

PA?.AMET£H     «2  1[N  =  2.69 

PARAMETER    R22rN=l.05 

PARAMETER     R23:n=0.105 

PARA/'iE'^ER    R2*IN=0.69 

PARAMETER     R2SIN=l3rt0.0 

PARAMETER  R26lN=2.o9 

PARAMETER  R27IN=l,05 

3APAMETER  R23IN-0.105 

PARAMETER     R29I,N=0.c9 

PAfiAMETE.=     RIV-c.3'5£5 

P  AR  AMET  £R     T  IJ= I  .  1 23  i 

PARAMETER     Z~i*:z.-^>- 

FREQl=l .00 

F  =  E02=0  .00274. 

GM£GAl  =  2»J.  l*15927-sFKC01 

aM£GA2=  2*3.1  AL3-;27*FREa2 

1_1=  l/(  C'C.Md  JA  l:»»2  ) 

L2=  l/(C'>C.ME3A2»»2  ) 

CLK1  =  (  (  I  Qi«*2)/  (2*C  J  )-f  (  (L  1/2)  »  (  =^23  I  N  *  «2  )  ) 

CLK2=(  (  IQ2*»2  )/(  2«C)  )  +  (  (L2/2)  »  ( .'^2  3  I  N«  »2  )  ) 

XC=fi0  I.M/  (  JRC*NUTC«C1_.<  I  ) 

X  l=fil  IN/  (  JRC»NUTC*C1.)<1  ) 

X2=R2L^/S10C 1 

x3=n3iN/ciac I 

X-^=fi4.  I  M/CEAOC 

X==P5IN/( OEACC'TIO ) 

Xe=R6  I  -J/BA  IN 

X7=R7  I.J/  (NUTC*T  :0  ) 

X  a=  R  3  I  .N  /  R  I  V 

X';iR9X;j/(jRC«NUTC*CL/<l) 

X10=R1.)  IN/(  JKC'=NUTC*a./<2) 

X  I  1=R1  I  IN/  (  Jr;C*NUTC*Cl_;<2) 

X  12  =  P 12IN/a  ICC2 

X 12=R IJ IN/5 IJC2 

X  1-3=R1')  IN/(  JRC*NUTC*CL.<2) 

X20=P20lN/(jPC*iai) 

X2  1=R21 rN/( JRC'IOl ) 

X22=r<22  IN/I  31 

x2J=C:-\EjA  I  '<«2 

X2  4.=fi2i  IN/  ( j^;c*  :oi  ) 
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S  l.C=R23IN/X23 

X25=R2SIN/( JRC*I G23 

X26=R2  6 IN/ (jaC» ta2J 

X27=R27IN/I 22 

X2a=0.MEGA2^»2 

X2'9  =  R2='  IN/  C  JBC»  102  ) 

S2C=R2a IN/X23 

TM=  14.. 77 

TC=1  .0 

GMeGAY=2.0»3  .1 4  15  92T/TY 

0«EGAM=2.0*3«  L^15<527/T.M 

OMEGAO=2.0*3-.l«  15927/ TC 

DYNAMIC     pj^£j,^^^26£5-f2.l3e5«CQS(aMEGAY*TIM£) 
TID£=l.l3Ee+5.6S£5-»2IN(CMEGA,vi»TIME) 

NOSORT 

SCL-AR=SIN(aMEGAO*T  I  ME) 

IF      CSGLAR     .LT.     3.0)      SCUAR=0.3  ^^^-.v 

JR  =  5UN/  (  1  .0  +  XC»NUT*Cl_CCi<l.-t-X  L0*NUT»CL.CC.<2) 
SUN=(9  75g>63  65*SIN(CMECiAY»T  IME)  )*SC1-AH 

SORT 

JP  l=5UN/(  1.0-fX20=»A82(ai)) 

JS2=SUN/( 1.0>X25«Ac3(Q2) ) 

XPRC  l  =  J R l»Aa3 ( CI  ) 

XFRC2=JR2*AdS( 02) 

Q0CT1=X21  «XPP CI  -x2*^X?fvCl-X 23*5  1 

a  1=INTGRL(  IQl  .GCCTl  ) 

S  1=  INTGRLC  SlC.a  I  ) 

FLCV»23=X23*S1  ^,     ^^ 

QCCT2=X26»XPRC2-X29*XPRC2-X23«S2 

Q 2= INTGRLC I0  2.GaCT2) 

S2=INTGR1.(  3  2C.02  ) 

FLCl*2a=X23«S2 

XPQT  l=(Ql*«2)/(2*C) 

XKINI=(L1/2)*FLCW23*«2 

CLCCKl  =  X?CTl-fXKIM 

XPCT2=(   ■iZ'^-^Z)  /  (  2»C} 

X<IN2=  (  L2/2  )  *F^L.-'2°-'*2 

CL'-C.<2=XPCT2>XK  [N2 

3wCTl  =  X  1«JR=»NUT=CI_CC:<L-X2  =  ^I-:  1-X3»-J  IM 

diJi=rJT^r^Lcaicci,3caTi  ) 

3  0CT2=X1 1 «jH-\Lr«CLCCK2-X I2«d I32-X13=3:a2 
0  IC12=I.MrC;-U(^'I  CC2,33Lr3)  

0L;GT  =  X2»d  lo  I  *X  1  2  »d  I  3  2- X  4  aCCALi-X  5*  OeAD*TI;^c 

jeAO=irjTCRL(  C-I^CC  .  CCoT  )  _      ,        ,     .  .  , -r  .  t  r  - - 

NDCT  =  X8*ft  IVE?>XJ*'3ICl^-XlJ«3t32*-X3«X.4*0i:A0-X7aNjT*TI^w 

-XQ^NUT  "JRrtCl-CCK  1-XI'9«JR*NUT«CLGC)<2 

NLT=  IMTGRI_(NUTC  .NOCT) 

FLCVmO— X0«Jf5»NUT*CLL:O<l 

FLC*  10=  XI  jajR»NUT»CLCC.<2 

FLCV.1_I  =  X2*'JI  CI  >X3*3I0  1 

FLL;'*l..:-=Xl2«3IC2-»-Xl3*d:02 

PU1=  I.nTGSLC  I°U  1  .FLO  *0  ) 

PU2=  I. NT  GPL  (  [PU2  iFl.C'.*lQ  ) 

?LI=INTGPI.(IPLI,FL0*L1) 

PI-2=  IN  rGPL(  IPL2.rUC*i-2  ) 

T£SMINAJ__^^^^^^  F  INTIM=3  65.:;.0  .JE1_MIN=2  .  OE-A  5  .  C^IL  r  =  0  .  0  I  ,  P"  C£L=  I  0  .  0 
METHOD  SIMP 
Pr.INT  U  ICl  .3  IC2.3EAC.NUT 

ENO 
ST3P 
END J  03 
/» 
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